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Part A 

Classical Solutions in Einstein’s Gravity 



Chapter 1 


General Introduction 

1.1 Cosmological Implications of the Theory of varying Speed 
of Light 

The constants of nature are not constant at all. During the early phase of the universe, 
the constants of nature were varying signihcantly, although now-a-days this variation is 
quite insignificant. Assuming their constancy at all times requires massive extrapola¬ 
tion, with no observational basis. Could the universe come into being riding the back 
of wildly varying constants? Several constants of nature have been stripped off their 
status in theories proposed in the past. Physicists have long entertained the possibil¬ 
ity of a varying gravitational constant G (e.g., Brans-Dicke theories [Brans and Dicke, 
1961]), a varying electron charge e [Bekenstein et al, 1982, 2002] and more generally 
varying coupling constants. Indeed with the advent of string theory, the variation of 
these constants seems to be fashionable. In sharp contrast, the constancy of the speed of 
light has remain scared and the term heresy is occasionally used in relation to Varying 
Speed of Light Theories (i.e., VSL Theories) [Barrow, 1999]. The reason is clear: the 
constancy of c, unlike the constancy of G or e, is the pillar of special relativity and 
thus of modern physics. Varying c theories are expected to cause much more structural 
damage to physics formalism than other varying constant theories. 

The hrst varying constant was the speed of light as suggested by Kelvin and Tait [1874], 
which was thirty years before Einstein’s proposal of special relativity. Thereby a varying 
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c did not shok anyone, as indeed c, nnlike G played no special role in the formalism of 
physics. This is to be contrasted with the state of affairs after 1905, when Eddington 
commented “A variation in c is self-contradictori/' [Eddington, 1946]. This astonishing 
statement does a disservice to the experimental testability and scientihc respectability 
of the theory of relativity. The Eddington’s statement has now been transmuted into 
“asking whether c has varied over cosmic history is like asking whether the numbers of 
liters to the gallon has varied ” [Duff, 2002]. The implication is that the constancy of the 
speed of light is a logical necessity, a dehnition that could not have been otherwise. For 
centuries the constancy of the speed of light played no role in physics and presumably 
physics did not start being logically consistent in 1905. Furthermore, the postulate of 
the constancy of c in special relativity was prompted by experiments (including those 
leading to Maxwell’s theory) rather than issues of consistency. Seen from an angle, even 
in a world where all seems to vary and nothing is constant, it is always possible to dehne 
units such that c remains a constant. So varying c is tautological and is tied to the 
dehnition of a system of units. 

In classical electromagnetism the speed of light is only constant in vacuum and it 
‘varies’ in dielectric media. The changing c breaks Lorentz invariance and conservation 
of energy. The geometry of Universe is not affected by a changing c. One can allowed a 
changing c to do the job normally done by superluminal expansion. For changing c, the 
gravitational laws should be modihed. The basic assumption is that a variable c does 
not induce corrections to curvature in the cosmological frame and therefore, Einstein’s 
equations, relating curvature to stress energy are still valid. The reason behind this 
postulate is that c changes in the Local Lorentzian frames associated with cosmological 
expansion. The effect is a special relativistic effect and not a gravitational effect. There¬ 
fore curvature should not be related with the variation of c. 


The good news for experimentalists is that once the theoretical choice is made, the 
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different theories typically lead to very different predictions. For instance, dilaton theo¬ 
ries violate the weak equivalence principle but maintain the Lorentz invariance whereas 
VSL theories maintain weak equivalence principle [Magueijo et al, 2002; Moffat, 2001] 
but break Lorentz invariance. These differences have clear observational implications. 
For instance the “STEP” satellite could soon rule out the dilaton theories capable of 
explaining the results given by Webb et al [http://einstein.stanford.edu/ STEP/]. Some 
VSL critics fail to notice that there are theories for which a varying c is actually a di¬ 
mensionless statement. This does not include the results of Davies et al [2002], where a 
dimensionless statement on varying c is achieved only because a constant G was assumed. 

Even after the proposal of special relativity in 1905 many varying speed of light theo¬ 
ries were considered, most notably by Einstein himself [1911]. VSL was then rediscovered 
and forgotten on several occasions. For instance, in the 1930’s VSL was used as an alter¬ 
native explanation for the cosmological redshift. This theory conflicts with hne structure 
observations [Stewart, 1931; Bue, 1932; Wold, 1935]. None of these efforts relates to re¬ 
cent VSL theories, which are hrmly entrenched in the success of the hot big bang theory 
of the universe. According to the concept, hrst modern VSL theory was Moffat’s ground 
breaking researches [Moffat, 1993], where spontaneous symmetry breaking of Lorentz 
symmetry leads to VSL and an elegant solution to the horizon problem. 

The basic dynamical postulate is that Einstein’s held equations are valid even when 
c 7^ 0, with minimal coupling (i.e., with c replaced by a held in the relevant equations) 
in this particular frame: 

Gfj,u — ( 1 . 1 ) 

This is inspired by the statement of Maxwell’s equations in dielectric media. The 
hrst strong assumption is that the held c does not contribute to the stress-energy tensor. 
More importantly, this postulate can only be true in one frame. Hence the action remains 
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as shown by Barrow et al [1998, 1999] 


5 






/ c\R + 2A) 
\ IGttG 



+ Rc 


( 1 . 2 ) 


The dynamical variables are the metric any matter held variables contained in 
the matter Lagrangian Cm and the held c itself. The Riemann tensor and Ricci scalar 
are to be compnted from at constant c in the nsual way. This can be trne in one 
frame. Therefore Einstein’s eqnations do not acquire new terms in the preferred frame. 
Minimal coupling at the level of Einstein’s equations is at the heart of the model’s ability 
to solve the cosmological problems. It requires of any action principle formulation that 
the contribution Cc must not contain the metric explicitly and so does not contribute 
to the energy-momentum tensor. 


Like inhation [Guth, 1981; Albretch et al, 1982; Linde, 1982, 83], modern VSL theories 
were motivated by the cosmological problems— the hatness, horizon, entropy, homogene¬ 
ity, isotropy and cosmological constant problems of big bang cosmology ([Albretch et 
al, 1999; Magueijo et al, 2000] for a review and [Avelino et al, 2003] for a dimensionless 
description). At its most basic VSL was inspired by the horizon problem. It does not 
take much to see that a larger speed of light in the early universe could ‘open up the 
horizons’ [Moffat, 1993; Albretch et al, 1999]. More mathematically, the comoving hori¬ 
zon is given by r/j = | which indicates that in the past Vh must have decreased so as to 
causally connect the large region at present. Thus 

d c 
d c 

that is, either one can get accelerated expansion (inflation) or a decreasing speed of 
light or a combination of both. This argument is far from general: a contraction period 
(a < 0, as in the bouncing universe) or a static start for the universe (a = 0) are contra¬ 
dictory examples. However the horizon problem is just a warm up for the other problems. 
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The VSL cosmological model first proposed by Albretch and Magueijo [1999] making 
use of a preferred frame, thereby violating the principle of relativity. In this respect a 
few remarks are in order. Physicists don’t like preferred frames but they often ignore 
the very obvious fact that they have a great candidate for a preferred frame: the cos¬ 
mological frame. Modern physicists invariably choose to formulate their laws without 
reference to this preferred frame. The VSL theory proposed by Albretch et al [1999], 
makes a radically different choice in this respect: it ties the formulation of the physical 
laws to the cosmological frame. Equation (1.1) can only be true in one frame; although 
the dynamics of Einstein’s gravity is preserved to a large extent, it is no longer a geo¬ 
metrical theory. However, if c does not vary by much, the effects of the preferred frame 
are negligible, of the order where n is the rank of the corresponding tensor in 

general relativity. Inertial forces in these VSL theories may turn out to be a strong 
experimental probe and are currently under investigation. 


Since Friedmann equations are still valid even when c 7 ^ 0. Thus the Friedmann 
equations are 


SttG 


kc^ 


P 


(1.4) 


and 


d 

a 



( 1 . 5 ) 


where pc^ and p are the energy density and pressure respectively, k = 0, ±1 is the spa¬ 
tial curvature and a dot denotes a derivative with respect to cosmological time. If the 
universe is flat {k = 0 ) and radiation dominated {p = |pc^), we have as usual a oc 
As expected from minimal coupling, the Friedmann equations remain valid, with c being 
replaced by a variable in much the same way that Maxwell equations in media may be 
obtain by simply replacing the dielectric constant of the vacuum by that of the medium. 
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Combining the two Friedmann equations (1.4) and (1.5) now leads to: 


„a / p\ 3kcc 


( 1 . 6 ) 


i.e., there is a source term in the energy conservation equation. This turns out to 
be a general feature of this VSL theory. It is to be noted that in general relativity 
stress-energy conservation results directly from Einstein’s equations, as an integrability 
condition via Bianchi identities. By trying this theory to a preferred frame, violations of 
Bianchi identities must occur [Besset et al, 2000] and furthermore the link between them 
and energy conservation is broken. Hence one expect violations of energy conservation 
and these are proportional to gradients of c. In the brane world realization of VSL 
[Youm, 2001], these violations of energy conservation are nothing but matter sticking to 
or falling off the brane. 


In the context of preset c{t), two scenarios were considered: phase transitions and 
Machian. In phase transitions, the speed of light varies abruptly at a critical temperature 
as in the models of Moffat [1993] and Albretch et al [1999]. This could be related 
to the spontaneous breaking of local Lorentz symmetry [Moffat, 1993]. Later Barrow 
considered scenarios in which the speed of light varies like a power of the expansion factor 
i.e., c cx a"^, the so-called Machian scenarios. Taken at face value such scenarios are 
inconsistent with experiment (Barrow et al [2000] for an example of late time constraint 
on n). Such variations must therefore be conhned to the very early universe and so 
the c-function considered by Barrow [1999] and Barrow et al [1999] should really be 
understood as 

where oq is the scale at which VSL switches off. The problem with putting in ‘by hand’ 
a function c{t) is that the predictive power of the theory is severely reduced [Moffat, 
2001 ]. 
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Now one may consider more quantitative regarding conditions for a solution to the 
cosmological problems. Now starting with the horizon problem, one can consider the 
phase transition scenario [Moffat, 1993; Albretch et al, 1999] for the time tc when the 
speed of light changes from c~ to . The past light cone intersects t = tc aX a. sphere 
with comoving radius r = c+(?7o ~ Vc) where rjQ and rjc are the conformal times now 
and at tc respectively. The horizon size at tc, on the other hand, has comoving radius 
Th = c~rjc- If ^ then r r^, meaning that the whole observable universe today 

has in fact always been in causal contact. This requires 


c~ 1 1 T+ 


( 1 . 8 ) 


where Zeq is the redshift at matter radiation equality and T+ and T+ are the universe 
and Planck temperature after the phase transition. If k, this implies light travel¬ 
ing more than 30 orders of magnitude faster before the phase transition. It is tempting, 
for symmetry reasons, simply to postulate that c~ = oo but this is not strictly necessary. 


Considering now the flatness problem, let pc be the critical density of the universe: 


Pc = 


( 1 . 9 ) 


SttG \a^ 

i.e., the mass density corresponding to the flat model {k = 0 ) for a given value of A 
Now dehning deviations from flatness in terms of e = — 1 with Q = the equations 

(1.4), (1.5) and (1.6) are simplihed to 


e = e(l + e)(2 + 37 )- + 2e - (1.10) 

a c 

where 7 — 1 = ^ is the equation of state (7 = 1, 4/3 for matter/radiation). In 
the standard big bang theory e grows like in the radiation era and like ‘a’ in the 
matter era, leading to a total growth by 32 orders of magnitude since the Planck epoch. 
The observational fact that e can be at most of order 1 nowadays requires either that 
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e = 0 strictly or that an amazing fine tuning must have existed in the initial conditions 
(e < 10“^^ at t = tp). This is the fiatness puzzle. 

As equation (1.10) shows, a decreasing speed of light < 0) would drive e to 0, 
achieving the required tuning. If the speed of light changes in a sharp phase transition 
with I - |S> then a decrease in c by more than 32 orders or magnitude would suitably 
flatten the universe. But this should be obvious even before doing any numeries, from 
inspection of the non-conservation equation (1.6). Indeed if p is above its critical value 
(as in the case for a closed universe with k = 1) then equation (1.6) tells that energy is 
destroyed. \i p < pc (as for an open model, for which k = —1) then energy is produced. 
Either way the energy density is pushed towards the critical value pc- In contrast to 
the big bang model, during a period with ^ < 0 only the fiat, critical universe is stable. 
This is the VSL solution to the fiatness problem. 

VSL cosmology has had further success in resolving other problems of big bang cosmol¬ 
ogy that are usually tackled by inflation. It solves the entropy, isotropy and homogeneity 
problems [Moffat, 1993; Albretch et al, 1999]. Also it solves at least one version of the 
cosmological constant problem [Albretch et al, 1999] ([Moffat, 2001] for a discussion of 
the quantum version of the Lambda problem). Further VSL cosmology has a quasi-fiat 
and quasi-Lambda attractor [Barrow et al, 1999] (i.e., an attractor with non-vanishing, 
but also non-dominating Lambda or curvature). 

These cosmological implications of VSL have led to much further work. The stability 
of the various solutions was demonstrated by Szydowski et al [2002] and Biesiada et al 
[2000] using dynamical systems methods [Cimento et al, 2001; Gopakumar et al, 2001]. 
The role of Lorentz symmetry breaking in the ability of these models to solve the cos¬ 
mological problems was discussed by Avelino and Martins [1999]. Also combinations of 
varying c and varying G have been studied by Barrow et al [Barrow et al, 1999]. Barrow 
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[2002] has cast doubts on the ability of this model to solve the isotropy problem if the 
universe starts very anisotropic. It was also found that if c falls fast enough to solve 
the flatness and horizon problems then the quantum wavelengths of massive particles 
states and the radii of primordial black holes will grow to exceed the scale of the particle 
horizon. However this statement depends crucially on how the Planck’s constant h scales 
with c in VSL theories. 

The relation between VSL and the second law of thermodynamics was investigated 
recently by Barrow [2002], Youm [2002], Coule [1999] and Cimento et al [2001]. In par¬ 
ticular Cimento et al [2001] found that if the second law of thermodynamics is to be 
retained in open universes, c can only decrease, whereas in flat and close models it must 
stay constant. A similar discussion in the context of black holes can be found in Davies 
et al [2002]. It is not surprising if the second law of thermodynamics were violated in 
VSL models with hard breaking of Lorentz invariance. After all the first law is violated 
in these models. However a deviation from hrst principles remains elusive. 

It is nevertheless interesting that several observational puzzles can be solved with VSL. 
With supplementary the redshift dependence in hne structure constant a could be seen 
as a result of a variation of c. Another puzzle was the observation of rare very high energy 
cosmic rays, in conflict with standard kinematic calculations based on special relativity 
which predict a cut-off well below the observed energies. This could present the hrst 
experimental mishap of special relativity and evidence for some VSL theories. Finally, 
the recently observed accelerating universe may be well explained with the variation of c. 


1.2 Quintessence 

For the last few decades, it is generally believed that after a small period of inhationary 
era, when the universe was accelerated, the universe has been evolving in a big bang 
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scenario (deceleration of the universe) and most of the present day observations are in 
accord with the Standard Cosmological Model (SCM). The SCM can give a satisfactory 
explanation to other observational properties of the present universe (e.g., primordial 
nucleosynthesis, extragalactic sources redshift, cosmic microwave radiation). But diffi¬ 
culty has started when data from high redshift type la Supernovae ( 2 ; ~ 1) observations 
[Perlmutter et al, 1998, 1999; Riess et al, 1998; Garnavich et al, 1998] suggest an accel¬ 
erated universe at the present epoch. So there must be some matter field which is either 
neglected or unknown responsible for this accelerated universe [Bachall et al, 1999]. This 
type of matter held is called quintessence matter, shortly Q-matter and the problem is 
called quintessence problem. According to Caldwell et al [1998] and Sahni et al [2000] 
this Q-matter can behave like a cosmological constant by combining positive energy 
density and negative pressure related by the equation of state p = —p [Turner, 1997; 
Ostriker et al, 1995]. 

However this does not constitute the only possibility: among all the other candidates, 
the missing energy can be associated to a dynamical time-dependent and spatially ho¬ 
mogeneous or inhomogeneous scalar held 0 evolving slowly down its potential V (0). The 
resulting cosmological models are known as quintessence models [Caldwell et al, 1998; 
Turner et al, 1984, 1997]. In these models, the scalar held can be seen as a perfect huid 
with a negative pressure given by p = (7 — l)p (0 < 7 < 1). In what follows one shall 
focus on equations of state with a constant 7 and disregard the case of time varying 
equations of state invoked in some quintessence models, called tracker models, to solve 
the cosmic coincidence problem [Zlatev et al, 1999; Cooray et al, 1999]. 

Going on with the analysis, very interesting question is whether it is possible to con¬ 
struct a successful common scheme for the two cosmological mechanisms involving rolling 
scalar helds i.e., quintessence and inhation. This perspective has the appealing feature 
of providing a unihed view of the past and recent history of the universe, but can also 
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remove some weak points of the two mechanisms when considered separately. Indeed, 
inflation could provide the initial conditions for quintessence without any need to £x 
them by hand and quintessence could hope to give some more hints in constraining the 
inflation potential on observational grounds. 

A homogeneous and isotropic universe is characterized by the Friedmann-Robertson- 
Walker (FRW) line element 

ds^ = dt^ — a^{t) ( — siri^Odcf)^ (l-H) 

— kr^ J 

where a is the scale factor and k is the spatial curvature index (i.e., /c = 0 , ± 

1). Assuming the universe is hlled with perfect normal matter plus quintessence fluid 
corresponding to some scalar held governed by Klein-Gordon equation, the over all 
stress-energy tensor of the cosmic huid without the dissipative pressure has the form 


= (p + v)upiiy + pg^,y, = -1 (1.12) 

where p = and P = P^ + P^- Here p^ and p^ are the energy density and 

pressure of the matter whose equation of state is p^ = ( 7 ™. — l)p^ with adiabatic index 
in the interval 1 < 7 m < 2. Likewise p^ and p^, the energy density and pressure of the 
minimally coupled self-interacting Q-matter held ip i.e.. 


P, = ^ (VO, p, = Ip - 


(1,13) 


are related by an equation of state similar to that of the matter, viz., p^ = ( 7 ^ — l)p^, 
so that its adiabatic index is given by 

ijj^ 


Ik,2 


(1.14) 


2^2 + V{lP) 

where for non-negative potential V{ijj) one has 0 < 7 ^ < 2. The scalar held can be 
properly interpreted as Q-matter provided 7 p < 1 [Zlatev et al, 1999; Steinhardt et al. 
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1999] . As usual an overdot means derivative with respect to cosmic time. In general 
varies as the universe expands and the same is true for 7 ^, since the massive and 
massless components of the matter fluid redshift at different rates. 

The Friedmann equation together with the energy conservation of the normal matter 
fluid and quintessence (Klein-Gordon equation) are 

+ “7 = q(dm + Pv) (1-15) 

CL O 

Pm + ^HjrnPm = 0 (1-16) 

^ + + = 0 (1.17) 

d(p 

where H = ^ denotes the Hubble factor. Introducing 17 ^ = = 

— and H with = 3H^ as the critical density, the set of equations 

(1.15)-(1.17) can be recast as [Ellis et al, 1997] 

fim + ^0 + Hfc = 1 (1-18) 

H = H(H-l)(37-2)i7 (1.19) 

hp=[2 + (37 - 2 )H - 3-fp]npH ( 1 . 20 ) 

where 7 is the average adiabatic index given by 

7 'Im^m T Ap^p ■ 

The combined measurements of the Cosmic Microwave Background (CMB) temper¬ 
ature fluctuations and the distribution of galaxies on large scales seem to imply that 
the universe may be flat or nearly flat [Turner, 1997; Ostriker et al, 1995]. Hence the 
interesting solution at late times of equation (1.19) is H = 1 (i.e., /c = 0) and discard the 
solution H = 0 as incompatible with observation. The solution H = 1 is asymptotically 
stable for expanding universes {H > 0) provided that the condition < 0 holds in a 
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neighbourhood of = 1 and this implies 7 < 2/3. Hence the matter stress violates the 
strong energy condition p + 3p < 0 and so the deceleration parameter 

ad 

?=-^<0 ( 1 . 22 ) 

(because, ^ = — |(p + 3p) > 0). Hence as a consequence the universe accelerates its 
expansion. 

From theoretical point of view a lot of works [Caldwell et al, 1998; Ostriker et al, 
1995; Peebles, 1984; Wang et al, 2000; Perlmutter et al, 1999; Dodelson et al, 2000; 
Faraoni, 2000] have been done to solve the quintessence problem and possible candi¬ 
dates for Q-matter are cosmological constant (or more generally a variable cosmological 
term), a scalar held [Peebles et al, 1988, 2002; Ratra et al, 1988; Ott, 2001; Hwang 
et al, 2001; Ferreira et al, 1998] with a potential giving rise to a negative pressure at 
the present epoch, a dissipative huid with an effective negative stress [Cimento et al, 
2000 ] and more exotic matter like a frustrated network of non-abelian cosmic strings or 
domain wall [Bucher et al, 1999; Battye et al, 1999]. 

Unfortunately, most of the helds (Q-matter held, tracker held) work only for a spatially 
hat {k = 0) FRW model. Recently, Cimento et al [2000] showed that a combination of 
dissipative ehects such as a bulk viscous stress and a quintessence scalar held gives an 
accelerated expansion for an open universe {k = —1) as well. This model also provides 
a solution for the ‘coincidence problem’ as the ratio of the density parameters corre¬ 
sponding to the normal matter and the quintessence held asymptotically approaches a 
constant value. Bertolami and Martins [2000] obtained an accelerated expansion for the 
universe in a modihed Brans-Dicke (BD) theory by introducing a potential which is a 
function of the Brans-Dicke scalar held itself. Banerjee and Pavon [2001] have shown 
that it is possible to have an accelerated universe with BD theory in Friedmann model 
without any matter. 
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Chapter 2 


Brans-Dicke Cosmology in an 
Anisotropic Model when Velocity of 
Light Varies 

2.1 Prelude 

In recent years theories with varying speed of light (VSL) has been attracted considerable 
attention [Albretch et al, 1999; Barrow et al, 1998, 1999, 2000; Magueijo, 2000, 2001, 
2003] due to its ability to solve the so-called cosmological puzzles - the horizon, flatness 
and Lambda problems of big-bang cosmology. All previous attempts to overcome these 
difficulties invoke the basic inflationary form, where the observable universe experiences 
a period of ‘superluminal’ expansion. Here, one has to modify the matter content of the 
universe such that ordinary Einstein gravity becomes repulsive and halts the exponen¬ 
tial expansion. One can however resolve to theories with varying speed of light as the 
alternative method to solve the above mentioned cosmological puzzles. Here, instead of 
changing the matter content of the universe, one has to change the speed of light in the 
early universe. The universe is assumed to be radiation dominated in the early stage 
and the matter content is assumed to be the same as in the standard big-bang (SBB) 
model. So the geometry and expansion factor of the universe go in accordance with the 
SBB model. 

So far some works have been done with VSL by Magueijo and co-workers [Albretch et 
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al, 1999; Barrow et al, 1998, 1999, 2000; Magueijo, 2000, 2001] and others [Jacobson et 
al, 2000; Landan et al, 2000; Besset et al, 2000; Kiritsis, 1999; Harko et al, 1999] spe¬ 
cially in isotropic cosmological models (which also inclnde BD-theory). Very recently, 
Magneijo [2001] has investigated, the possibility of black holes formation by stndying 
spherically symmetric solntions with VSL. 

In this chapter, we stndy Brans-Dicke Cosmology in anisotropic Kantowski-Sachs 
space-time model, considering variation of the speed of light. We address the flat¬ 
ness problem considering both the cases namely, (i) when the Brans-Dicke scale held 
(j) is constant (ii) when (j) varies, specially for radiation dominated era pertnrbatively 
and non-pertnrbatively. For radiation dominated era we give the exact solntion to the 
hatness problem. We stndy solntions to the lambda, qnasi-lambda and qnasi-hatness 
problems and examine their asymptotic behavionr. 

2.2 The Basic Equations 

For anisotropic Kantowski-Sachs (KS) model with metric ansatz 


ds'^ = —c^dt^ + a^{t)dr‘^ + h‘^{t)dVt\ , 


the Brans-Dicke (BD) held eqnations with varying speed of light are 



2 


( 2 . 1 ) 



( 2 . 2 ) 


and the wave eqnation is 



(2.3) 
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Here the velocity of light c is an arbitrary function of time, u is the BD coupling 
parameter and the BD scalar held is 0 = From the above held equations we have the 
non-conservation equation 


D + ('’ + 1 


cc 


(2.4) 


Now p = \pc^ is the equation of state in the radiation era for which the general 


solution 0 is (from equation (2.3)) 


, f dt 


(2.5) 


The above held equations in BD-theory have been formulated in Jordan frame. To 
switch over to Einstein frame we shall have to make the following transformations: 


di = a/ G(j) dt, d = \/G(j) a, b = \/G(j) b, a = ( a; + 


1/2 


ln(G0), 


P = (^0) p, p= p , 


and the above held equations become (c is treated as constant) 


( 2 . 6 ) 


I/I rii 


^ + 22^ = -JttG \p+% 


a 


/2 


(2.7) 


and 


'/2 


^d'b' ^ a- 

■ I + 2—— + — — SuGp H—— , 

b a h 2 


d' rfi'\ , SttG 


a” + \ - + 2^\ a' = 
a b 


a/ 6 + den 


P 


3p 


( 2 . 8 ) 


(2.9) 


with the prime “ ' ” and “ . ” represent the diherentiation with respect to t and t 
respectively. 


One can interpret these held equations as standard KS equations with constant G and 
a scalar held a is added to the normal matter . The scalar held behaves like a “stih” 
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perfect fluid with equation of state 


Pa = Pa 


IGttG 


( 2 . 10 ) 


If the velocity of light is constant, then in Einstein frame total stress-energy tensor is 
conserved but there is an exchange of energy between the scalar held and normal matter 
according to the following equation 


p' + 





\/6 + iuj V / 

( 2 . 11 ) 


On the other hand, if the 
conservation” equations 


velocity of light varies then we have two separate 


“non- 


P 


“'- + 2-'’' 

a b 


, P 


a 


\/Q + Auj 


P 2 


3p \ cc 

' + 


A'lrGb'^ ’ 


and 


Pa + 


a' V 

- + 2 - 
a b 


. , Pa 

Pa + — 


a 


3p 
P~ — 


( 2 . 12 ) 


(2.13) 


a/6 -l- 4a; 

Thus standard KS model held equations with constant G for VSL are also the held 
equations for BD theory in Einstein frame but here one must add to normal matter the 
scalar held energy and pressure. Therefore, the total energy density and pressure are 
given hj pt = p + Pa and Pt = P + Pa respectively. 


2.3 Perturbative Solutions to the Flatness Problem 


For a possible solution to the hatness problem in VSL theory, let us hrst study solutions 
when there are small deviations from hatness. The critical energy density pc in a BD 
universe is given by the equation 



+ 2-J 

a b 


^Tipc 

0 



0 ^ A 
0 


2 


(2.14) 
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In Einstein frame the critical energy density for normal matter is 


Pc = 


8tiG 


b ) 0, b 


a 


12 


Pc 

G202 


Hence, if we define the total critical energy as 


Pa Pc A' Pa ) 


then from (2.15) we have 


pa. 


SttG 


+ 2 - 


a b 


Let us dehne a relative flatness parameter as 

Pt Pa 


— 


Pc 


whose evolution equation is given by 


(2.15) 


(2.16) 


(2.17) 


( jy \ 

+ ^t) ( 'J— + 2(7 — 1)Y I +2— it 

\ a b J c 

Here 7 — 1 = pjpc? = p/pc? , is the equation of state, with 7 constant. Since 

£ 


(2.18) 


l + {2uj + 3)07327r0pe 
in Jordan frame, the deviation from flatness is given by 


(2.19) 


6 = 


P- Pc 


-< e (2.20) 

pc + {2uj + ?>)cj? 18271 cj) pc 

where e = (p—pc) / Pc measures natural deviations from flatness. Now J, the adaptation 
parameter satishes the differential equation 


i = 5(1 + i) h- +2(7-1)7 5 ( 37 -2)t) + 2 L 

'a 0 2 CD c 


( 2 . 21 ) 


If 6 is assumed to be very small compare to unity i.e., J <C 1, then neglecting the 
square term in ( 2 . 21 ), one can integrate the above to get 


( 2 . 22 ) 
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where is the integration constant. Hence we have 

1 ^_£o_ (2a; + 3)0^ 

e a'^62(7-i)0(37-2)/2g2 327rpc£ 

In BD theory G oc 0“^ , so to solve the flatness problem 0 should decrease in the early 
universe and the first term on the right hand side must dominate over the second one. 

In particular, in radiation era (7 = 4/3) the above equation (2.23) simplifies to 


1 ^ ^0 _ ( 2 a; + 3)a^ 

e 32Ti(j)a‘^b^pce 

where we have used equation (2.5) to obtain cf . If a is small and c = cq 
then [Barrow et al, 1999] it is still possible to solve the flatness problem for n < 

—1 — . Furthermore, since pc oc l/(a 6 ^)'^/^ at late times, so the second term 

is always very small compare to the first and hence may be neglected. 


2.4 Non-Perturbative Solutions with (p = Constant: Matter 
and Radiation Dominated Era 

When 0 = 00 (constant) we have 5 = £ in (2.20) and the differential equation for 6, 
( 2 . 21 ) can be written as 

e = e{l + e)(^- + 2(7 - l)M + 2- e (2.25) 

\ a 0 I c 

Now, if we assume a power-law form in the scale factors of the velocity of light then 
the above equation becomes 

i = + e) (0 + 2(, - 1)^) + 2e (^) + 2(„ - l)t) (2.26) 

In order to get an exact analytic solution to the above equation, we assume £ <C 1, so 
we can neglect £^-term. We then have the solution 


S — £()(l 


7+2n^/(3n—2)^2(7—l)+4n(n—l)/(3n—2) 


( 2 . 27 ) 
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with Sq as the integration constant. 

Furthermore, without assuming any restriction on the parameter £, if we assume the 
two arbitrary constants n and 7 to be equal then also we have an exact analytic solution 
for 


1 


(2.28) 


(57-2)7(37-2) 


+ A [a'>62(7-1)] 


where A is an integration constant. 

2.5 Exact Solutions to the Flatness Problem: Radiation Dom¬ 
inated Era 

To obtain the exact solutions to the flatness problem we assume that the variation of 
light is governed mathematically by the relation [Barrow, 1999] 


c = Co 


n 


(2.29) 


However, in the Einstein frame the above relation reduces to c = Cq oAmA 1 ) 
and equation (2.12) becomes 


l/(3n-2) 



(2.30) 


The equation of state p = (7 — l)pc^ simplihes the above equation further to 



(2.31) 


An exact integral can be obtained (assuming n = 7 ) as 



(57-2)7(37-2) 


(2.32) 


Hence for the radiation dominated era, we have in Jordan frame 


= B + 


(2.33) 














CHAPTER\^ Brans-Dicke Cosmology in an Anisotropic Model when Velocity of Light Varies23 



Fig 2.1 ; Here the flatness parameter e in (2.27) has been plotted against t for the radiation 
era (7 = 4/3) choosing different values for the constant n. The scale factors here are assumed 
to be in simple expanding form. 



Fig 2.2 : This figure shows the variation of the flatness parameter from (2.28), where there 
is no restriction on e. This figure is also drawn for the radiation era and assuming as before 
the simple expanding form for the scale factors. The variation of e for different choice of the 
arbitrary constant A is shown in the figure. 
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where we have used = p{ab‘^y^^ . 

Moreover, in the Jordan frame we have equation (2.29) for the velocity of light and 
we can write 


c 

c 


n 

3n — 2 


n —h 2[n 
a 




n ( 

+ 2; 


(2.34) 


Thus, if 0 is a decreasing function then |c/c| < \ [na/a +2{n — l)b/b]/{3n — 2) \ and we 
can therefore conclude that a varying speed of light in the early universe (with weaker 
gravity) helps to solve the flatness problem. 

Alternatively, for general equation of state p = (7 — l)pc^ , if we make the transfor¬ 
mation 




(2.35) 


and assume the variation of the velocity of light as 


C = CqX 


nl2 


then the equation of continuity (2.4) becomes 




, a ^b\ ncicpx'" "x 
p + 7P I - + 2- ) = 


(2.36) 


Integrating this equation gives 


pa 


T'627 = b + 


8n{n -1- 1) 


_i \-| n+l 


(2.37) 


(provided n 7^ —1) which gives the evolution for general 7 . 


2.6 Solutions to the Quasi-Flatness Problem 

2.6.1 When BD Scalar Field 0 is Constant 

In this section we shall investigate whether it is natural to have evolution which asymp¬ 
totes to a state of expansion with a non-critical density . Let us start with the continuity 
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equation for constant 0 and assume as before the equation of state to be 


P 

pc^ 


= 7-1 


1 


and the velocity of light to be 


c = 


Co 


n/(3n—2) 


then the integral gives 


^ 47r(57 - 2) L J 


(It is to be noted that in order to obtain the integral one has to assume n 
Now substituting this value of p in equation (2.2) we have 


ah Stt B 

^ ah cj) {ah'^)^ 


+ 


jclf, 2 (57-2)7(37-2) 


47r(57 — 2) 


+ 


{ah^y 


-Icf, 


In particular, for radiation era the above equation simplifies to 


2 l( = 


8nB 


ah (f){ah‘^Y/^ 7 

If n is the density parameter, then it can be defined as 

n SttGp 


5^ 


n - 1 cV62 

We note that the ratio between the two terms on the right hand side is 
constant for a quasi-flat open universe. However, for 7 = 4/3, we get 


n 


StiGB 


(2.38) 

7)- 


(2.39) 

(2.40) 

(2.41) 

almost a 


n-1 7 c2(a26)W9 


(2.42) 
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For expanding universe a, b grows with time. So asymptotically the second term on 
the R.H.S will be negligible and we have 



which leads to an open universe with hnite hi value today. 


2.6.2 When BD Scalar Field Varies 

For the radiation era (7 = 4/3), let us dehne 

y = z = (j){a^by^^, c = cqz^' 

So from the held equation ( 2 . 2 ) we get (using equation (2.5)) 


y' z' z'^ 

2 --y 

y z z^ 


327rp 


{ab^f/^ - 


4c^ z 





Now, from equation (2.38) the expression of p takes the form 


(2.43) 


(2.44) 


piab^f/^ 


BmcQz‘^^^^ 
47 r( 2 m + 1) 


In addition the critical density pc is obtained from the differential equation 


(2.45) 




J2 


32/ 


y z z^ 

Hence the expression for the density parameter is 


2u;\ 


(2.46) 


or equivalently. 


pc 


Q 


+ ¥f-(i + ¥) 


y z 


21LSL _ _ (1 + ^) si 

y z z‘^ \ 3 / V 


p 


SnB 


n - 1 p- Pc 
Thus, if 2m + 1 < 0 then as t —>■ 00 , z 


^2 ^2m+l 
CqZ 


+ 


2m 


2m + 1 
00 so we have hi 


(2.47) 


1. But if 2m + 1 > 0 


then as f —> 00 , hi —>■ —2m. So asymptotically, we have a quasi-hat open universe. 
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2.6.3 General Asymptotic Behaviour 

In this section, we consider general asymptotic behaviour when 2m + 1 > 0 i.e., for 
quasi-flat open universe. The equation (2.44) can be approximated by 

2y'z' - ^ ~ r ^2(m+i) 
z 

which has a first integral 

y _ [T ;22(m+l) 

Thus, if we assume 2 ; ~ 7 - 2 /!^oo^ then y ~ and 

/ 1-2/O00 \ 

(2.48) 

For the case of hloo < 1, we have 0 —>• 0o, a(r) ~ and 6(r) ~ in the limit 

as r —> cx), or equivalently, 

as t —> cxo. 

In addition, when hloo = 1, we get a I) as expected for flat radiation 

asymptote. 


2.7 The Lambda and the Quasi-Lambda Problems 


In this section, we shall examine the effect of variation of speed of light when a cos¬ 
mological term is introduced in the BD held equations. In fact, the incorporation of 
a cosmological term is equivalent to introduction of a vacuum stress which obeys an 
equation of state 


PA = 


Pa 


with 


Ac2 


Thus equation of continuity can be generalized to 
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, , (a b\ cc 

(p + Pa)+ + 

Also the field equation ( 2 . 2 ) can now be written as 

ah 

- + 2-- = SvrGp - Tj + Ac' (2.49) 

m ab ¥ 

Now, using the above expression for vacuum stress-energy and using 7 — 1 = p/pc? as 
the equation of state for the matter and assuming 

c = co 

the above equation of continuity can be integrated (assuming n = 7 ) to give 


B 


P = 


7C^ 


[a^62(7-i)] 3^=2 _ 2) 1 


(a6')^d[a^6'A-i)]3|^ 

(2.50) 


{a¥p 87rG'(57 — 2) {a¥y 


SttG a¥ 


On substituting this value of p in the above field equation (2.49), we have 


i? ah 


»„GB 2^,4 

(a62)7 + (57 - 2) ( alf )-< 




27 /( 37 - 2 ) 


-A<;(( 57 - 2 )^^ j ( 2 , 51 ) 

If we define the generalized density parameter as 




87tG{p + pa)&^ 


then using the density parameter O from equation (2.41) we have 


_ O _ 87rG{p + pA)b^ 


(2.52) 
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If we now substitute the value of p from equation (2.50) we obtain 

^ _ 27 (47-3)7(37-2) SnGB ^^^^ 2 ( 7 -!) j-(47-2)7(37-2) 

1 5'y 2 Cq 


+ M? 


-A (37 - 2) [a^62(7-l)j 2)7(37 2) j ^^^2^7^[a752(7-l)]277(37-2) (' 2 _ 53 ) 

when 7 > 1. The 2nd, 3rd and 4th terms are negligible compare to the hrst term, 
hence we have 

fl - 1 57 - 2 L J ^ ^ 

So, for large a, b (i.e., a,b ^ 00 ) 

n 

—- —>• cx) i.e., f2 —>• 1. 

Furthermore, if A = 0 then 


^ ^ 27 r 7^2(7-l)l ( 47 - 3 ) 7 ( 37 - 2 ) SttGB (47-2)7(87-2) 

- 1 57 - 2 L J cl ^ ^ 

So for the radiation era (7 = 4/3) 


n 

n - 1 


8 

14 


Cq 


(2.56) 


As a, 5 —>■ CX 3 , f 2 /(fl — 1 ) —> 8/14 i.e., 12 = —4/3. This asymptotic behaviour is same 
as in Sec 2.6.1 where for quasi-flatness problem the BD scalar held is assumed to be 
constant. 


2.8 Discussion 


In this chapter, an extensive analysis of the BD solutions for anisotropic cosmological 
model has been done where there is a variation of the velocity of light with time we have 
assumed the velocity of light to be in power-law form in the sacle factors. The hatness 
problem has been discussed in details. Here perturbative, non-perturbative and exact 
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solutions of flatness problem have been obtained. The graphical representation for non- 
perturbative solution in the radiation era has some interesting feature. Quasi-flatness 
problem with general asymptotic bahaviour has also been discussed and in most cases 
we obtain an open universe. Finally, we have also studied the lambda and quasi-lambda 
problem in VSL and the asymptotic behaviour is very similar to that in the flatness 
problem. 


Chapter 3 


The Cosmology in a Perfect or 
Causal Viscous Fluid with Varying 
Speed of Light 

3.1 Prelude 

Scalar dissipative in cosmology may be treated via the relativistic theory of bulk viscos¬ 
ity [Maartens, 1995; Zimdahl, 1996]. The causal and stable thermodynamics of Israel 
and Stewart provide a satisfactory replacement of the unstable and non-causal theories 
of Eckart and Landau and Lifshitz. If viscosity driven inflation occurs (leaving aside for 
the moment various questions about the hydrodynamical consistency of such models) 
then this necessarily involves non-linear bulk viscous pressure [Maartens, 1995]. A non¬ 
linear generalization of the Israel-Stewart theory is more satisfactory model of viscosity 
driven inflation. We hnd exact solutions in a flat and non-flat FRW universe for simple 
and consistent thermodynamic co-efficients and equation of state. Provided the ther¬ 
modynamic parameters satisfy certain conditions, the model admits thermodynamically 
consistent inflationary solutions. 

There has been a recent trend of investigation where one specihes a particular form of 
the expansion of the universe and then seeks for the matter or any other held distribution 
that would lead to such a temporal behaviour of the scale factor. The inhationary uni¬ 
verse scenario can solve some of outstanding problems [Guth, 1981] of standard big-bang 
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cosmology. Ellis and Madsen [1991] have considered an FRW model with a minimally 
conpled scalar held along with a potential and perfect hnid in the form of radiation. 
For brief bnt comprehensive reviews we refer to the papers by Maartens [1995], Hiscock 
and Lindblom [1983] and Lindblom [1996]. In recent commnnication, Zimdahl [1996] 
discnssed how the trnncated theory is indeed a good limit to the fnll cansal theory in 
certain physical sitnations. Mandez [1996] presented a similar work where attempts had 
been made to hnd 0 and V both fnnctions of time t and thns V = V (0) for an expo¬ 
nentially expanding nniverse. 

In this chapter, we examine the behavionr of 0 and V for power-law inhation of scale 
factor R{t) and velocity of light c{t) which varies with t [Albretch et al, 1999; Barrow 
et al, 1999]. We generalize the investigation of Ellis and Madsen [1991] and Banerjee 
[1998] to inclnde a hnid and we do not ignore the presence of the hnid in the calcnlation 
for the expression for the scalar held 0 and the potential V in terms of time t and hence 
obtained V = F(0) for a variety of scale factor and velocity of light. 


3.2 The Model in Gravitational Field Equations 


The energy momentnm tensor for a bnlk viscons hnid and minimally coupled scalar held 
0 are given by 

= {pc^ +P + T^)Vf,V^ + (p + T^)g^u (3.1) 


and 




-0,„0'“ + F(0) 


(3.2) 


where p is the energy density, p the thermodynamic pressure, vr the bulk viscous 
pressure, F(0) the scalar potential and the four velocity satisfying the condition 
= Sq and = —1. 
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Consider, the line-element in FRW space time, 


ds^ = -c^dP + R\t) 


dr"^ 


1 — kr"^ 

The gravitational held eqnations are given by 


-I- r^{d6^ -I- sin^Odcj) ) 


(3.3) 


3 



Vf‘ + ¥(? + p 


(3.4) 


and 

2// + 3//^ + ^ = -ly + Vc^-hJ^ (3.5) 

n‘^ 2 

where, Hnbble constant H = ^ and k is the cnrvatnre parameter. (We have chosen 
the nnit SttG = 1 ). 

The ‘conservation eqnation’ for the hnid is 


p + 3H = 4^ p (3.6) 

where we have assnmed the velocity of light c is an arbitrary fnnction of time only. 
The eqnation of state for the hnid distribntion is 

P = (7 - l)pc^ (3.7) 

where 7 is the constant adiabatic index of the hnid . 

The cansal evolution equation for the bulk viscous pressure is given by 


TTT + TT = —3pH -rvr 3H -|- 


T i] T 

T p T 


(3.8) 


where r is the ralaxation time, p is the co-efhcient of bulk viscosity and T is the 
temparature of the hnid. For r = 0, (3.8) reduces to the non-causal equation vr = —3pH. 
In (3.8), e = 0 gives the truncated theory (extended irreversible thermodynamics [FIT]), 
while e = 1 gives the full Israel-Stewart-Hiscock theory. It can be noted that in full Israel- 
Stewart-Hiscock theory the derivations from equilibrium are small so we can assume 
| 7 r| -Cp = {c)pP [Zimdahl, 1996; Maartens et al, 1997]. 
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Using equation of state, the field equations (3.4) and (3.5) becomes 


= -2H + 


2kc^ 71 

”r2“ 77 


(3,9) 


and 


U = 


H + 3H^- 


7 


P + 


TT 


+ 


2k 


(3.10) 


The power-law form of scale factor R and the velocity of light are assumed as 


R = Rot^, c = Cot™ (3-11) 

where, Rq, cq, n are positive constants and m is real constant with the restriction n > 1 
(for accelerating and expanding universe). 

For power-law expansion, (3.9) and (3.10) reduces to 


^2 _ 2 n ^ 2/cCQ^2(m-n) 

fir 


TT 


7P 




(3.12) 


and 


U = 




3n^ — n f 2 — j 


P + 


TT 


2clR^ 


+ 


2k 

RlR^ 


(3.13) 


Since > 0, this implies. 


n -fp ^ > n 

R Rf 2 2clR^ - 


for the model to be consistent. 


(3.14) 


3.3 Solutions for Perfect Fluid 

For a perfect fluid, tt = 0, and the conservation equation (3.6) gives 

p = (3.15) 


where, po is an integration constant. 
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In this case, consistent condition (3.14) becomes 


ff_ I _ 7P0 ^4m-37n ^ r, 

Rf 2 

For k = 0, the model is valid for 


(3.16) 


t < 


\ l/(4m—37n+2) 
2n 1 

7P0 ) 


, if 4m — 37/1 + 2 > 0 


OR 


/ \ l/(4m—37n+2) 

t>(^] , if 4m- 37/1 + 2 < 0 

and for k = +1, —1 the model is also valid for restricted period of time. 

In absence of any fluid ( i.e., po = 0) fhe condition (3.16) satishes for A; = 0, +1 for all 
time t and for k = —1, condition (3.16) satishes for 

, , ,+i?gV/2(m-n+l) 

t < I -^1 , if m > n — 1 


OR 


, . ,+i?gV/2(m-n+l) 

t > I j , if m < n — 1 


For perfect huid, (3.12) reduces to 


02 = ^ + _ ^pQt4m-37n 

t Rq 


(3.17) 


Generally, it is very difficult to integrate this equation to 0 = 0(t). So we are assuming 
two cases : 

(i) A; = 0 ; In this case (3.17) integrates to 


M +log 




7P0 j-4m-l-2— 8777 . _ ^ 




(3.18) 


where 0o is an integration constant. 
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The potential V is given by 




-'0 I 


{3n^-n ^ ^ , 2 m- 3 ^n 


^2m+2 


2 J 


PoC 


For simple case, n = l,m = 1/4 and for radiation 7 = 4/3, we have 


(3.19) 


-00 = - 2 ^ 2 ^/!-^ + 72 log 


2V2jl-^ + l 


272^/1-^-1 


and 


1 / = 




1-^ 

6t 


(3.20) 


(3.21) 


(ii) k ^ 0 : We shall discuss for m = n — 1. In this case (3.17) yields to 


[(4 - 37)n 


2](0-0o) = v^ 



JPo 

2B 


^(4-37)71-2 _j_ jQg 


^1 - mt(4-37)n-2 _ 1 ■ 
- mt(4-37)n-2 ^ X 

(3.22) 


and (3.13) gives 

r = 1|2B + 3n(n - l)|i + (^1^) ( 3 . 23 ) 

where B = n + 

Hq 

Similarly, we can choose, 2(m — n)= Am — S'yn OR 4m — 3'yn = —2 and we have same 
form of 0 and V as well as (3.22) and (3.23). 


3.4 Solutions for a Fluid with Bulk Viscosity 

A causal viscous fluid with bulk viscous stress is given by equation (3.8) and we shall 
discuss the model for e = 1. For the sake of simplicity, we take the thermodynamic 
quantities to be power function of fluid density p i.e., 

T] = apP, T = [fp\ r = ^ = (3.24) 
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where a, [3, 3, /i, r are positive constants. 

On differentiating eqnation (3.6) and nsing eqnation of state, we have 

p + 3hr (qp + J) + 3hf (^7P + ^ - = 4 - ^p^ (3.25) 

Then eqnation (3.8) rednces to (using equations (3.6) and (3.25)) 


P 


1 + 3n + 


2m(7 — 1) 


7 




9 n ^7 „ , 6mn\ p 5 i . 

-^ 

2 j / P a 


+ ( 3 n 7 — 4m) 


6 

a t 


27 


27 


P 




P 


-j-2m+2 


= 0 


(3.26) 


where we have assumed the standard thermodynamic relation between temparature 
and density for barotropic fluid, 


^ 7-1 

T ~ p T- 


(3.27) 


Now comparing T between (3.24) and (3.27) we obtain 

7 — 1 

r = - - (3.28) 

7 

To complete the solution for p, we need some cases: 


Case / : p 7 ^ 1 and p 7 ^ 1 — m. 

For p 7 ^ 1, the solution of the equation (3.26) is of the form 

2 m+l 

P = Pot r-i (3.29) 

with the condition 

(2m + l)[ 6 n 7 + 4 m (7 — 1) + 1] + (p — l)( 9 n^ 7 ^ — Ihm^q — 12mn) = 0 (3.30) 

where the constant po is given by 


Po 



QarP 


{ 2 m + 1 

V P-1 


+ 377-7 — 4m 


1 

p-1 


(3.31) 
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Condition (3.14) alongwith equations (3.6) and (3.29) yields 


n , ^.2(m-n) ,P^( ‘^rn + l 
Rq 6n \ /i — 1 

This condition is satisfied for all values of A;, {k - 
of time. For fc = 0, the model is consistent for 


— 4mJ t >0 (3.32) 

: 0, +1, —1) only for restricted period 


^ ~ _po[2(a - l)(/i - 1) - a] 

where a = + 2. 

ii-i 

Using equation (3.6), the equations (3.12) and (3.13) reduces to 


(3.33) 


^ ^ ,2(m-n) _ do 

R Rq 3n 



2m + 1\ 
F-1 ) 


2m+l 

t f'-i 


(3.34) 


and 


V 


1 




3rR — n 


El 

6n 


^4m — 6n + 


2m + lA 

) 


t ''-1 


2k 


(3.35) 


Two cases are considered as follows: 

(i) /c = 0 : Equation (3.34) can be integrated to yield, 


0 — 00 



2a/ 1 — + log 


Vl - AC -1 
a/1 - AC + 1 


(3.36) 


where, a = + 2, A = ^ ^4m — ^ind 0o is an constant of integration. 

Now (3.35) can be written as 


U = 




+ «L(4m -6„ + a- 

R 6n 


(3.37) 


It is difficult to write t = f(0) from (3.36) and from (3.37) we cannot write V = U(0) 


analytically. 
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V 



Fig 3.1 : Here the potential function V has been plotted against cf from equations (3.36) and 
(3.37) for radiation era (7 = 4/3) with po = l,co = 1, and cpo = 1. 

(ii) k 0 : In this case, we shall discuss the model for m = n — 1. After integration 
(3.34) we have 


~ 00 — 


Vw 


nA 

1 - —C + log 




1 - + X 


(3.38) 


kcl 


where B = n + a = + 2 and A = — a — 2). 

Also (3.35) reduces to 


1 / = 


-^0 L 


2B + 3n{n — 1) — ( nA + —— 

3n 


1 


(3.39) 


Similarly, if we choose = 2(m — n), then we have same form of 0 and V as well 
as (3.38) and (3.39) respectively. 


Case II: fi = \ — m. 

(i) A; = 0 : Equation (3.34) integrates to yields 


\/2n{l — A) logt 


(3.40) 
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V 



Fig 3.2 : The potential function V in equation (3.42) has been plotted against 4> for radiation 
era (7 = 4/3) with n = 1, m = —3/16, po = 1, cq = 1, and ipo = 1. 


and equation (3.35) becomes 


1 / = 


V/i 


(3.41) 


^ 2 {m+l) 

where Vq = [3n^ — n + — 3n — 1)] and A = + 1). 

From equations (3.40) and (3.41), we have the potential V is an exponential function 
of the scalar field (j) i.e., 


2(m+l)(<^-<(>0) 

V = V/e (3.42) 

(ii) k 0 : For k ^ 0, two cases arises : 

(a) m = 77 . — 1 : Equation (3.34) integrates to 


and equation (3.35) yields 


'2n{l-A) + 


2kCn 


\ogt 


(3.43) 


^ = ^0 + 


2k 


Thus V can be expressed as a function of 0 in the form 


(3.44) 


2k 

^ = ( ^0 + ) Exp 


2 n (0 - 0 o) 


2 n(l -A) + ‘^ 


(3.45) 
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ip) m ^ n — 1 ■. Equation (3.34) integrates to 


(j) — po 


y^2n(l - A) + ^t2(m-n+l) 
2{m — n + 1 ) 


+ \/2n{l - A) logt 


and equation (3.35) gives 


(3.46) 




t2(m+l) ^2^2n 


(3.47) 


In case ( 6 ), it is very difficult to eliminate t from (3.46) and (3.47) to the form 
E = E( 0 ). 

Case III: = 1. 

In this case relaxation time becomes a constant. 

(i) m = —I : In this case, the solution of equation (3.26) is of the form 


-3n7-22^ _25j_, 

P = Pot '=0 6 “ 


(3.48) 


with the condition : 


^ 9n^a 
3n7 H- — 


OR 




+ 7(2 — 3n) — 2 + 9n 7 + 6 n — 47 = 0 


_(2+3n7)+^ 

P = Pot 


(3.49) 


(3.50) 


with the condition: 


27 rPa^ — 6 nacQ{'y + 1) + 2cg(7^ — 7 + 1) = 0 


(3.61) 


The expression for p in (3.48) yields a complicated expression for cjp from equation 
(3.12) which cannot be analytically integrated. 
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Now we shall discuss for the solution (3.50). Using (3.50) the equations (3.12) and 
(3.13) reduces to 


2n ^ 


2kcQ ^ / 3na \ ,-( 2 + 3 n 7 )+= 


i?2^2n+l 


+ Po 


\ r 2 

V ^0 


-7 t 


and 


U = 


-0 L 


3n^ — n f If 3na \ \ -( 2 + 3717 )+^^ 
- 1 - Po f 1 + 7 1-^-7 ) )^ 


+ 


2k 


(a) k = 0 : Equation (3.52) integrates to the form 


\f 2 fn 


X 


2 a /1 + wt^ + log 


a/1 + — 1 


Vl + ujt^ + 1 


where x = 3n — 7 j and w = ^x. 
and (3.53) yields 


(3.52) 


(3.53) 


(3.54) 


U = 


-0 L 


3n^ — n 


t 


- (po + nw)t 


X—1 


(3.55) 


It is very difficult to eliminate t from (3.54) and (3.55) to the form V = V{(j)). 

{b) k ^ 0 : In this case, we shall discuss for 3nQ; = 7 Cq. Equation (3.52) integrates to 


— 00 — 




1 - 2 n 
and (3.53) reduces to 


„ 3ka , „ 






..1/2 xl 2P 

V = -2 (3n^ - n - Po)- + 


t 


In this case also, we cannot express V in the form V = V{(j)). 
(ii) m = —1 : The solution of equation (3.26) is of the form 


(3.56) 


(3.57) 


(3.58) 


with the condition; 


+ 2 y [{ 3 n — 2)7 — 2] + 9n^7^ — I67 + 12n = 0 


(3.59) 
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where z= and ±yz = (1 ^ ; 2)[2 + 7 ( 3 n - 1)] - ( 3 n 7 + 4). 

The expression for p yields a complicated expression for 0^ from equation (3.12) which 
cannot be analytically integrated. 

3.5 Discussion 

In this chapter, we have considered perfect fluid distribution with a scalar held 0 having 
potential V{(j)) for isotropic space-time model, considering variation of the velocity of 
light. We have considered the power-law inhation and assumed a polynomial form (in 
time) for the velocity of light. In this huid distribution V cannot be expressed as a 
function of (j) explicitely. For simple case (n = 1, m = 1/4) and radiation dominated 
era (7 = 4/3), the graphical representation of V with 0 shows that V is monotonically 
decreasing function of 0 (Fig 3.1). Also it can be shown that 0 increases with c. In 
case / and II of viscous huid distribution we have obtained solutions for diherent values 
of k and arbitrary parameters. In one of the case, where co-efhcient of bulk viscosity 
r] ~ the potential function V decreases exponentially with cj) (Fig 3.2), but for 

other cases, V cannot be expressed as a function of 0 analytically. In case III, where 
1] ^ p, the relaxation time becomes a constant and ^ < 0 i.e., velocity of light decreases. 
In this case also V cannot be expressed as a function of 0 explicitely, but it can be shown 
that V is monotonically decreasing function of 0. So for all cases of perfect and viscous 
huid distribution the potential function V is monotonically decreasing with (j). 




Chapter 4 

Quintessence Problem and 
Brans-Dicke Theory 

4.1 Prelude 

From theoretical point of view, a lot of works [Ostriker er al, 1995; Peebles, 1984; Wang 
et al, 2000; Caldwell et al, 1998] has been done to solve this quintessence problem and 
possible candidates for Q-matter are Cosmological Constant (or more generally a variable 
cosmological term), a scalar held (which was introduced by Ratra and Peebles [1988]) 
with a potential giving rise to a negative pressure at the present epoch, a dissipative 
huid with an effective negative stress [Cimento et al,2000] and more exotic matter like 
a frustrated network of non-abelian cosmic strings or domain walls [Bucher et al, 1999; 
Battye et al, 1999]. In these works, it is generally assumed that Q-matter behaves as a 
perfect huid with barotropic equation of state and most of them have considered spa¬ 
tially hat model of the universe (only the work of Chimento etal [2000] has been done for 
open model of the universe). Recently, Banerjee etal [2001] have shown that it is possible 
to have an accelerated universe with BD-theory in Friedmann model without any matter. 

In this chapter, the quintessence problem has been discussed in anisotropic cosmo¬ 
logical models for Brans-Dicke theory. We have shown that it is possible to have an 
accelerated expanding universe without any quintessence matter due to the Brans-Dicke 
held. In order to solve the held equations we have assumed power-law form of the scale 
factors. Here we have found that the anisotropic scalar characterizes the accelerating or 
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decelerating nature of the universe. Finally flatness problem has been discussed in this 
context using conformal transformation in Jordan frame. 


4.2 Field Equations and Solutions 


Consider anisotropic space-time model described by the line element 


ds'^ = —dt^ -I- a^dx'^ -|- b'^dQl 


(4.1) 


where a and b are functions of time alone : we note that 

{ dy"^ -\- dz'^, when k = 0 (Bianchi I model) 

d9‘^ + sin^Odcf)^, when k = +1 (Kantowaski-Sachs model) 

dO"^ + sinh'^Odcj)'^, when k = —1 (Bianchi III model) 

Here k is the curvature index of the corresponding 2-space, so that the above three 
types are described by Thorne [1967] as Euclidean, open and semi closed respectively. 

Now, in BD theory, assuming a perfect fluid distribution as only matter held, the held 
equations for the above space-time symmetry are 



1 


(3 -|- 2u;)0 


[(2 -I- u})p^ + 3(1 -I- a;)pj.] — a; 



(4.2) 


M , 9 -- = ^--^ 
b j ab (f) 6 ^ 



0 2 


and the wave equation for BD scalar held is 


(4.3) 


+ (4.4) 

Here and are the density and hydrostatic pressure respectively of the huid 
distribution, obeying the barotropic equation of state 


Pf = hf - 1)P/ 
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( 7 j. being the constant adiabatic index of the fluid causality demands 0 < 7 ^ < 2 ) and 
as usual u is the BD coupling parameter. Now the above held equations, via the Bianchi 
identities lead to the energy conservation equation 

d/ + j {Pf + P/) = 0 (4.5) 

At present, the universe is considered as matter dominated (i.e., hlled with cold matter 
(dust) of negligible pressure) so considering = 0, we have from equation (4.5) (after 
integration) 

P, = y ( 4 - 6 ) 

where V = alP is the volume index at the present instant and is an integration 
constant. Also using = 0 and equation (4.6), the wave equation has a hrst integral 


1/0 


3 + 2o; 


+ Cq 


(4.7) 


We now assume a power-law form of the scale factors keeping in mind that we must 
have an accelerated universe to match the recent observations. So we take 


a{t) = Opt", h{t) = 


(4.8) 


and consequently 

r = {V„ = a^l) (4.9) 

where (ap, 6 p) are positive constants and (a,/3) are real constants. Thus the matter 
density and the BD-scalar held takes the form 

P/= {pf=Po/K) (4-10) 


and 


Pit 


2-a-2P 


{3 + 2uj){2-a-2f3) 


+ 


K(l-«-2/3) 


(4.11) 
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The field equations will be consistent for the above solutions (eqs.(4.8)-(4.11)) provided 
we have the following restrictions on the parameters 


Cq = 0, /? = 1 and 


(4.12) 


either 


or 


a_l, ij_-2(l + — 


a)=-2, j2=a-l 


(4.13) 


(4.14) 


It is to be noted that both the cases coincide for fiat model of the universe and we 
have the solution 


a = Opt, b = b^t, Pj = p^t (j) = p^t 2ci; + 3 = —1 (4-15) 


so for k ^ 0, the solution in first case is 


Qj — ^ 0^1 ^ — ^ 0^1 Pf — Pi^ 1 ^ — 


and 


, 4k \ 

2- + 3 = -|l + 5j^). 


Pit 


-1 


2(n + 3 


(4.16) 


(for fc < 0, 6^ > 4/3) while for the second case the expression for the geometric and 
physical parameters are 


a = a^t, b = bgt, Pj = p^t 


= — r", 2a; + 3 = -1 
a 


(4.17) 


and we have a restriction 


k 


The deceleration parameter has the expression 

' a — 1 


q = 


a + 2 r 


(4.18) 


(Note that a = 1 corresponds to the first two case) 
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Thus we always have an accelerated model of the universe (in the third case i.e., 
eq.(4.17)) (except for —2 < a < 1) as predicted by recent Supernova observation. Hence 
for the solutions represented by equations (4.15) and (4.16) g = 0 i.e., at present the 
universe is in a state of uniform expansion and the conclusion is identical to that of 
Banerjee et al [2001]. Furthermore, the solution corresponds to equation (4.17) is new 
as we have a negative deceleration parameter. This solution is valid for closed (or open) 
model of the universe for a > 1 (or a < 1). Therefore, it is possible to have an ac¬ 
celerated universe today by considering anisotropy model of the universe and hence the 
quintessence problem may be solved for closed, open or flat type of model of the universe. 

4.3 Conformal Transformations : Flatness Problem 

In cosmology, the technique of conformal transformation is often used (for mathematical 
simplihcation) to transform a non-minimally coupled scalar held to a minimally coupled 
one [Faraoni et al, 1999]. Usually, in ‘Jordan Conformal frame’ the scalar held (also 
BD scalar held) couples non-minimally to the ‘Einstein frame’ in which the transformed 
scalar held is minimally coupled. In the last section, we have developed the BD theory in 
Jordan frame and to introduce the Einstein frame we make the following transformations: 

drj = dt, d = a, b = b, tf: = In (j), (4.19) 

As a result the held equations (4.2)-(4.4) transformed to 



(4.20) 



(4.21) 


and 



(4.22) 
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where ' = 

drj 

These equations are the well known fields equations for the anisotropic cosmological 
models (described here) with a minimally coupled scalar field "0 (massless). This scalar 
field behaves like a ‘stiff’ perfect fluid with equation of state 


Pth 




/2 


(4.23) 


IbvrG 

In Einstein frame, the total stress-energy tensor is conserved, but the scalar field and 
normal matter change energy according to 


{pf +Pf) - - i 


+ ( "^ + +P^) 


= -Yfe-3p,) (4.24) 


Thus combining the two energy densities we have from the above equation 


p' -|- S'yHp = 0 


(4.25) 


Here H = ^ is the Hubble parameter in the Einstein frame and 7 is the 

average barotropic index defined by the relation 


yfi = + 7 ^fi^ (4.26) 

where 

n = ^ (4.27) 

is the density parameter. 

From equations (4.21) and (4.25) after some algebra, we have 

n ' = - 1) [777^ + 2(7 - i ) Hh \ (4.28) 

where Ha = ^ and 77^ = |- . 

This evolution in hi shows that hi = 1 is a possible solution of it and for stability of 
this solution, we have 

2 


(4.29) 
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for the solutions given in equations (4.15) and (4.16) and the restriction is 
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7 < 


2 

a + 3’ 


(4.30) 


for the solution (4.17). 

Since the adiabatic indices do not change due to conformal transformation so we take 
7 ^ = 1 (since = 0) and 7 ^ = 2. Hence from (4.26) and (4.27), we have 


7 = — - - 

' 12 ,+ 12 , 

/ -0 

Now due to upper limit of 7 (given above) we must have the inequalities 

7 <4171 

OR 


(4.31) 


2 (q^ (4.32) 

■' a ^ 

according as 7 is restricted by (4.29) or (4.30). Also from the held equation (4.21), 
the curvature parameter 12 = —k/l? vanishes for the solution 12 = 1 , provided we are 
restricted to a = 1. Therefore, depending on the relative magnitude of the energies 
of matter and the BD-scalar held (as in isotropic case) it is possible to have a stable 
solution corresponding to 12 = 1 and hence the hatness problem can be solved. 


4.4 Discussion 

In this chapter, we have considered three anisotropic cosmological models namely, Bianchi 
III {k < 0), axially symmetric Bianchi I (A; = 0) and Kantowski-Sachs {k > 0) space- 
time. We have shown that the anisotropic character is responsible for getting an accel¬ 
erated model of the universe. In fact, for the present model the shear scalar is given 
by = |(q! — 1)^, so the deceleration parameter (eq. (4.18)) is proportional to 
and it shows how anisotropy characterizes the accelerating or decelerating universe. In 
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other words, we can say that anisotropic nature of the universe has an effect on the 
quintessence problem. The problem for negative coupling constant uj is same as in 
isotropic case and there is problem in big-bang nucleosynthesis scenario as claimed by 
Banerjee etal [2001]. In this case the modihed version of BD-theory (where the coupling 
parameter a; is a function of the scalar held) is very similar to that for isotropic case, so 
we have not presented in this chapter. 


Chapter 5 


A Quintessence Problem in 
Self-interacting Brans-Dicke Theory 

5.1 Prelude 

The quintessence proposal faces two types of problems [Steinhardt et al, 1999]. One 
of these problems (referred to as hne tuning problem, eliminate classically by Ratra 
and Peebles [1988]), is the smallness of the energy density compared to other typical 
particle physics scales. The other problem known as the cosmic coincidence problem is 
that although the missing energy density and matter density decrease at different rates 
as the universe expands, it appears that the initial condition has to be set so precisely 
that the two densities become comparable today. Quintessence has been proposed as 
that missing energy density component that along with the matter and baryonic density 
makes the density parameter equal to 1. A special form of quintessence held called the 
‘tracker held’ has been proposed by Ratra and Peebles [1988] to tackle this problem 
[Also Zlatev et al, 1999]. Since in a variety of inhationary models scalar helds have 
been used in describing the transition from the qnasi-exponential expansion of the early 
nniverse to a power law expansion, it is natural to try to understand the present accel¬ 
eration of the universe which has an exponential behaviour too, by constructing models 
where the matter responsible for snch behaviour is also represented by a scalar held 
[Strarobinsky, 1998; Saini et al, 2000]. Inverse power law is the other potential [Peebles 
et al, 1998; Ferreira et al, 1987] that has been studied extensively for qnintessence mod¬ 
els, particularly, for solving the cosmic coincidence problem. Recently, Bertolami and 
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Martins [2000] obtained an accelerated expansion for the nniverse in a modified Brans- 
Dicke (BD) theory by introdncing a potential which is a fnnction of BD scalar held itself. 

Similar to the previous chapter, we have studied quintessence problem in this chapter 
for Brans-Dicke theory but with a scalar held which is self-interacting. Assuming the 
power law form for the scale factors and Brans-Dicke scalar held, we have studied the 
accelerating or decelerating nature of the universe for various values of the parameters 
involved. Though for the present universe we have considered dust models only but it 
has been shown in principle that perfect huid with barotropic equation of state can be 
considered to study the quintessence problem. Lastly, in addition to the non-decelerating 
solution, hatness problem has been solved potentially. 


5.2 Field Equations and Solutions 

The Brans-Dicke theory is given by the action [Brans and Dicke, 1961] 


S 




cf)R — — 0 -|- Cm 

(P 


d^x 


(5.1) 


where 0 is the BD scalar held, uj is the dimensionless constant BD parameter and Cm 
is the Lagrangian for all other matter helds. We have chosen the units SttG = c = 1. 

The matter content of the universe is composed of perfect huid and a scalar held as 
the quintessence matter. We assume that the universe is homogeneous and we consider 
an anisotropic space-time with line-element which is given in eq.(4.1). 

Now the BD-held equations are 


b yb j (p 2^ V / ^ 


a b ah {Pr^+P^) 1 

a b a b cp 2 



(5.2) 


(5.3) 






CHAPTER\^ A Quintessence Problem in Self-interacting Brans-Dicke Theory 


54 


, ^ab ^ k _ip^+ p^) (a , 

h] ^ ah^ ~ 6 \a^ b 


and the wave equation for the BD scalar held 0 is 


0 UJ 


(6.4) 




[(P„ - SpJ + (Pt -3p,)] 


(5.5) 


p^ and are the density and the pressure of normal matter, p^ and are those due 
to the quintessence held given by 

= \'>P^ + V (^), ^ 

where V (0) is the relevant potential. 

The wave equation for the quintessence scalar held 0 is 


7 ; dV{ijj) 

From the above held equations, we have the matter conservation equation 


(6.7) 


Assuming that at the present epoch, the universe is hlled with cold matter (dust) with 
negligible pressure, so using pm = 0, the conservation equation (5.8) gives 


pm — 


Pi 

ab'^ 


(5.9) 


where pi is an integration constant. 

Now we assume, the power law form of scale factors a, b and the BD scalar held 0 are 


= aif“, b = bit^, 0 = (fit^ 


( 6 . 10 ) 


where Oi, bi, 0i are positive constants and a, 0, 6 are real constants with a + 20 > 3 
(for accelerating universe). 
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From the field equations (5.2), (5.3) and (5.4) using (5.10) we have the expression for 
as 

^2 ^ ^5-213 — P^f-a-2i3— 2a{a — l) —ujS'^ — S{6 — 1) — {a — 2/3)6]t^~^ (5.11) 

U-^ CI\U-^ 

From (5.5) the potential V is given by 

^ ~ + 3)(a + 2/9 + (5 - 1)5 + 2/9^ - 2Q;(a; - 1) 

2b-^ lQj\b-^ 4 

-0^5^ - 5(5 - 1) - (a - 2j3)5]t^-^ (5.12) 

The wave equation (5.7) for the quintessence scalar field if can be written in the form 

After integration (5.13) the expression for V is 

T/_ ^01 (2a + 2/? + 5) ^ 5 _ 2/3 Pi ^-a- 2 p , 0i(2a + 4/? + 5 - 2) 

bf (b- 2/3) 2a,6? + 2(S- 2) ' ' ' 

+5(5 - 1) + ojS^ + (a- 2/9)5 - 2/9^] (5.14) 

Now, the consistency relations of the constants coming from the two identical equations 
(5.12) and (5.14) for the potential V are 

4a+ 2/9+ 35 = 0 (5.15) 

5(5 - 2)(a + 2/9 + 5 - l)(2o; + 3) = (4a + 8/9 + 35 - 6)[2a(a - 1) + 5(5 - 1) 

+ud^ + (a - 2/9)5 - 2/9^] (5.16) 

From these relations for an accelerating universe (a + 2/9 > 3), there are two possibil¬ 
ities: one in which k = 0 and the second where k ^ 0. 
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Case I: k = 0 

In this case, consistency conditions are (5.16) and 

a + 2f3 = c (5.17) 

where c is an constant (> 3). 

For solving (5.16) and (5.17), we may choose 6 = —2 and we have two possible 
solutions: 

(i) a = f3=l 6 =-2 

or 

(ii) a = 4 — c, /3 = —2 + c, 6 = —2 

For both solutions (i) and (ii) the expression for (eq.(5.11)) becomes 





2(2o^ + 3)0it-^ 


This indicates that a; < —3/2 as -0^ cannot be negative. 

For c = 4, one has 2|2u; + 3|0i > and the deceleration parameter q 
In this case, equation (5.18) integrates to 





(5.18) 


1 

4' 


(5.19) 


where 


— 2{2ijj T 3)01 — 



and the relation between V and 0 becomes 
V = (5.20) 


where Vi being a constant, related to the constants ai,pi etc. The model works for 
all time 0 < t < cx) with the condition 2|2a; + 3|0i > is satished. The value of u is 
related to the other constants as follows: 


2(2a; + 3)0i H- pr 

aibi 



for solution (i) 


2(2a; + 3)0i 




0,1), for solution (ii) 


and 
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For c 7 ^ 4, the model does not work for the whole range of time 0 < t < cx). In this 
case, the deceleration parameter is g = | — 1. 

If c > 4, then the rate of acceleration is faster than q = —^ and from equation (5.18), 
"0^ > 0 restricts the validity of the model for t > ti where 


t = 

^ 2\2uj + 3|ai6f0i 


1 

c—4 


( 6 . 21 ) 


Further if 3 < c < 4, then the universe will expands with an acceleration but with a 
rate less than q = —| and as before from equation (5.18) for real the model works 
upto the time ^2 where 

Pi 


h — 


2\2uj + 3|ai&^0i 


1 

c—4 


(5.22) 


For q = — 2 , the present age of the universe can be calculated from (5.4) as 


to — 


2 - 20. - — + ^ 

201 01 


1/2 




where H„ = - and Hh = t- 


For large ui limit. 


to = 


\/—2lo 


[{H^)o + 2{HMHb)o]"^" 

where oj is obviously a negative quantity. 


Case II : k ^ 0. 

In this case, consistency conditions are (5.15), (5.16) and 

a + 2/3 = 3 


Solving these three equations we have the following solutions: 
(i) a = 1, P = 1, 6 = —2 for all values of oo. 


ii) ^ = = 5 = -2± 


2uj+3 


provided for oj > 3/2 or a; < —3/2. 


(5.23) 


2ai-3 
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For the solution (i), the model works for a limited period of time, 0 < t where 


ti = 201 


k 


- {2uj + 3) 


ai^i 

Pi 


(5.24) 


For an open universe i.e., for k = —1, {2uj + 3) < 0 and |2a; + 3| > ^. 

For a closed universe i.e., for k = +1, {2uj + 3) > 0 and |2a; + 3| < ^. 

For the solution (ii) the model works for the time, where t satishes the equation 


2{2uj + 3)01 


It 


2a; + 3 


2u;-3 2a; - 3 


T2 


(6.25) 

Qj^O-^ 0-^ 


provided for uj > 3/2 or a; < —3/2. 


5.3 Flatness Problems and its Solutions 

One important aspect of this model is that potentially it can solve the flatness problem. 
Now we make a conformal transformation [Faraoni et al, 1999] as 




(5.26) 


In this section, we make the following transformations: 


dr] = y/^dt, a = y/0a, b = b, 0 = In 0, = 0 p = 


-2 

10 'f' P 'll} 1 


Pn,=<l> Vm. \ 


(5.27) 


So the held equations (5.2) - (5.4) transformed to (after some manipulations) 


b'\^ a" a! y k _ (3 + 2a;) 00'^^ 

^1 —2-2—^ ^ = (p +P,)H- 

b a ab b^ ^ 


(5.28) 


b'\ a'b' k , , (3 + 2a;) /0' 


(5.29) 
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and 


h" a" k _a'h' fb' 


where ' = 

dr) 


The BD scalar field in new version is given by 


(3 + 210 ) f (j)' 

p^ = —^—b:) 


We define the dimensionless density parameter hi as 


(5.30) 


(5.31) 


a = n +si., + si, = -E- 

m ' </. ' 4 , o j^2 


(5.32) 


where p = is the total density and hlj are defined accordingly. 

Using (5.28)-(5.31), and combining the energy densities, we have the equation for the 
conservation for the total energy, 


p + ?>H{p + p) = 0 


(5.33) 


Here H = | -\- 2~j is the Hubble parameter in the Einstein frame and 7 is the net 

barotropic index defined as 


7 SI = 7 „ Sl„ + 7 , $1, + 7 „ SI, 


(6.34) 


From equations (5.29) and (5.33), we have the evolution equation for the density 
parameter as 

- 1) [-fHa + 2(7 - l)Hb\ (5.35) 

where - and = I • 

“a ° b 

The individual 7 ^ ’s are defined by the relation pi = ( 7 * — l)pj. So the ratios ^ remain 
same in both frames. For our choices of matter, = 0 and = p^, so we have 7 ^ = 1 
and 7 ^ = 2. The other index 7 ,^ is related by the equation 


+Pv, ^ 

P 4 , ^ 


(5.36) 
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It has been shown that 7 ^ is varies with time. The equation (5.35) indicates that 
= 1 is a possible solution and this solution determines that () _ <0. This solution 


H 


is stable for expanding universe {H > 0) if 7 < | with the relevant condition 


+ 4S2* < (2 - (5.37) 

From the held equation (5.29), the curvature parameter = —kj I? vanishes for the 
solution = 1. So for BD-scalar held it is possible to have a stable solution correspond¬ 
ing to = 1 and hence the hatness problem can be solved. 


5.4 Discussion 

We have investigated the nature of the potential relevant to the power law expansion 
of the universe in a self-interacting Brans-Dicke (BD) cosmology with a perfect huid 
distribution for anisotropic cosmological models. We have considered a non-gravitational 
quintessence scalar held if with a potential V = V{ip). This scalar held in BD-theory 
is shown to give rise to an accelerated expansion for the present dust universe (where 
we have taken = 0). It is to be noted that at early stages of the evolution of the 
universe pm is non-zero. But if we take barotropic equation of state Pm = (7 — l)Pm, 
then equation (5.9) is modihed to pm = 

We have presented a class of solutions describing non-decelerating universe for both 
hat {k = 0) and curved space-time (/c 7 ^ 0). For = 0, there are two possible solutions 
for diherent choice of the parameters. In both the solutions, the parameter uj must be 
negative (in fact 2a; -|- 3 is negative) to make the quintessence held real. The validity 
(on the time scale) of the solutions depends on the parameter c (dehned in eq.(5.17)). 
For c = 4, the model works for all time 0 < f < cx), while for c > 4, the model works for 
{ti, 00 ) where ti is given by (5.21). From this model we cannot predict the geometry of 
the universe before ti. Similarly for 3 < c < 4 the model is valid for t < t 2 i.e., from 
the begining to the time t 2 - Here also we cannot predict the nature of space-time of the 
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universe after t = t 2 - 

In curved space-time {k ^ 0) there are also two possible solutions for different choice 
of the parameters. The time interval over which the solutions are valid are given in 
equations (5.24) and (5.25). For close model {k > 0) the coupling parameter uj is 
restricted by the inequality 0<2tn-|-3<l/6f while for open model oj satishes —l/h\ < 
2a; -|- 3 < 0. From the above solutions, we note that the coupling parameter uj may be 
positive or negative (with some restrictions). Hence for all solutions the choice of uj is 
not in agreement with the observations. 

Finally, for non-decelerating solution it can also potentially solve the flatness problem 
(without any restriction on the parameters) and it has been shown that 0 = 1 could be 
a stable solution in this model. 


Chapter 6 

A Quintessence Problem in 
Brans-Dicke Theory with Varying 
Speed of Light 

6.1 Prelude 

In this chapter, quintessence problem has been studied as in the last two chapters for 
the Brans-Dicke theory but considering time variation of the speed of light. Here also 
without any quintessence matter, an accelerated form of the universe is possible. For 
power law expansion of the universe, we have shown that the time variation of the ve¬ 
locity of light is also in the power law form. Further assuming exponential form for 
the scale factors and the Brans-Dicke scalar held, the velocity of light is also in the 
exponential form and the deceleration is negative dehnite without any dependence of 
the parameters involved. Similar to the previous chapters the hatness problem can be 
solved automatically. 


6.2 Field Equations and Solutions 

We consider the line-element of anisotropic space-time model 

ds"^ = —c^dt^ + a^dx"^ + d^d^ll ( 6 . 1 ) 


where a, b are functions of time only and dQ], is given in chapter 4. 
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Now, the BD-field equations with varying speed of light are 


a Stt 

a 6 (3 + 2u})(j) - 


(2 + u!)p, + 3(1 + O’)—^ 


— a; 





( 6 . 2 ) 


kc^ I 9 ^^ I ^ 

b) ab 0 62 ya bj (j) 2 ' ' 

and the wave equation for the BD scalar held (j) is 


(6.3) 


, a ^b\ ■ Svr / 


3P, 


(6.4) 


Here the velocity of light c is an arbitrary function of time, u is the BD coupling 
parameter, p^ and are density and hydrostatic pressure respectively of the huid 
distribution with barotropic equation of state 


Pf = ilf - 1)P/ 

( 7 ^ being the constant adiabatic index of the huid with 0 < 7 ^ < 2 ). 
From the above held equations, we have the ‘non-conservation’ equation 


Pf 


a b 

+ 2-, 
a b 


, Pf 


kcc 

Anh^ 


( 6 . 6 ) 


As at the present epoch the universe is hlled with cold matter (dust) with negligible 
pressure, so using = 0, the equation (6.4) and (6.5) gives 


and 




Svr 

3 2a; 



d 

dt 


{ab^Pj) 


kacj) d 2 
Svr dt 


( 6 . 6 ) 


(6.7) 


Now eliminating p^ between ( 6 . 6 ) and (6.7), we have 
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Since the scale factors, scalar field and velocity of light are time dependent, so we 
shall consider the following cases to obtain exact analytic form for the variables assum¬ 
ing some of them in polynomial form or in exponential form. 

Case I : The power law form of scale factors a, b and BD scalar held 0 are assumed as 

a{t) = b{t) = 0(t) = (6.9) 

where are positive constants and a, /3, /i are real constants with the restriction 

+ 2/9 > 3 (for accelerating universe). 

If we use the relations (6.9) in (6.8), we have 

c = (6.10) 

where the positive constant is restricted as 

2kC(3 — 11 

b‘^fi{a + 2/3 + — l){a + 2(3 + — 2) = —— (6-11) 

Also from (6.6) the expression for is 

P, = (6.12) 

with 

3 I 2co^ 

Po = 0oP(« + 2^ + P - 1) 

For consistency of the held equations the restriction between the parameters is 
{^ + 2(3 + p - ^ =— - ^ -^[(2-t-a;)(l-/?) - (a-/r)(3-F 2a;)]- (6.13) 


Case II : The exponential form of scale factors and BD scalar held are assumed as 

a{t) = OgC"*, b{t) = 0(t) = (p^e^^ 


(6.14) 
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where are positive constants and a,l3,pL are real constants. 

From (6.8), we have 

c = 

where the positive constant Cp will satisfy 

AikScP' 

(3 + 2a’)yu(cr + 2/3 + /i)^ = — ° 

^0 

From (6.6), the expression for is 

Pf = 

with Po = + 2/3 + p). 

From held eqnations the relation between the parameters becomes 

either a + 2/3 + /r = 0 


(6.15) 


(6.16) 


(6.17) 


(6.18) 


or (a + 2/3 + /i) 11 + | = /^(4 + 3a;) (6.19) 

In both the cases we have the deceleration parameter 

g = —1 < 0 

Case III : In this case, we have assnmed the power law form of scale factor b, velocity 
of light c and BD scalar held 0 to be 

b{t) = b^t^, c{t) = (j){t) = (j)^C (6.20) 

Using (6.20), we have from (6.8), 

( 6 . 21 ) 

To solve the diherential eqnation we assnme 6 = (3. It is to be noted that the above 
diherential eqnation can be solved for all k. Bnt for k = ±1 the solntions are not 
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Fig 6.2 



Fig 6.3 
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In Figs 6.1 - 6.5, we have shown the variations of q over t for various values of co and the values 
of parameters oq = ui = 1. We have taken w = —1.8,—1.5,—1,—0.87,—0.75 respectively in 
Figs 6.1 - 6.5. 
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Fig 6.6 


In Fig 6.6, we have shown the variations of q over t and uj in the range 0 < t < 10 and 
-2.5 < < 0. 


consistent to the field eqnations. So we consider only k = 0. In this case the explicit 
solntion is 

_ ^+2 2(a; + l) ' 

a{t) = a^t + a^t 


2(a; + l) 

b{t) = b^t 3“+4 

2(a; + l) 

C{t) = Cgt 31^+4 


0(t) = 


^ ^ 2a;+3 

^f Alt 3a;+4 (aQ+a-^t) 


with ap,a^ as integration constants. 

The deceleration parameter has the expression 


( 6 . 22 ) 


3[(3cj+2)(aQ+ccj (3a;+4)£^] 

^ [(3a^+2)(aQ t)+a^ {3a2+4)t]^ ’ 

which is finite for all t and ui except t = — and we have the non-decelerated 

nniverse for —2 < uj < —2/?). Also we have the non-decelerated nniverse asymptotically 
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except for —1 < a; < —1/2. The variation of q over time has been shown graphically for 
different values of oj in figures 6.1 - 6.5. 


6.3 Conformal Transformation and Flatness Problem 

One important aspect of this model is that it can solve the flatness problem as well. To 
see this we make a conformal transformation [Faraoni et al, 1999] as 

9uv ~ (6.23) 

which enables us to identify the energy contributions from different components of 
matter very clearly. 

In this section, we have developed the BD theory in Jordan frame and to introduce 
the Einstein frame, we make the following transformations : 

dr] = \/0 a, d = a, b = b, = In 0, 


Pf = (^ P^=(^ \ 


(6.24) 


So the held equations (6.2) - (6.4) transformed to 


d" ^b” . f (3 +2a;) ^,2 


2j^-iT[p,+ 


(6.25) 


b'V _a'6' 


a' a He" „ _ (3 + 2a;) ,2 


b J abb 4 




and 


d' J 
a b 


Stt 




where ' = 

drj 


(6.26) 


(6.27) 


The scalar held ip (massless) behaves like a ‘stih’ perfect huid with equation of state 


P^ Pjb 


IQnG 


(6.28) 
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If the velocity of light is constant, then in Einstein frame total stress-energy tensor is 
conserved but there is an exchange of energy between the scalar held and normal matter 
according to the following equation 




dv- + ( -^ + 2 j ) (Pv- +p^) 


f _ 3p 
Pf-^ 


(6.29) 


On the other hand, if the velocity of light varies then we have separate ‘non-conservation’ 
equations 


a' y 


Pf 


^/ + ^ + + 7^ - “T + 


fj' 3p 


kcc' 

AttGP 


and 


a' y 


3p, 


+ ^ + 2 t {p^ + P^) - ^ { Pi. - 


(6.30) 


(6.31) 


Thus combining the two energy densities, we have from the above equations, the 
equation for the conservation for the total energy is 


p' + 3-fHp = 0 (6.32) 

Here H = | + 2—'^ is the Hubble parameter in the Einstein frame and 7 is the net 

barotropic index dehned as 

7 O = 7 ^ 0 ^ + 7 ^ 0 ^ (6.33) 

where 

LI = (6.34) 

is the dimensionless density parameter. 

From equations (6.26) and (6.32), we have the evolution equation for the density 
parameter as 

f2' = n(n-l)[7if, +2(7-1)^6] 

where Ha = ^ and Hh = t ■ 

^ a ^ 0 


(6.35) 
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This equation in ^2 shows that ^2 
this solution, we have 


1 is a possible solution of it and for stability of 



(6.36) 


Since the adiabatic indices do not change due to conformal transformation so we take 
7 ^ = 1 (since = 0) and 7 ^ = 2. Hence from (6.33) and (6.34), we have 


+ 2 ^ 2 ,/, 

7 = 

/ '0 

Now due to upper limit of 7 , we must have the inequality 

<4|nj 


(6.37) 


(6.38) 


according as 7 is restricted by (6.36). 

From the field equation (6.26), the curvature parameter = —kc^/ P vanishes for 
the solution (2 = 1. So for BD-scalar field it is possible to have a stable solution corre¬ 
sponding to (2 = 1 and hence the flatness problem can be solved. 


6.4 Discussion 

We have performed an extensive analysis of solutions to Brans-Dicke theories with vary¬ 
ing speed of light. For the power-law forms and exponential forms of the cosmological 
scale factors and scalar held in Cases I and II respectively, we have the velocity of light 
in the same form and we identihed the cases where the quintessence problem can be 
solved for some restrictions on parameters. In the Case //, we always have the acceler¬ 
ated universe. For Figs 6.1 and 6.2, q decreases to a hxed negative value asymptotically 
for uj = —1.8 and uj = —1.5 . For Fig 6.3, q is linear with time and always negative 
for UJ = —1 . For Fig 6.4, q decreases and then there is a singularity and after that 
q increases. In Fig 6.5, q increases to a hxed positive value asymptotically. From Fig 
6 . 6 , we conclude that the decelerated parameter is negative whenever —2 < a; < —2/3 . 
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In this range q decreases till t = and then increases to a hxed positive value 

asymptotically. 

Along with providing a non-decelerating solution, it can solve the flatness problem 
also. In fact, it has been shown that 0 = 1 could be a stable solution in this model. 



Chapter 7 


Does Cosmic No-Hair Conjecture in 
Brane Scenarios Follow from 
General Relativity? 

7.1 Prelude 

The idea of brane world scenarios may resolve the challenging problem in theoretical 
physics namely the unihcation of all forces and particles in nature. It is suggested that 
we live in a four dimensional brane embedded in a higher dimensional space-time. As a 
result, the fundamental higher dimensional Planck mass could be of same order as the 
electro weak scale and thereby one of the hierarchy problems in the current standard 
model of high-energy physics are resolved [Maartens, 2001; Chakraborty et al, 2002; 
Randall et al, 1999]. 

According to Randall and Sundrum [1999] it is possible to have a single massless 
bound state conhned to a domain wall or 3-brane in hve dimensional non-factorizable 
geometries. They have shown that this bound state corresponds to the zero mode of 
the Kaluza-Klein dimensional reduction and is related to the four dimensional gravita¬ 
tion [Chakraborty et al, 2002]. Hence all matter and gauge helds (except gravity) are 
conhned to a 3-brane embedded in a hve dimensional space-time (bulk) while gravity 
can propagate in the bulk. As a consequence, the gravity on the brane can be described 
by the Einstein’s equations modihed by two additional terms, namely (i) quadratic in 



CHAPTER^ Does Cosmic No-Hair Conjecture in Brane Scenarios Follow from General Relativity?! A 

matter variables and (ii) the electric part of the five dimensional Weyl tensor [Maartens, 

2000; Campos et al, 2001]. 


In general terminology, the CNHC [Gibbons et al, 1977; Hawking et al, 1982] states 
that “a// expanding universe models with a positive cosmological constant asymptotically 
approach the de-Sitter solution”. To address the question whether the universe evolves 
to a homogeneous and isotropic state during an inflationary epoch. Gibbons and Hawk¬ 
ing [1977] and then Hawking and Moss [1982] developed this conjecture. Subsequently, 
Wald [1983] gave a formal proof of it for homogeneous cosmological models (Bianchi 
models) with a positive cosmological constant. He assumed that the matter held should 
satisfy strong and weak energy conditions. 

In this chapter, we wish to extend Wald’s [1983] result in the brane world scenario and 
examine whether the new conditions can be minimize using those in general relativity. 

7.2 Cosmic No Hair Conjecture in Brane World 

According to Roy Maartens [2001], the Einstein equations on the brane can be written 
as 



( 7 . 1 ) 


where and are the two correction terms (local and non-local) in the energy 
momentum tensor. The local correction term has the expression 



( 7 . 2 ) 


while is the electric part of the 5D Weyl tensor in the bulk. Now the scalar 
constraint (initial value constraint) equation and the Roychoudhuri equation on the 
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brane has the form 




— A + k/T. 


4-*- ^ 






n^rP 




(7.3) 


and 

= —A + 

where is the unit normal to the spatial homogeneous hypersurfaces. In terms of 
the homogeneous hypersurface elements namely, the projected metric h^i/(= g^u + n^riu) 
and the extrinsic curvature VyU^) and using the Gauss-Codazzi equations the 

above two equations namely equations (7.3) and (7.4) become 

K^ = ZK + + ?>KlT^yn^^rC + Z4S^yn^^rP - SSf.yU^^n'' (7.5) 


T/^y - -Qf^yT n^rP + S^y - -g^,yS n^rp - S^^yU^rP (7.4) 


and 


k 


- a- k\ 




n^rC 





n^rC + S^y-n^rR 

(7.6) 


where the dot denotes the Lie derivative with respect to proper time, K is the trace of 
the extrinsic curvature, is the shear of the time like geodesic congruence orthogonal 
to the homogeneous hypersurfaces and ^^'>R is the scalar curvature of the homogeneous 
hypersurfaces. 

Using the idea of Wald and proceeding along his approach (for details, Wald et al 
[1983] and Chakraborty et al [2001]) one can hnd that for CNHC we must have 

(a) S^yU^^rN > 0 and > 0 (7.7) 

and 


( 6 ) 


E^yU^rP < 0 


(7.8) 


in addition to the weak and strong energy conditions for the matter held 
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Now if we use the expression (7.2) for S^i, in (7.7) then we get 


^-Th - - (T^X)(^X) > 0 (7-10) 

and 

^Ta-(r^X)(7^X)>0 (7.11) 

where a = and b = (T^u — nPn’' are positive due to (7.9). 

Also using the symmetry properties of it is possible to decompose it with respect 
to any time like observer u {u'^Ua = — 1) as [Maartens, 2000; Campos et al, 2001] 


S = - 


4 r 


) If + 2(m^Qi/) + 


with the properties 


Q^uP = 0 , PP = 0 , Pf.^u’' = 0 


If we consider the Bianchi models then due to the symmetry of the spatial geometry 
we may choose 

Q^i = Pfiu = 0 , 

and the scalar part namely U is termed as dark energy density as it 
momentum tensor that of a radiation perfect fluid. So the restriction (7.8) 
dark energy density should be always positive i.e., 

W > 0 (7.12) 


has energy- 
implies that 


As it is not possible to make any restriction on to satisfy inequations (7.10) and 
(7.11), so let us examine with some realistic model for the matter held. 
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7.3 Examples 

7.3.1 Perfect Fluid Model 

In this case the energy-momentum tensor has the form 

T^u = {p + p)n^n^ -f pgf,u, = -1 

with p and p as the energy density and isotropic pressure respectively. 

The weak and strong energy conditions demand 

a = p>0 and 26 = p-|-3p > 0 (7-13) 

Hence the inequations (7.7) (i.e., inequations (7.10) and (7.11)) take the form 

> 0 and p(3p-|- 2p) > 0 (7.14) 

which are always true. Thus for perfect fluid model CNHC is automatically satisfled 
in brane scenarios if it is valid in general relativity. 

7.3.2 General Form of Energy-Momentum Tensor 

The general form of the brane energy momentum tensor for any matter flelds (scalar 
held, perfect fluids, kinematic gases, dissipative fluids, etc.) including a combination of 
different flelds can be covariently written as [Maartens, 2001] 

= pn^Hu + ph^y + (7.15) 

Here the energy flux and the anisotropic stress are projected, symmetric and 
traceless that is 

= 0, = 0, = H^,^, = 0 

Thus for this form of energy-momentum tensor, the restrictions on now result 

\p" - > 0 


(7.16) 
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and 

lp{3p + 2p)-qX>0 (7.17) 

As for realistic matter, the energy density should be larger than the anisotropic stress 
and heat flux is very small in magnitude so the inequalities (7.16) and (7.17) are auto¬ 
matically satisfied. Hence the CNHC is satisfied for the above form of general energy- 
momentum tensor. 
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Part B 


Study of Some Collapsing Models 



Chapter 8 

General Introduction 


An outstanding problem in relativistic astrophysics and gravitation theory today is the 
hnal outcome of an endless gravitational collapse. What will be the end state of such a 
continual collapse which is entirely dominated by the force of gravity? According to the 
general theory of relativity, which is a classical gravity theory, after undergoing super¬ 
nova explosion in the last state of its evolution, a star with a mass much larger than the 
solar mass will contracted without limit, due to its strong gravity and form a ‘domain’ 
called space-time singularity. Space-time singularity is a very general phenomena, not 
only the gravitational collapse of stars of very large but also in physical processes in 
which the general theory of relativity plays an important role. In fact it was proved by 
Hawking and Penrose [1965, 1967, 1970] that the appearance of space-time singularity 
is generic i.e., space-time singularity appears for any space-time symmetry. However the 
singularity theorems of Hawking and Penrose only prove the causally geodesic incom¬ 
pleteness of space-time and say nothing about the detailed features of the singularities 
themselves. For example, one can not obtain information about how the space-time cur¬ 
vature and the energy density diverge in a space-time singularity from these theorems. 

Space-time singularities can be classihed into two kinds, according to whether or not 
they can observed. A singularity in space-time is said to be naked singularity when non¬ 
space like geodesics can escape from it in future direction that means the singularity is 
visible to the local or distant observer. If such singularity can reach the neighbouring 
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or asymptotic regions of space-time, the singularity becomes locally or globally naked 
respectively. A space-time singularity that can not observed is called black hole. Is such 
naked singularity formed in our universe? With regard to this question, Penrose [1969, 
1979] proposed the so-called Cosmic Censorship Conjecture. There are two versions of 
this conjecture. The Weak Cosmic Censorship Conjecture states that the singularities 
produced by gravitational collapse must always hidden in a black hole [Penrose, 1969]. 
If weak cosmic censorship conjecture fails then gravitational collapse can (generically) 
result in a naked singularity rather than black hole. If so, then it is quite possible that 
the formation of black hole would be a non-generic outcome of collapse . This conjecture 
implies the future predictability of the space-time outside the event horizon. A closely 
related conjecture, known as the Strong Cosmic Censorship Conjecture [Penrose, 1979], 
asserts that, generically, time-like singularities never occur, so that even an observer who 
falls into a black hole will never experience the singularity. It states that all physically 
reasonable space-times are globally hyperbolic. 

Most of the counter examples of Cosmic Censorship Conjecture present in the liter¬ 
ature either belong to Vaidya space-time [Vaidya, 1943, 1951] or to the Tolman-Bondi 
space-time [Tolman, 1934; Bondi, 1947]. The Vaidya metric describes the external 
gravitational held of a radiating star, whereas the Tolman-Bondi metric describes the 
gravitational held inside a spherically symmetric inhomogeneous distribution of dust 
matter. It has been shown by Lemos [1992] that the nature of naked singularities are 
same in Vaidya and Tolman-Bondi space-times. When the binding energy function in 
Tolman-Bondi metric goes to inhnity, the Vaidya metric and Tolman-Bondi metric rep¬ 
resent the same collapse. 


CHAPTER\^ General Introduction 


82 


8.1 Spherical collapse 

Before restricting matter fields, one can present the Einstein equation for a general 
spherically symmetric space-time. In spherically symmetric space-time, without loss of 
generality, the line element can be written as 


ds^ = ^ g2A(i,r)^^2 ^ gin^Qd^)'^) (8.1) 


One can choose a reference system which is comoving for matter helds. In this co¬ 
ordinate system the stress-energy tensor Tj^ that is the source of spherically symmetric 
gravitational held must be of the form 


TU — 


/ -p 0 0 

0 S 0 

0 0 n 

y 0 0 0 


0 \ 
0 
0 

ny 


( 8 . 2 ) 


where p(t,r), S(t,r) and n(t, r) are the energy density, radial stress and tangential 
stress respectively. If one considers a perfect huid, which is described by 


= (P + p)u^u^ + pg^u (8.3) 

then the stress is isotropic i.e., S = 11 = p. 

From the Einstein equation and the equation of motion for the matter helds, one can 
obtain [Harada et al, 2002] (using units, SnG = c = 1) 


m' = AirpR^R' 

(8.4) 

m = —AnTiR^R 

(8.5) 

R! = Ru + R'X 

(8.6) 

, R' 


(p + s)i.'-2 (s-n) - 

(8.7) 
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m=]^R{l- (8.8) 

where m = m(t, r) is the Misner-Sharp mass [Misner et al, 1964] and the prime and 
dot denote partial derivatives w.r.t. r and t respectively. 


8.1.1 Dust Collapse 

Tolman-Bondi-Lemaitre (TBL) [Tolman, 1934; Bondi, 1947; Lemaitre, 1933] class of 
solntions of Einstein’s eqnations represents inhomogeneous spherically symmetric dust 
models. One can restrict the matter content of the model to a dust fluid i.e., a pressure 
less fluid. Therefore E = If = 0. So from equations (8.4) to (8.8), one obtains the 
solutions 


, , F(r) 

m(r) = 2 

(8.9) 

1 / = 0 

(8.10) 

p/2 

2A _ ^ 

1 + /(r) 

(8.11) 

= /(r) + ltd 

(8.12) 

/ , F' 

(8.13) 

- Rm' 


where the arbitrary functions F(r) and f{r){> —1) are twice the conserved Misner- 
Sharp mass and the specific energy respectively. In equation (8.10) one can use the 
rescaling freedom of the time co-ordinate. This means that a synchronous comoving 
co-ordinate system is possible. 

For gravitational collapse of dust model, R(t, r) < 0 and without loss of generality 
rescale R such that i?(0,r) = r. So the solution of R (from eq. (8.12)) has the implicit 
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form [Dwivedi et al, 1992] 


t = 


v^L 


^3/2^ f _4 




(-¥)] 


(8.14) 


Here G{y) is a strictly real positive and bounded function which has the range —oo< 
y < I and is given by 


G{y) = 


^Vv 


2 

3 




for 0 < y < I 


for y = 0 


sink ^yj—y 


yi-i/ 

y 


for ?/ < 0 


(8.15) 


One should note that the cases / < 0, / = 0 and / > 0 correspond respectively to 
the ranges in terms of y given by0 <r/<l, |/ = 0 and |/ < 0. Also it should be noted 
that / < 0, / = 0 and / > 0 represents bound (elliptic), marginally bound (parabolic) 
and unbound (hyperbolic) models. 


The solution for marginally bound collapse (/ = 0) is explicitly given by 


^ {(„(r) - 


(8.16) 


where ts/(r) is given by 




^ ^3/2 

3v/F 


(8.17) 


which represents the shell focusing singularity [R = 0) time and the singularity occurs 
at different instants of time for different shells characterized by r-constant hypersurface. 
The central shell focusing singularity occurs at r = 0 and the corresponding time being 
^0 = ^s/(0). This brings in the essential difference from the analogous situation in the 
homogeneous space-time where the function tsf{r) is a constant and all the particles 
reach the centre simultaneously. 
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The total mass within radius r is given by the integral 

rr nr 

m{r) =4:71 pR^dR = dvr / F'{r)dr (8.18) 

Jo Jo 

or 2m(r) = F{r), since T’(O) = 0. This is because R = 0 when r = 0. The function 
F{r) thus represents the effective gravitational mass within radius r. For F =constant 
0 ), we hnd p = 0, so that the solution applies to empty space outside the boundary 
of the spherical distribution. In addition if / = 0 and tsf{r) = r are chosen then one 
may obtain Schwarzchild metric in a different co-ordinate system. The case F = 0 cor¬ 
responds to the absence of any gravitational field and by a suitable transformation the 
metric can be brought to Minkowskian form. 


An important property of the TBL solution is that the assignment of the arbitrary 
functions appearing in it over the interval from 0 to some ro completely determines the 
solution of the interior problem for the sphere of this radius [Landau et al, 1975]. It does 
not depend on how the functions are assigned for r > tq. Physical interpretation of the 
arbitrary functions becomes more clear when we write the equation (8.12) in the form 



iFjr) 

R 



which invites the interpretation of the quantities 


(8.19) 


E(r) = -/(r) 


T(t,r) = -R%r) 


V(t,r) = 


iF(r) 


m{r) 


( 8 . 20 ) 

( 8 . 21 ) 

( 8 . 22 ) 


R{t,r) 7?(t,r) 

as a conserved total energy at a co-ordinate radius r, the kinetic energy and potential 


energy respectively. 
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The TBL model has two kinds of singularities. The expression (8.13) shows that the 
density p diverges when either i? = 0 or i?' = 0. There are known as shell focusing 
and shell crossing singularities respectively. A shell focusing singularity {R = 0) on a 
shell of dust occurs when it collapses at or starts expanding from the centre of matter 
distribution. Now the moment t = tsf{r) corresponds to the arrival of the shell with a 
given radial co-ordinate r at the centre of matter distribution. Since the different shells 
of dust arrive at the centre at different times, there is always a possibility that any two 
shells of dust cross each other at a finite radius in course of their collapse. This is a 
situation when the comoving system breaks down and one encounters the shell crossing 
singularity (i?' = 0). On the hypersurface R = 0 [Seifert, 1979; Hellaby et al, 1984], 
both the matter density and Kretch ma un scalar diverges. But beyond the boundary 
surface it has negative density. So the region beyond it is unacceptable. It is therefore 
of interest to find condition which guarantee that no shell crossing will occur. A simple 
solution is to assure that the greater the radius of a shell, the latter it reaches the centre. 
In such a situation one must have > 0 to avoid shell crossing. 

In TBL model there are some important distinction between shell crossings where only 
one metric component goes to zero {prr) and the true big bang and recollapse surfaces 
where the angular components go to zero {geo and The shell crossing surfaces 

are always timelike whereas the bang or recollapse surfaces are always spacelike. The 
frequency shift of light coming from the shell crossing surfaces is hnite (red or blue) 
[Hellaby et al, 1984] and quite regular while that coming from the big bang surface is 
always infinite. The surface density a on shell crossing remains finite [Lake, 1984] while 
on big bang or recollapse surfaces a diverges. For these reasons it is believed that shell 
crossings are not serious physical singularities but they indicate the breakdown of the 
basic choice of comoving co-ordinate system in the Tolman-Bondi model. 


To investigate whether the naked singularities found by Eardley and Smarr [1979] 
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and Christodoulou [1984] go against the Cosmic Censorship Conjecture, Newman [1986] 
proposes an alternative formulation to the Cosmic Censorship Conjecture. He replaces 
the stability condition present in the original conjecture by the curvature strength criteria 
for the singularity. According to Tipler [1977] a curvature singularity is strong if any 
object hitting it is crushed to zero volume. Newman proposes that all naked singularities 
are in some sense gravitationally weak. The curvature strength of a singularity may 
be estimated in terms of the following limiting focusing conditions (LFC) [Clarke and 
Krolak, 1986]. Now it is convenient to study the strength of naked singularity, consider 
a null geodesic affinely parameterized by fi which terminates at the singularity where 
/i = 0. If the Ricci curvature along this geodesic grows as /i“^ or more rapidly near the 
singularity the later is said to be strong curvature in nature. Analytically, dehne 

T(/i) = RabK^K^ (8.23) 


where Rab is the Ricci tensor and is the four tangent to the geodesic. Then the 
following focusing conditions: 


the LFC 


Urn > 0 

/X —> 0 


(8.24) 


and the strong LFC 


Urn > 0 

/X —0 


(8.25) 


may be used to estimate the strength of the singularity. If only LFC is satisfied the 
singularity is essentially weak in nature. Fulhllment of strong LFC is equivalent to the 
termination of the geodesic in a strong curvature singularity in the sense of Tipler [1987]. 
Newman [1986] has simplified and generalized Christodoulou’s analysis to cover a wide 
class of Tolman-Bondi models including both the marginally bound and the time sym¬ 
metric cases. He has shown that both the shell focusing and shell crossing dust naked 
singularities that were found by Eardley and Smarr and by Christodoulou are not strong 
curvature singularities as defined by Tipler, Clarke and Ellis [1980]. 
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For dust, Newman has shown that naked shell focusing singularities are not strong 
curvature singularities, whereas Lake [1988] and independently Ori and Piran [1990] 
point out that the singularities found by Ori and Piran [1987] are strong curvature sin¬ 
gularities in the sense of Tipler. The addition of pressure in the spherical collapse of a 
perfect fluid should not lead to a stronger singularity than in the pressureless case. The 
problem is resolved by Waugh and Lake [1988] by considering the singularity structure 
of self-similar space-time, specihcally the marginally bound Tolman-Bondi case. By a 
clear and simple demonstration they showed that when the inhomogeneity parameter is 
greater than a critical value, the apparent horizon will always occur after the Cauchy 
horizon resulting a globally naked strong curvature shell focusing singularity. 

It may be interesting to point out that most of the space-times that have been studied 
to investigate the nature of singularity following a gravitational collapse are self-similar. 
It becomes important therefore to give a clear dehnition of self-similarity and also to 
discuss its advantages and disadvantages. Self-similarity is dehned by the existence of a 
homothetic killing vector held By a homothetic killing vector held we mean that in 
a local co-ordinate system 

+ ip,a = ^ ^5 'q/3 (8.26) 

So along each integral curve of that vector held all the points are similar to each 
other, apart from the scale factor that evolves linearly with the proper length along the 
integral curve in question. In specihc co-ordinate systems, self-similarity is manifested 
by a simple scaling relation for the metric functions. For instance, a spherical space- 
time is self-similar if there exists a radial area co-ordinate r and an orthogonal time 
co-ordinate t for which the metric functions satisfy 


gu{ct,cr) = gtt{t,r) 


(8.27) 


grr{ct,cr) = grr(t,r) 


( 8 . 28 ) 
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R{ct,cr) = cR(t,r) (8.29) 

for every c > 0. It is well known that Einstein’s eqnations admit self-similar solntions. 
In a spherical self-similar space-time it is easy to solve the eqnations of motion for radial 
nnll geodesics and this enables ns to stndy the cansal strnctnre of the space-time. The 
main difhcnlty arises when one tries to employ self-similarity to describe gravitational 
collapse in the asymptotic behavionr at inhnity. Self-similar solntions are not asymp¬ 
totically flat. In other words, in a self-similar solntion the mass of the collapsing object 
is inhnite. 


Now consider TBL model for bonnd, marginally bonnd and nnbonnd cases (/ < 0, / = 
0 and / > 0). For gravitational collapse sitnation of dnst, one reqnire R(t,r) < 0 and 
withont any loss of generality rescale R snch that 

/?(0,r) = 0 (8.30) 

Now partial derivatives of R like R' and R' are of crncial importance in the analysis. 
From eq.(8.14) solntion of R nsing the initial condition (8.30) have the form 


Vf 




So the shell focnsing singnlarity occnrs at the time 


(8.31) 


tsf{r) 


r3/2G'(-f) 


(8.32) 


The central singnlarity occnrs at r = 0, the corresponding time being to = F/(0). For 
the above three cases the expression for R', R' have the form 


1+ /3 -r] 



t 

r 


R' = iv-P)P 


R 


(8.33) 
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where 


k 


P A \ l + (3-r] 

2r 2rP2 [ yXTT 



t 

r 


(8.34) 


rF' r f 

R{t,r) = rP{t,r), rj = r]{r) = — {3 = P{r) = —, F{r) = rA(r) (8.35) 

R j 

At this point, the functions f{r), A(r) are analytic at r = 0 such that A(0) 7 ^ 0. It 
should be pointed out that when A(r) =constant and /(r) =constant, the space-time 
becomes self-similar in the sense of equations (8.27), (8.28) and (8.29) for any c > 0, 
whereas in general it is non-self-similar. The condition that A(0) 7 ^ 0 implies that 
tsf{0) = 0. It follows that the point r = 0, t = 0 corresponds to the central singularity 
on the t = 0 hypersurface where the energy density p blows up. 


Now introducing a new variable 

X = - (8.36) 

r 

the function P{t,r) = P{X,r) is given with the help of (8.31) and (8.32) by 

where tsf{r) = r0(r) = rG (—y) . 


The tangents ^ and 77* = ^ to the outgoing radial null geodesic with p as the 

affine parameter satisfy 


dK^ 

dp 


R' 

+ f 




0 


(8.38) 


dt _ K* _ P! 
dr~K^~ y/TTJ 


(8.39) 


The purpose is to hnd whether these geodesics terminate at the central singularity 
formed at r = 0, t = ts/(0) in the past. The exact nature of this singularity t = 0, r = 0 
could be analyzed by the limiting value of X = 7 at f = 0, r = 0. If the geodesics 
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meet the singularity with a dehnite value of the tangent then using equation (8.39) and 
L’Hospitals rule, one can have 

Xn = Urn X = Urn - = Um ^ = lim , 

r dr vl+/ 

t^O t^O t^o (8.40) 

r —>-0 r —>0 r — i >0 r —>0 

Now use the notation Aq = A(0), Pq = /9(0), /o = /(O) and Q = Q{X) = P{X,0), 

which is a function of X alone given by equation (8.37) 

with ©0 = ©(0) and Qo = Q{Xo). One can simplify equation (8.40) with the help of 
equations (8.31), (8.32), (8.33) and (8.41) as 


(1 ~ Po)Qo + 


Po 


.-y/Ao + /o 




^ + /o - Xo Vl + /o = 0 (8.42) 


If the equation (8.42) has real positive root then the singularity could be naked. In 
order to be the end point of null geodesics at least one real positive value of Xq should 
satisfy equation (8.42). Clearly if no real positive root of the above is found the singu¬ 
larity t = 0, r = 0 is not naked. The existence of a real positive root of equation (8.42) 
is a necessary and sufficient condition for the singularity to be at least locally naked. 
Such a singularity could be globally naked as well. The details of this will depend on 
the global features of A(r). 


Now for simplicity, consider /(r) = 0 and Aq 7 ^ 0. From equations (8.32) and (8.41) 
one may have 

e„ = Q(x) = of (jf - (8.43) 

Equations (8.42) and (8.43) imply that 


(3©o - Xo)=' = 27X3©2 (©o - Xo) 


(8.44) 












CHAPTER\^ General Introduction 


92 


This has real and positive roots if 


©0 > g + 5 a/3^ (8.45) 

and therefore the singnlarity is naked for which eqnation (8.45) is satished. The eqna- 
tion (8.44) has two real positive roots. Snch a sitnation is analyzed by Joshi and Dwivedi 
[1992] and it is shown that in this case the singnlarity will be globally naked. 


For the TBL model nnder consideration, using equations (8.12) and (8.15) the expres¬ 
sion of T in (8.23) becomes 


T = = 


F' 


r2p2ji> 


(8.46) 


Using L’Hospital’s rule and equations (8.33), (8.34), (8.38) and (8.39), one can see 
that 


lim = 

/i —> 0 


(l + /o)V2(ho-l)2Q^ 


> 0 


(8.47) 


where ho is a related with the constants /o, Aq, Pq, Qq and Xq. 


So equation (8.47) satisfies the strong limiting focusing condition (8.25). For this rea¬ 
sons, the naked singularity at t = 0, r = 0 is strong curvature singularity. 


Now to determine naked singularity in TBL model, one should discuss another ap¬ 
proach of geodesic study by Barve, Singh, Vaz and Witten [1999]. From equations (8.16) 
and (8.17) the solution of R for marginally bound case (/ = 0) is 


^ 3/2 


. 3/2 _ t 

2 


(8.48) 


It follows from the equation (8.13) that the function F{r) becomes hxed once the 
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initial density distribution p(0,r) = p{r) is given i.e., 


F{r) 


I 


p{r)r‘^dr 


If the initial density p(r) has a series expansion 


(8.49) 


p{r) = po + pir + +. (8.50) 

near the centre r = 0, the resulting series expansion for the mass function F{r) is 


F(r) = FqT^ + Fir‘^ + F 2 r^ + F-^r^ +. (8.51) 

where Fj = and j = 0,1, 2, 3,.... 

The central singularity at r = 0 forms at the time 

^0 = tsf{r) = (8.52) 

From equation (8.17) the singularity time at r, near r = 0 approximately written as 

tsfir) = to- r™ (8.53) 

3Fq 

Here F^ is the hrst non-vanishing term beyond Fq in equation (8.51). The hrst non¬ 
vanishing derivative in the series expansion in (8.50) should be negative and the density 
functions decreases as one moves out from the centre. So F^ < 0. From equation (8.53) 
it is to be seen that tsf{r) > to 


To examine the singularity at t = to, r = 0 is naked i.e., there are one or more 
outgoing null geodesics which terminate in the past at the central singularity, only 
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radial null geodesic is considered. Let the geodesic have the form 

t = tQ + a r°‘ (8.54) 

to leading order in (t, r) plane where a, a > 0. Now comparing (8.53) and (8.54) one 
may have a > m and if a = m then a < 

Now for TBL metric, an outgoing null geodesic must satisfy the equation 

^ = R! (8.55) 

dr 

Consider hrst a > m\ In this case, equations (8.48), (8.52), (8.54) and (8.55) near 
r = 0 give the results 


aar“ ^ = ( 1 + 


2m\ 

\ 2Fo 
2/3 


TP \ 2/3 , 

TTI \ 2771 


'f 3 


(8.56) 


which implies a = 1 + ^, a = (—^] • So for a > m, m takes values 1, 2 
only. Hence the singularity is naked for m = 1 and 2 i.e., for the models pi < 0 and 
pi = 0, p 2 < 0. There is at least one outgoing geodesic given by (8.54) which terminates 
in the central singularity in the past. If m > 3 then the condition a > m cannot be 
satished and the singularity is covered. 


Consider next that a = m. In this case also equations (8.48), (8.52), (8.54) and (8.55) 
give the result 


mar"" ai/F^ - 


(2m + 3) F„ 




6 -/fi/ V SFj 

which implies m = 3 and a satishes the biquadratic equation 


■\ ^ 2m 

'f 3 


(8.57) 



- ‘TIP - '^ 3 -^3 


0 


0 


0 


(8.58) 
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By substituting b = ^, ^ this equation can be written as 

463(36 + 0 - + = 0 (8.59) 

The singularity will be naked if this equation admits one or more positive roots for 6 
which satishes the constraint 6 < —| and also ^ satishes the condition ^ < —25.9904. 
Now to investigate whether there is an entire family of radial null geodesics which ter¬ 
minate at the naked singularity, a solution for the geodesics correct to one order beyond 
the solution (8.54), is assumed 


t = to + a r'^ + d r°‘^^ (8.60) 

where d, /3 > 0. In this case, a, a are same as in before and d is totally arbitrary. 
There will be an entire family of outgoing null geodesics terminating at the singular¬ 
ity provided fd is non-negative. It is essential that (3 be non-negative, otherwise these 
geodesics will not lie in the space-time. 

The normal vector to the boundary of any hypersphere R — Rq = 0 is given by 
Ifj, = {R, R', 0, 0). In order that is a null vector and the boundary of this hypersphere 
is a null surface, so 

Ritah{r),r) = F{r) (8.61) 

where t = tah{r) is the equation of apparent horizon which marks the boundary of 
the trapped region. If the apparent horizon starts developing earlier than the epoch of 
singularity formation then the event horizon can fully cover the strong gravity regions in¬ 
cluding the hnal singularity which will thus be hidden within a black hole. On the other 
hand, if trapped surfaces form sufficiently later during the evolution of collapse then it is 
possible for the singularity to communicate with outside observers. That means, naked 
singularity or black hole appears if tah{r) > R or tah(r) < to respectively near r = 0. 
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From equations (8.48) to (8.53) and (8.61), apparent horizon time is given by 

tah{.r) = to - \FQr'^ -+. (8.62) 

For m < 3, equation (8.62) shows that tahir) > to and so naked singularity forms. For 
m = 3, naked singularity or black hole forms if F 3 < —2Fq^‘^ or F^ > —2Fg^^ respec¬ 
tively. For m > 4, equation (8.62) shows that tah{f) < to and so a black hole forms. 


The expression for shear scalar is given by 


^2m 

As expected, the shear diverges as one approach to the singularity. 


= 


I R' R 


3 \ R R 


m^F^ 

12Fo 


1 - 3Fl/h + ?Fot2 


(8.63) 


8.1.2 Perfect Fluid Collapse 

The perfect fluid matter model is one of the most natural ways to introduce matter 
pressure. If the pressure is bounded from above, there can appear a shell crossing naked 
singularity [Muller Zum Hagen et al, 1974] which is gravitationally weak. There can 
appear a strong curvature naked singularity even with unbounded pressure. 

The examples presented before are all based on collapse of pressureless matter based 
on TBL model and one may reasonably think that these naked singularities will not be 
covered by an event horizon when pressure is introduced. Also from physical point of 
view the assumption that matter is dust like is never admissible near the singularity, 
since p ^ 00 and one should use the ultra-relativistic equation of state p = ^ p. However, 
it appears that the general character of the limiting laws of compression are to a large 
extent independent of the equation of state of the matter [Lifshitz et al, 1960, 1961]. 
Ori and Piran [1987] have investigated the self-similar collapse of a perfect fluid with an 
equation of state p = (7 — l)p. They have shown that if the equation of state is soft 
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enough i.e., 7 — 1 -C 1, a naked singularity may appear. In another work, Ori and Piran 
[1990] used the self-similarity to model spherically symmetric gravitational collapse of 
compact objects with a reasonable equation of state. It is shown that a p = (7 — l)p 
equation of state is compatible with self-similarity. They show that in a signihcant part 
of the space of self-similar solutions, there is a globally naked central singularity from 
which null geodesics emerge to infinity. 

For 1 < 7 < 1.4, they found numerically that there is a discrete set of self-similar 
solutions that allow analytic initial data beyond a ‘sonic point’. These solutions can be 
labeled by the number of zeros in the velocity held of world lines of constant circumfer¬ 
ential radius relative to the huid element. There exists a pure collapse solution among 
these self-similar solutions which they call the general relativistic Larson-Penston solu¬ 
tion. They showed that a central naked singularity forms in this Larson-Penston solution 
for 1 < 7 < 1.0105. They also showed that there are naked singular solutions with oscil¬ 
lations in the velocity held for 1 < 7 < 1.4. The results of this work were conhrmed and 
extended to 1 < 7 < 25/16 [Foglizzo et al, 1993]. This naked singularity is ingoing null. 
It was shown that this naked singularity satishes both the limiting focusing condition 
and the strong curvature condition for the hrst null ray [Waugh et al, 1988, 1989]. 

In order to judge the necessity of self-similarity assumption, Harada [1998] numerically 
simulated the spherically symmetric and adiabatic gravitational collapse of perfect huid 
with the same equation of state p = (7 — l)p without this assumption. Since null co¬ 
ordinates were used in these simulations, he could detect naked singularities not relying 
upon the absence of an apparent horizon. The result was that naked singularities develop 
from rather generic initial data sets for 1 < 7 ^ 1 . 01 , which is consistent with the 
result of Ori and Piran obtained using the self-similarity assumption. In fact, through 
numerical simulations by Harada and Maeda [2001], it was found that generic spherical 
collapse approaches the Larson-Penston self-similar solution in the region around the 
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centre, at least for 1 < 7 < 1.03. This finding is supported by a linear stability analysis 
of self-similar solutions. Although the hnal fate of the generic spherical collapse of a 
perfect fluid with larger values of 7 is not known, Harada [2001] analytically showed 
that the Larson-Penston self-similar solution is no longer stable. Recently Ghosh et al 
[2002, 2003], Joshi et al [2002] and Goswami et al [2002, 2003] gave the solutions for 
perfect fluid and they have shown that perfect fluid halts collapse. 


8.1.3 Other Examples on Spherical Collapse 

Gritical behaviour in gravitational collapse was discovered by Ghoptuik [1993]. He found 
so-called critical behaviour, such as a scaling law for the formed black hole mass. He also 
found that there is a discrete self-similar solution that sits at the threshold of black hole 
formation which is called a critical solution. Similar phenomena have been observed in 
the collapse of various kinds of matter helds, for example, axisymmetric gravitational 
waves [Abrahams et al, 1993], radiation fluid [Evans et al, 1994] and more general per¬ 
fect fluids [Neilsen et al, 2000, 2001]. One interesting example of spherical collapse is 
a spherical self-gravitating system of counter rotating particles i.e., an Einstein cluster. 
One should mention several additional examples of naked singularities in type I and type 
H matter. If one consider imploding “null dust” into the centre, the space-time is given 
by the Vaidya metric [Vaidya, 1943, 1951]. For this space-time, it has been shown that 
naked singularity formation is possible from regular initial data [Ghosh et al, 2000, 2001; 
Beesham et al, 2003]. There have been many analysis with a set of assumptions and the 
conditions for the appearance of a naked singularity are written down in terms of the 
energy density, radial stress and tangential stress [Goncalves et al, 2001, 2002]. In these 
analysis, in general the conclusion is that naked singularities are possible for generic 
matter helds that satisfy some energy conditions. In this approach, this conclusion is 
very natural because there remains great freedom in the choice of matter helds, even if 
some energy condition is imposed. 
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8.2 Non-Spherical Collapse 

There is very little progress in studying non-spherical collapse. The quasi-spherical dust 
collapse models given by Szekeres metric [Szekeres, 1975] were analyzed by Szekeres him¬ 
self [1975]. It has been shown that shell focusing singularities are possible in the solution 
and the conditions necessary for the appearance of naked singularities are very similar 
to those in the case of spherical dust collapse [Joshi and Krolak, 1996; Goncalves, 2001]. 
The global visibility of this singularity was also examined [Deshingkar et al, 1998]. Now 
the basic difficulty is the ambiguity of horizon formation in non-spherical geometries and 
the influence of gravitational radiation. Schoen and Yau [1983] proposed a sufficient cri¬ 
terion for the formation of trapped surfaces in an arbitrary space-time but it fails to say 
anything about the conditions which lead to the formation of naked singularities. This 
problem has been restated by Thorne [1972] in the form of a conjecture known as hoop 
conjecture which states that “horizon form when and only when a gravitational mass M 
gets compacted into a region whose circumference in every direction is C < AnM”. He 
analyzed the causal structure of cylindrically symmetric space-times and found remark¬ 
ably different nature from spherically symmetric space-times. The event horizon cannot 
exist in cylindrically symmetric space-time. Therefore, if an inhnitely long cylindrical 
object gravitationally collapses to a singularity, it becomes a naked singularity, not a 
black hole covered by an event horizon. Then it is natural to ask what happens if a very 
long but hnite object collapses. The hoop conjecture was derived from such a thought 
experiment. It should be noted that the hoop conjecture itself does not assert that a 
naked singularity will not appear. It is a conjecture regarding the necessary and suffi¬ 
cient condition for black hole formation. 

Many examples of naked singularity formation are obtained from the analysis of spher¬ 
ically symmetric space-time. These examples are consistent with the hoop conjecture. 
The circumference C corresponding to a radius R centered at the symmetric centre is 
2tiR. For a central naked singularity of a spherically symmetric space-time, the ratio of 
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C to the gravitational mass M within R can be estimated as 


lim 


2itR 

AttM 




r 


0 


r 


0 


r 


0 


where is a constant and 1 < a < 3. In this case the condition C > AnM coincides 
with the condition that the centre is not trapped. Therefore the appearance of a central 
naked singularity in spherical gravitational collapse does not constitute a counter exam¬ 
ple to the hoop conjecture. 

As for the non-spherically symmetric case several studies support the hoop conjecture. 
A pioneering investigation is the numerical simulation of general relativistic collapse of 
axially symmetric stars [Nakamura et al, 1982]. From initial conditions that correspond 
to a very elongated prolate fluid with sufficiently low thermal energy, the fluid collapses 
with this elongated form maintained. In the numerical simulation, black hole was not 
found. This result suggests that such an elongated object might collapse to a naked 
singularity and in this case it would provide evidence supporting the hoop conjecture. 

Shapiro and Teukolsky [1991, 1992] numerically studied the evolution of a collisionless 
gas spheroid with fully general relativistic simulations. They found some evidence that 
a prolate spheroid with a sufficiently elongated initial conhguration and even with a 
small angular momentum, might form naked singularity. They also found that when a 
spheroid is highly prolate, a spindle singularity forms at the pole and the spindle might 
be a naked singularity. Subsequently, there were attempts namely, gravitational radia¬ 
tion emission in aspherical collapse [Nakamura et al, 1993], analytical studies of prolate 
collapsing spheroids [Barrabes et al, 1991; Pelath et al, 1998] and others [Harada et al, 
1998; Iguchi et al, 1999, 2000] to prove or disprove the conjecture. Interestingly, all of 
them either conhrmed or failed to refute the conjecture. 
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It should be noted that there are problems in the hoop conjecture to prove it as a 
mathematically unambiguous theorem, such as how one should dehne the mass of the 
object and the length of the hoop. Although nobody have found a counter example to 
the hoop conjecture, it is necessary to obtain a complete solution to these problems. 


Chapter 9 


Naked Singularities in Higher 
Dimensional Spherical Gravitational 
Collapse 

9.1 Prelude 

The Cosmic Censorship Conjecture [Penrose, 1969, 1979] is yet one of the unresolved 
problems in General Relativity. It states that the space-time singularity produced by 
gravitational collapse must be covered by the horizon. However the singularity theo¬ 
rems as such do not state anything about the visibility of the singularity to an outside 
observer. In fact the conjecture has not yet any precise mathematical proof. Several 
models related to the gravitational collapse of matter has so far been constructed where 
one encounters a naked singularity [Waugh et al, 1988; Lemos, 1991, 1992; Joshi et al, 
1999; Ilha et al, 1997, 1999; Rocha et al, 2000]. 

It has recently been pointed out by Joshi et al [2002] that the physical feature which 
is responsible for the formation of naked singularity is nothing but the presence of shear. 
It is the shear developing in the gravitational collapse, which delays the formation of 
the apparent horizon so that the communication is possible from the very strong gravity 
region to observers situated outside. Joshi et al [2002] have analyzed in details how 
the presence of shear determines the growth and evolution of an inhomogeneous dust 
distribution represented by four dimensional Tolman-Bondi metric and have attempted 
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to clarify the nature of singularities as the hnal outcome. 

The objective of this chapter is to fully investigate the situation in the background of 
higher dimensional space-time with both non-marginally and marginally bound collapse 
[Banerjee et al, 1994; Ghosh et al, 2000, 2001; Deshingkar et al, 1999]. As it is only 
to be expected, in one way or another, these works all deal with propagation of null 
geodesics in the space-time of a collapsing dust [Barve et al, 1999]. In this context we 
mention that Ghosh and Beesham [2000, 2001] have also studied dust collapse for (n-|-2) 
dimensional Tolman-Bondi space-time for marginally bound case (/ = 0), considering 
the self-similar solutions. They have concluded that higher dimensions are favourable 
for black holes rather than naked singularities. Also recently, Ghosh and Banerjee [ 2002 ] 
have considered non-marginal case (/ ^ 0) for dust collapse in 5D Tolman-Bondi model 
and have shown that the degree of inhomogeneity of the collapsing matter is necessary 
to form a naked singularity. 

In this chapter, we show that in non-marginally bound collapse (i.e., for / 7 ^ 0), the 
naked singularity may appear in any dimensional space-time, but for marginally bound 
collapse the naked singularity may appear only when the space-time has dimensions 
upto hve. In more than hve dimensions and marginally bound case tah < to that is the 
apparent horizon at any point in the region forms earlier than the central shell focusing 
singularity, which indicates that the shell focusing singularity hrst appearing at r = 0 
remains hidden behind the apparent horizon and so gives rise to a black hole always. 
Also we study gravitational collapse for non-marginally bound case considering non self¬ 
similar solutions and it is generalization to higher dimension of the Dwivedi and Joshi 
[1992]. The geodesic equations can not completely solved due to the presence of the 
complicated hypergeometric function. So we present numerical results which favour the 
formation of naked singularity in any dimension. 
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9.2 Basic Equations in Higher Dimensional TBL Model 

The {n + 2) dimensional Tolman-Bondi type metric is given by 


= e'' — e^dr"^ — R^dVtl^ (9.1) 

where v, A, R are functions of the radial co-ordinate r and time t and dVt^ represents 
the metric on the ?7,-sphere. Since we assume the matter in the form of dust, the motion 
of particles will be geodesic allowing us to write 

= 1 (9.2) 

Using comoving co-ordinates one can in view of the field equations [Banerjee et al, 1994] 
arrive at the following relations in {n + 2) dimensional space-time 


e^ = 


R 


n 


and 


R^ = f{r) + 


1 + /(r) 
F{r) 


R 


1-1 ’ 


(9.3) 


(9.4) 


where the function /(r) classihes the space-time [Joshi, 1993] as bound, marginally 
bound and unbound depending on the range of its values which are respectively 


f{r) < 0, /(r) = 0, /(r) > 0. 


The function F{r) can be interpreted as the mass function which is related with the 
mass contained within the comoving radius r. 

The energy density p{t, r) is therefore given by 


p(f,r) 


nF'{r) 

2R^R' 


(9.5) 


The above models are characterized by the initial data specihed on the initial hyper¬ 
surface t = ti, from which the collapse develops. As it is possible to make an arbitrary 
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relabeling of spherical dust shells by r — g{r), without loss of generality we £x the 
labeling by the choice R(ti,r) = r, so that initial density distribution is given by 

, , nF' , , 

«(’-) = ^ ( 9 ' 6 ) 

Now the curve t = tgf{r) dehnes the shell-focusing singularity and is characterized by 

Ritsf{r),r) = 0 (9.7) 

It indicates that the shell focusing singularity occurs at different r at different epochs. 
So the collapse in this case is not simultaneous in comoving co-ordinates and is in 
variance with that in the homogeneous dust model [Ilha et al, 1997, 1999; Chatterjee et 
al, 1990], 

Further, within the collapsing cloud, the trapped surfaces will be formed due to un¬ 
bounded growth of the density and these trapped surfaces are characterized by the 
outgoing null geodesics. We observe that the normal vector to the boundary of any 

hypersphere R — Rq = 0 is given by = (77, 7?', 0, 0, 0,.). In order that is a null 

vector and the boundary of this hypersphere is a null surface we must have 

R{tah{.r),r) = [F(r)]^ (9.8) 

where t = tah{f') is the equation of apparent horizon which marks the boundary of the 
trapped region. 

To characterize the nature of the singularity we shall discuss the two possibilities 
namely, (7) tah < 4/(0) and (77) tah > 4/(0). The hrst case may correspond to forma¬ 
tion of black hole while the second one may lead to naked singularity. If the apparent 
horizon will form earlier than the instant of the formation singularity, then the event 
horizon can fully cover the strong gravity region and also the singularity. As a result, no 
light signal from the singularity can reach to any outside observer and the singularity is 
totally hidden within a black hole. On the other hand, in the second case the trapped 
surfaces will form much later during the evolution of the collapse and it is possible to 
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have a communication between the singularity and external observers. 


9.3 Marginally Bound Case 

In this case (i.e., /(r) = 0) equation (9.4) can be integrated out easily to give 


_ n+l n-i-1 Tl \ / ,—r , , 

H 2 = r 2 -y^F(r) (t — tj) 


(9.9) 


where we have used the initial condition R{ti, r) = r. 

Combining (9.8) and (9.9) along with the initial condition R{ti, r) = r we arrive at 
the relation 

9 9 

(9.10) 


W N 2 n + l ^-1/2 2 , , 1 

r) — ti = - r 2 F - 

^ n+l n+l 


Since our study is restricted to the region near r = 0, the mass function F[r) should 
vanishes exactly at r = 0 and it is possible to express F{r) as a polynomial function 
of r near the origin. One should at the same time keep in mind that initially the 
central density is regular. Further from physical considerations one may argue that 
p'(r) should vanish exactly at r = 0, but is negative in the neighbouring region. All 
these consideration lead to the following expressions for F{r) and the initial energy 
density pj(r), 

F(r) = + 0 ^”+^ + . (9.11) 


and 

Pi{r) = po + pir + p2r‘^ + . (9.12) 

with pj = j = 0,1,2,... . Obviously the initial central density is 

given by Pc = Since p'(r) < 0 in the region r 0 one must have +2 < 0. 

Now in view of (9.7) and (9.9) one can write 




2 ^++i )/2 

(n + 1) |f„r”+> + fjr”+3 + 


(9.13) 
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If we now denote to = ts/(0) as the instant of first shell focusing singularity to occur 
at the centre r = 0, the relation (9.13) can be used to obtain the following relation 


to ti 


(n 


m 


1/2-' 


which in combination with the relation (9.10) yields 


(9.14) 


tah{r) - to 


J: _^ 2 , 

n + l ^ 


1 Fo 
n + l 


+ 


1 Fj 4 
n + l +3/2 ^ 


2 

n+l 



n+l 

J-n-l -|- 


2 +2 3 ti -1 

- T n—1 

9 1 n —2 ' 

n"' — 1 p— 

-^0 


+ 



1 


Fo 


An-2 
f n — 1 



(9.15) 


The relation (9.15) is valid in general (n + 2) dimensional space-time and reduces to the 
4 dimensional expression (n = 2) in the region near the centre as given in the paper of 
Joshi et al [2002]. In 5D space-time the naked singularity appears only if +2 7 ^ 0 and 
further IF 2 I > 2 + 9 . Otherwise for either +2 = 0 or IF 2 I < 2+q there must occur a black 
hole, because in that case tahi'f') < to that is, the horizon appears earlier than the shell 
focusing singularity at r = 0. Only in a very special case of IT 2 I = 2Fq the occurrence 
of naked singularity or black hole will depend on the co-efficients of higher powers of r. 
For general (n + 2) dimensions the close examination of (9.15) reveals that if there is to 
exist a naked singularity we must have 


n+l 
n — 1 


> 2 , 


(9.16) 


which means n < 3. It is interesting to note that has the value 2 when n = 3 

and then it decreases monotonically with the increasing number of dimensions. So for 


n + l 


space-time with dimension larger than hve the term dominates leading to 


a black hole. Our conclusion is that for the space-time having larger than hve dimen¬ 
sions the existence of the naked singularity is prohibited or in other words the shell 
focusing singularity is fully covered by the horizon in such cases. However, if we relax 
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the restriction namely, p'(r) = 0 at r = 0 then Fi will be non-zero and this will leads to 
the possibility of naked singularity in any dimension. 

9.4 Non-Marginally Bound Case 

In this case (i.e., /(r) ^ 0), integrating equation (9.4) with the help of initial condition 
R{ti, r) = r, we have the form of R as 





F 


(9.17) 


where 2-^1 is the usual hypergeo metric function with a = \ + Using equations 
(9.7) and (9.8) separately in equation (9.17) we have 



(9.18) 


and 



1 

— , a, a 1 , 



The only new input in the present case is the function /(r), which may be expressed as 
a power series in r near the centre r = 0. 


We assume 


f{r) = for^ + fir^ + f2r^ + 


(9.20) 


The form (9.20) is chosen because /(r) vanishes as r —0, which is demanded by the 
regularity condition at r = 0 [Harada et al, 2000, 2002]. 


From equation (9.18) we have the central shell focusing singularity time 



(9.21) 
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Therefore, equation (9.19) and (9.21) give (after using (9.11) and (9.20) in (9.19)) 


tah{r) -to = 


3/2 ® ® + 2 , —+ 

(3n - 4)^0^^ 2 Fo 


{n + 


n, n T 1, —T-- 

2 Fo^ ^n-l)i3n-A)foF^^\fo + Fo) 


(4 - 3n)flF^ 2^i[ 7:, a, a + 1, - 


/o 


+{-2/oFo((3 - 2n)fl + 2(n - l)/ 2 Fo) - 4(n - l)flF,{h - F^) + 4(n - l)flF 2 


+ (3n — 4)/^^Fo^} 2 Fi[-, a + 1, a + 2, —^] 


+ 0{r^) (9.22) 


From the above relation we note that near r = 0, tahi^) > or < 0 according as 

/i 2-^1 [-, a + 1 , a + 2 , —^] < or > 0 . 

/ no 

Though the hypergeometric function has arguments depending on n yet 2 F\ is always 
positive. So the above restriction simply implies fi < or > 0. Now according to Joshi 
etal, if a comoving observer (at fixed r) does not encounter any trapped surfaces until 
the time of singularity formation then it is possible to visualize the singularity otherwise 
the singularity is covered by trapped surfaces, leading to a black hole. Thus /i < 0 
leads to naked singularity while /i > 0 is the restriction for black hole. However, the 
above conditions are not sufficient, particularly for non-smooth initial density prohle. 
Therefore, we may have locally naked singularity or black hole depending on the sign of 
/i- 


Calculation of shear near r = 0: 


We still consider the non-marginally bound case /(r) 7 ^ 0. The shear scalar in the 
{n + 2 ) dimensional spherically symmetric dust metric so far discussed in this section 
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may be estimated by the factor a = 2 (n+i) ~ Rj- Using (9.3) and (9.4) one gets in 

turn 


a = 


n 


2 (n + l) 




n 


[{RF' 


8 (n + 1) 


(n + l)R'F} + R^-\Rf - 2R'f)] 

(F + 

(9.23) 


Since at the initial hypersurface {t = U) we have chosen R{ti,r) = r, so the initial 
shear cXj is of the form 


a = 


n — 2 [{rF' 

8 (n- 1) 


(n — 1)F} + r"" ^{rf — 2/)] 

(F + 


(9.24) 


Thus using the power series expansions (9.11) and (9.20) for F{r) and /(r) in the 


above expression for cxj we have 


n — 2 


(Ti = 


fn + YZ=2 ^ifm + Um)U 


a/(/o + Fq) + fir + Ylim=2 '^ifm + Fm)F 


n — 2 


8{n — 1) fo + Fq [_ 


hr + 2(/2 + F 2 ) - 


fl 


2{fo + Fo) 


+ O(r^) 


(9.25) 


For marginally bound case /(r) = 0 i.e., /o = /i = /2 = ••• = 0. So the shear in this 
case becomes 

Z^m =2 


a = 


n 


8 (n + 1) ^V2 


1 + R^Foit - uy - (n + l)F^^\t - U) 


(9.26) 


and the initial shear aX t = ti is given by 


(Ji = 


n 


'TZ=2^Fm^ 


8(n + 1) 


, 1/2 


We note that the initial shear (9.25) vanishes when /i = 0, (F 2 + / 2 ) 


(9.27) 

(U 3 + h) = 


. = 0 and hence even if the initial shear is zero the dust distribution may be inho¬ 
mogeneous. It is interesting to note that the dependence on r of the initial shear does 
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not depend on the number of dimensions of the space-time. In fact the expression (9.27) 
coincides exactly with the value of the initial shear at calculated earlier by Joshi et al 
[2002] for AD space-time (n = 2). In view of what has been discussed so far the expres¬ 
sion (9.27) reveals that the existence of the naked singularity is directly related with the 
non-vanishing shear in four and hve dimensions. But in larger dimensions the central 
singularity seems to be covered by the appearance of apparent horizons. Also the above 
expression (9.25) for cxj reveals that the existence of the naked singularity is not directly 
related to the non-vanishing of the shear as it does in the marginally bound case. In fact 
the shear in view of (9.26) increases and goes to inhnity at this epoch, which is expected. 


9.5 Geodesic Study in Marginally Bound Case 


Since the radius of the spherical shell R shrinks to zero at the time tsf{r) so from (9.13) 
we have (chosen = 0) 


tsf{r) 


o "+1 

2 r 2 

n + 1 ^F{r) 


(9.28) 


Now the Kretchmann scalar 


F'2 p pi p2 

K = [n(n — 1) -|- 1]— -- — 2n^{n — 1)—-—-—- -|- n^in^ — 1)—-—- (9.29) 

I- ' > I J^2nJ^i2 ' ' J^2n+lJ^i2 ^ ' J^2n+2 ^ > 

diverges at f = tsf{r) i.e., R = 0. Thus it represents the formation of a curvature 
singularity at r. In fact the central singularity (i.e., r = 0) forms at the time 

=\ I Tu (**■“) 

V (^ + l)Po 

The Kretchmann scalar also diverges at this central singularity. 

Now if we use the expansion (9.11) for F{r) in equation (9.28) then near r = 0, the 
singularity curve can be approximately written as (using (9.30)) 


Frr 


(n 


m. 


3/2 


«•/(>■) = to 


(9.31) 
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where m >2 and Fm is the first non-vanishing term beyond Fq. Thns tsf{r) > to as 
Fm < 0 for any m > 2. 

To examine whether the singnlarity at t = 0 is naked or not, we investigate 

whether there exist one or more ontgoing nnll geodesics which terminate in the past at 
the central singnlarity. In particular, we will concentrate to radial null geodesics only. 


Let us start with the assumption that it is possible to have one or more such geodesics 
and we choose the form of the geodesics (near r = 0 ) as 

t = to + ar^", (9.32) 

to leading order in t-r plane with a > 0, a > 0. Now for t in the geodesic (9.32) should 
be less than tsf{r) in (9.31) for visibility of the naked singularity so on comparison we 
have 

F 

a > m and a < - — — 5777 . (9.33) 

(n + 1 )fP 

Also from the metric form (9.1), an outgoing null geodesic must satisfy 


dt 

dr 


= R' 


(9.34) 


But near r = 0, the solution (9.9) for R simplihes to 


R = r 




2 

n+1 


Thus combining (9.32) and (9.35) in equation (9.34) we get 


(9.35) 


aar"-^ = 


1-^,/Fo {to + ar") - + ar“) 


1 - (1 + ) (to + ar“) 


n —1 
n+1 


(9.36) 


Now if there exists a self consistent solution of this equation then it is possible to have 
at least one outgoing radial null geodesic that had started at the singularity i.e., the 
singularity is naked. In order to simplify the above equation we shall use the restrictions 
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in equation (9.33) in the following two ways: 


(i) a > m : 

The equation (9.36) becomes (in leading order) 


aar“-^ 



2m \ 
n+lj 


Fm 

2Fo 


2 

n+1 


2m 

n+1 


(9.37) 


which implies 

P \ 2/(n+l) 

+) (9.38) 

Thus for a > m we have m < an^l n > 2. 

n—1 — 

For n = 2, m may take values 1 and 2 and we have pi = Fi = 0 and p 2 < 0, +2 < 0. 
As a result, a is real and positive from equations (9.33) and (9.38). Moreover, these re¬ 
strictions are already assumed in the power series expansion for Pi{r) so that the initial 
density gradient is negative and falls off rapidly near the centre. But for n > 2, m = 1 
is the only possible solution for which no real positive solution of ‘a’ is permissible from 
equations (9.33) and (9.38). Hence with the restriction a > m, we have a consistent 
solution of equation (9.36) only for four dimensional space-time i.e., it is possible to have 
(at least) null geodesics terminate in the past at the singularity only for four dimension 
and we can have naked singularity for n = 2. 


ot — 1 -|- 


2m 

n+1 


and a = 


{ii) a = m : 

In this case equation (9.36) simplihes to 


mar 


m—1 


Fm 

2+0 


_ 1 — n 



{2m + n + l) Fm 
2{n + 1) +0 


n -|- 1 
2 



a 


2m 

n+1 


(9.39) 


A comparative study of equal powers of r shows that m = and a depends on Fm 


and + 0 - Here for n = 2, m = 3 and this situation is already discussed by Barve et al 
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[1999]. For n = 3, we have m = 2 and from (9.39) we get 

262(46 + 0 + (26 + 0 ^ = 0 ( 9 . 40 ) 

where 6 = and ^ 

We note that for real 6 , we must have 6 < — | i.e., a < - 372 , which is essentially 

the restriction in (9.33). It can be shown that the above cubic equation has at least 
one positive real root if e < -(11 + 5^5). Thus, if +2 < -(11 + 505)^2 we have at 
least one real positive solution for a which is consistent with equation (9.36) (or (9.39)). 
Further for n > 3, we can not have any integral (positive) solution for m and hence 
equation (9.39) is not consistent for n > 3. So, it is possible for (at least) radial null 
geodesics which initiate from the singularity and reach to an external observer without 
get prevented by any trapped surface for n <3. Therefore, naked singularity is possible 
only for four and hve dimensions and for higher dimensions [n > 4 ) all singularities are 
covered by trapped surfaces leading to black hole. 


Further, to examine whether it is possible to have an entire family of geodesics those 
have started at the singularity, let us consider geodesics correct to one order beyond 
equation (9.32) i.e., of the form 

t = to + ar" + (9.41) 

where as before a, d, a and [3 are positive constants. Thus equation (9.36) is modihed 
to 

aar"-^ + {a + = 


n+l 

2 


■a/F]) (to + + dr"+^) — + dr"+^) 

n—1 


1 - ^^/Fo (^1 + (to + ar^ + dr^+f^) 

So for a > m (retaining terms upto second order) we have 


n—1 

n+l 


(9.42) 


1 + 


2m 


2+0 


2 

n+l 


rn+i + _Dr" 


(n —l)m 
n+l 


aar“-^ + (a + ;5)dr“+^-^ = 


n + l 


(9.43) 
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As before, we have the values of a and a in equation (9.38) and 



n + 1 


r\ I (3—n)m 
^ n+1 


D 


with 





— n — 1 


(9.44) 


As we have similar conclusion as before for a > m so we now consider the case a = m. 
But if we restrict ourselves to the hve dimensional case then m = a = 2 and we have 
the same cubic equation (9.40) for b. Now fd can be evaluated from the equation 



(9.45) 


and we must have /d > 0, otherwise the geodesics will not lie in the real space-time. 
As there is no restriction on d so it is totally arbitrary. This implies that there exists 
an entire family of outgoing null geodesics terminated in the past at the singularity for 
four and hve dimensions only. 

9.6 Geodesic Study in Non-Marginally Bound Case 

In the present discussion we are concerned with the gravitational collapse, we require 
that R{t,r) < 0. As it is possible to make an arbitrary relabeling of spherical dust 
shells by r ^ g{r) without loss of generality, we hx the labeling by requiring that on the 
hypersurface t = 0, r coincides with the radius 


R{0, r) = r 


(9.46) 


So the initial density is given by 



(9.47) 


Now from eq.(9.5) it can be seen that the density diverges faster in 5D {n = 3) as 
compared to 4D (n = 2). For increasing dimensions the density diverges rapidly. Hence 
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there is relatively more mass-energy collapsing in the higher dimensional space-time 
compared to the 4D and 5D cases. 

In this case (i.e., /(r) 7 ^ 0), equation (9.17) can be rewritten as 


t = 


(n 


1 )v/F . 


n+l „ rl , fr 

r 2 2Fi[-,a,a + l, -: 


n —1 


„n+i 1 fK 

R 2 2 -^ i [2 5 ® + 1 )- 'p 


n —1 


(9.48) 


Let us assume [Dwivedi, 1992] 


F{r) = r” ^A(r) 


a = a(r) = 

> 

13 = f3{r) = ^ 

R{t,r) = rP{t,r) 

So using equations (9.4), (9.48) and (9.49), we have the following expressions 


(9.49) 


R' = 


n — 1 


If 1 fpn-l 1 r 

P{f3 — a) + — lP a, a + 1,- —]~P ^ a, a + 1, —y] 


fpn-i [a - {3 P n - 1)P 2 + 

{a{n + l)-2(3}\ll + L^p - V 


1 + 


(9.50) 


and 


k = - 

r 


If 1 fpn-l 1 r 

I 2 -P'i[-, a, a-I-1,- - —]—P 2 2 ^ 1 )-, a, a-f 1, ——] {a(n-|-1) — 2/3} 


n+l 


, /P-i , (a-i3 + n-l)P-^ 
-a\ l-\ -;-^- 


A 


i + f 


(9.51) 


When A(r) = constant and /(r) = constant, the space-time becomes self-similar. Now 
we restrict ourselves to functions /(r) and A(r) which are analytic at r = 0, such that 
A(0) > 0. From eq.(9.5) it can be seen that the density at the centre (r = 0) is hnite at 
any time t, but becomes singular at t = 0 . 
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We wish to investigate whether the singularity, when the central shell with co-moving 
co-ordinate (r = 0) collapses to the centre at the time t = 0, is naked .The singularity is 
naked if there exists a null geodesic that emanates from the singularity. Let ^ 

be the tangent vector to the radial null geodesic, where /i is the affine parameter. Then 
we derive the following equations: 


dK^ 

d^i 


R' 

+ f 




0 


(9.52) 


dt _K^ _ R! 

Let us dehne X = ^, so that P{t,r) = P{X,r). So eq.(9.48) becomes 


(9.53) 


2 r 1 f n+l 1 fpn-1 'j 

" bWyt {"■^‘' 2 ’ “■ “ + 1> -^1 - 2 n[ 2 . a. a + 1 . ^11 (9-54) 

The nature of the singularity can be characterized by the existence of radial null 
geodesics emerging from the singularity. The singularity is at least locally naked if there 
exist such geodesics and if no such geodesics exist it is a black hole. If the singularity is 
naked, then there exists a real and positive value of Xq as a solution to the equation 


^0 


lim - = lim ^ = Urn 

r dr vl+ 

t—>-0r—f—>0r—>0 t—i>0r—*>0 


(9.55) 


Dehne Aq = A(0), Oq = a(0), /o = /(O), and Q = Q{X) = P{X, 0). Now from equation 
(9.50) it is seen that /d(0) = n — 1. We would denote Qo = Q{Xq), the equations (9.54) 
and (9.55) reduces to 


Xn = 


[n + 1) 


a, a + 1, ~ Qo 

/ Aq 


n+1 ]_ 

^ Cl, Cl + 1, 


foQ 


n—1 

0 


(9.56) 


and 
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Xa = 


( n / Ilf I 1 ) Qo a, a + 1, 

(n — l)vl +/o + 1 t 2 Ao 


n —1 


1 —n 

~ Qo^ 0 , 0 + 1 ,— 


{Q!o(n + 1) 


2n + 2} 



/oQr^ 

Ao 


+ 


1 —n 

ttoQo" 



/oor^ 

'^0 



+ Qo{n - 1 - ao) 


(9.57) 


TABLE: Positive values of Xq by eliminating Qq from equations (9.56) and (9.57) for 
different values of the parameters /o, Ao and ao in various dimensions. 
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/o 

Ao 

ttO 


Positive roots (Xq) 






4D 

5D 6D 

7D 

8D 

lOD 

12D 

16D 

-.033 

.034 .2 6.75507, 
1.14476 

4.32936, 3.15991, 
1.16321 1.1879 

2.44524, 

1.22434 

1.92883, 

1.29065 

.669541, 

.415056 

— 

— 

-.1 

.333 

.2 3.15662, 
1.24003 

1.72737, - 

1.48602 

— 

— 

— 

— 

— 

-.301 

.302 

.001 2.07773, 
1.65625 

— — 

— 

— 

— 

— 

— 

-.5 

.51 

.001 

— — 

— 

— 

— 

— 

— 

.1 

.1 

2 1.667, 

.806578 

1.1946 

.858211 

— 

— 

— 

— 

— 

.1 

.1 

10 .463679 

.465426 .466489 

.466499 

.464947 

.45438 

.429651 

.35278 

.1 

1 

.2 

— — 

— 

— 

— 

— 

— 

.1 

5 

5 

— — 

— 

— 

— 

— 

— 

1 .' 

001 

2 .996033, 

.414324 

.968869, .91379, 
.41436 .414412 

.848544, 

.414489 

.783787, 

.414605 

.669541, 

.415056 

.577489, 

.416214 

— 

10 .001 . 

2 .316014, 

.275327 

.309739, 

.275856 

— 

— 

— 

— 

— 

10 

.1 

5 .089435 

.089441 .089435 

.089409 

.08935 

.089061 

.088347 

.084352 

10 

1 

3 .128192, 

.128966, .19718, 
.41436 

.169063, 

.414412 

.142103, 

.414489 

— 

— 

— 


From the table we conclude the following results: 

(i) For same value of the parameters (/o, Aq, Oo) the possibility of positive real root is 
more in 4D and then it decreases with increase in dimensions. 

(ii) If Aq is positive but close to zero then it is possible to have a naked singularity. 

(iii) For same values of /o and Aq, the naked singularity in higher dimensions will less 
probable as we increase the value the parameter oq. 
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The strength of singularity [Tipler, 1987], which is the measure of its destructive 
capacity, is the most important feature. A singularity is gravitationally strong or simply 
strong if it destroys by crushing or stretching any objects that fall into it. Following 
Clarke and Krolak [1986],we consider the null geodesics affinely parameterized by p and 
terminating at shell focusing singularity r = t = fi = 0. Then it would be a strong 
curvature singularity as defined by Tipler [1980] in equations (8.23) and (8.25). 

It is widely believed that a space-time does not admit analytic extension through a 
singularity if it is a strong curvature singularity in the above sense. Now equations (8.23) 
and (8.25) can be expressed as 


lim = 

p —0 


lim 

/i —0 


nF' ( 

2rn-2pn-2jll R J 


(9.58) 


Using L’Hospital’s rule and using equations (9.50) and (9.51), the equation (9.58) can 
be written as 


lim = 2ra(n — l)AoQoXo(Ao-|-/ oQq ^)^/l + fp 

' ^ [(’^ ~ l)'^o-^oVl + /o — {n — l)Ao(5o — «o/oQo] 


(9.59) 


Thus along the radial null geodesics, the strong curvature condition is satisfied and 
hence it is a strong curvature singularity. So the formation of naked singularity will be 
less probable as we increase the dimension of space-time. 


9.7 Discussion 

In this chapter, we have studied spherical dust collapse in an arbitrary (n -|- 2) dimen¬ 
sional space-time. We have considered both non-marginal and marginally bound cases 
separately. For non-marginal case we have seen that naked singularity may be possible 
for all dimensions (n > 2). However, to get a definite conclusion about naked singularity 
we should study the geodesic equations for non-marginal and marginally bound cases. 






CHAPTER\B Naked Singularities in Higher Dimensional Spherieal Gravitational Collapsel21 


But we can not proceed further due to the presence of the complicated hypergeometric 
functions. Therefore, no dehnite conclusion is possible for non-marginal case. 

On the other hand, for marginally bound case we have dehnitely concluded that naked 
singularity is possible only for n < 3 by studying the existence of radial null geodesic 
through the singularity. Here, we should mention that this result depends sensitively on 
the choice of the initial conditions. In fact, if we do not assume the initial density to 
have an extremum value at the centre (i.e., pi 7 ^ 0 ) the naked singularity will be possible 
in all dimensions. We should mention that the naked singularity described above is only 
a local feature, it is not at all a global aspect i.e., it violates the strong form of CCC. 
The numerical results which favour the formation of naked singularity in any dimension. 


Chapter 10 

Study of Higher Dimensional 
Szekeres’ Cosmological Model 

10.1 Prelude 

Usually, cosmological solutions to Einstein’s field equations are obtained by imposing 
symmetries [Heckmann et al, 1962] on the space-time. One of the reasonable assumptions 
(in an average sense) is the spatial homogeneity. But when we consider cosmological 
phenomena over galactic scale or in smaller scale (detailed structure of the black body 
radiation) then we should drop the assumption of homogeneity i.e., inhomogeneous so¬ 
lutions are useful. 

Szekeres [1975] gave a class of inhomogeneous solutions representing irrotational dust 
for the metric of the form (known as Szekeres metric) 

= dt^ — e^°'dr‘^ — e^^{dx^ + dy^) 

Subsequently, the solutions have been extended by Szafron [1977] and Szafron and Wain- 
wright [1977] for perfect fluid and they studied asymptotic behaviour for different choice 
of the parameters involved. Later Barrow and Stein-Schabes [1984] gave solutions for 
dust model with a cosmological constant and showed the validity of the Cosmic ‘no-hair’ 
Conjecture. 


In this chapter, we give inhomogeneous solutions for (n -|- 2)-dimensional Szekeres 
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space-time with perfect fluid and a cosmological constant. An asymptotic study of par¬ 
ticular solutions are also given. 

10.2 Basic Equations 

The metric for the [n + 2)-dimensional Szekeres space-time is in the form 

n 

ds^ = dt^ — ( 10 . 1 ) 

i=l 

where a and (3 are functions of all the {n + 2 ) space-time variables i.e., 
a = a(f, r, Xi,..., a;„), /? = r, Xi,..., Xn). 

The Einstein equations for the perfect fluid with a cosmological constant is of the form 

Gf,u = + {p + p)Uf,Uy - (10.2) 

where p and p are energy density and isotropic pressure measured by an observer 
moving with with the fluid, A is the cosmological constant and is the fluid flow 
four vector. Since u = d/dt, the flow lines are geodesics and the contracted Bianchi 
identities imply that pressure is a function of t only i.e., p = p{t). As there is no 
restriction on the energy density so p is in general a function of all the (n -|- 2 ) variables 
i.e., p = p{t,r,xi, ...,Xn) and hence no equation of state is imposed. 

Now from the non-vanishing components of the held equations (10.2) for the above 
metric ( 10 . 1 ), we have 

1 • " f 1 

nd/j + -n{n - 1)/?^ - 6“^^ lal. +-{n - l){n - 2)j3l. + {n - 2)a^.^^, -f 

i=l 1 

+ {n - l)PxiXi} + |na'/3' - ^n{n + l)/3'^ - n/3"| = A + p (10.3) 
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1 ... 1 (1 

-n(n + 1)/?^ + n$- -n{n - ^ -(n - l)(n - 2)/?^, + (n - 1)A 


2 = 1 




= A — p 


(10.4) 


+ a + {n — l)aj3 ^ — n{n — 1)0^ + [n — l)f3 + e l{n — l)a' 0 - n{n — 1)0'^— 


n . 

(n - 1)0'} - 


+ -jin - 2)(n - 2,)(3l^ + a:,.:,. + (n - 2)/?^.^. + (n - 


-e + ^(n-l)(n-2)/32,| =A-p (10.5) 

axj {-oixi + I3xi)+ Pxj (ttxi + {n- 2)(5^.) - - (n - 2)j3xiXj =0, {i 0 j) (10.6) 

a0-00-0 = O (10.7) 

-aaxi + pax;^ - - {n - l)/5^, = 0 (10.8) 


ax.P' -Px. = 0 (10.9) 

where dot, dash and subscript stands for partial differentiation with respect to t, r 
and the corresponding variables respectively (e.g., Pxi = §^) with i,j = 1,2, ...,n. 

From equations (10.7) and (10.9) after differentiating with respect to x* and t respec¬ 
tively, we have the integrability condition 


Px.0^ = 0 , 

Thus, \i 0 0 Q we must have Pxi = 0, i = 
consider the following possibilities 

(0 P' 0 0, pi) 0 = 0, 


i = l,2,...,n (10.10) 

: 1, 2,..., n. In the following section we shall 

Pxi = 0 (i = 1,2, ...,n) 


to get solutions of the field equations. 
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10.3 Solutions to the Field Equations 

In this section we shall solve the held equations (10.3)-(10.9) using the above restrictions 
separately. 

Case T. /3' 7 ^ 0 


Here due to the restrictions (3xi = 0, i = 1,2, ...,n we have from the held equations 
(10.7) and (10.9), the form of the metric coefficient as 

= i?(f,r) (10.11) 

and 

e^ = R' + Ru' (10.12) 

Now substituting these forms for the metric coefficient in equation (10.4) we have the 
diherential equations for R and u as 


.■1 .1 n — ^ 

RR+ (n- 1)R^ + -{p{t) - A)R^ = ^/(r) (10.13) 

In 2 

and 

n 

+ =nf{r)-n (10.14) 

with /(r) as arbitrary function of r alone. 

Equation (10.13) can be integrated once to have the hrst integral 


+ f(r) + (10.15) 

n{n +1) it” ^ n J 

where F{r) is another arbitrary function of r (appears due to integration). 

Also from (10.14) the solution for v will be 


e 


= A{r) X] + X] Bi{r)xi + C'(r) 
2 = 1 2 = 1 


— U 


(10.16) 
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with the restriction 


Y, Bf - AAC = /(r) - 1 


(10.17) 


2=1 


for the arbitrary functions A{r), Bi{r), i = 1,2, ...,n and C{r). 
Further, to solve R completely let us consider p as a polynomial in t as 

p{t) = pof-"^ 


(10.18) 


(po and a are positive constants) and we have the general solution for R as 


n+1 

R— = 


according as ^ is an integer, non-integer and a = 2. Here Ci and C 2 are arbitrary 
functions of r and we have chosen 


1 (n + l)po 1, r- - 


It is to be noted that to derive the above solution we have chosen A = 0 = f{r). 
However, for non-zero A (but /(r) = 0) the solution is possible only for a = 0 and 2 as 


n+1 

R— 


C\,Cos{tf’-i^{pii - A} + - A)}, when a^O 

+ whena =2 


( 10 . 20 ) 


with C = +l-+i^. 


Further, if we consider the dust model (i.e., p(t) = 0) then the above solution (10.19) 
simplihes to R(^+^IR oc t. Hence for the usual 4D (i.e., n = 2) the scale factor R grows 
as as in the usual Friedmann model. 
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Now, the physical and kinematical parameters have the following expressions 


n F' + {n + l)Fn' po 
2 R^{R' + RF) F 


( 10 . 21 ) 


RR + (n + l)RRv' + nRR' 
R{R' + RF) 


( 10 . 22 ) 


a2 = 


n 


8{n + l)(n — 1)^ 


2R^-\Rf - 2R!f) + {n- 1){RF' - {n + l)R!F) 
RR^{R' + RF) 


(10.23) 


Case II: /?'=/ja,. = 0, i = l,2,...,n 


In this case from the held equations we have the form of the metric functions 

= R{t) (10.24) 

and 

e" = R{t) p(r, Xi,X 2 ,..., a;„) + p(t, r) (10.25) 

Then as before from the held equation (10.4) we have similar diherential equations in 
R and u as 

•■1 -1 T) — 1 

RR+-{n- 1)R^ + -{p{t) - A)R^ = -^K (10.26) 

and 

n 

E - 2)4 + 2^x,xJ = nK (10.27) 

i=l 

with K, an arbitrary constant. 

Here we take the solution for v in the form 


= PY,xj + Y,Qi^i + S 

i=l i=l 


(10.28) 
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where the arbitrary constants P,Qi {i = 1, 2,n) and S are restricted as before 


5^Q?-4PS = A' 

i=l 

Now to determine the function rj, we have from the held equation (10.6) 


(10.29) 


’'v) ^ Q 
dxidxi 

and then from the held equation (10.5) we have the solution 


(10.30) 


with M(r), nj(r) {i 


n 


n 


e '"rj = u{r) x^ + Vi{r)xi + w{r) 

i=l i=l 

1,2, ...,77.) and w{r) as arbitrary functions. 


(10.31) 


Further, to obtain the function /i we use the following combination of the held equa¬ 
tions namely, 

n 

^GS-G; = (n-l)(A-p(t)) (10,32) 

i=l 

and the resulting diherential equation in fi is 


Rfl + {n — l)Rfi -|- /i 


R + -{p{t)-A)R = g{r) 
n 


with 


g{r) = (n - 1) 


2{uS -|- wP) — ViQi 

i=l 


(10.33) 


(10.34) 


For explicit solution if we choose p{t) as in the previous case (eq.(10.18)) then the 
explicit form for R is same as in equation (10.19) except here Ci and C 2 are arbitrary 
constants and we have chosen iF = A = 0. Similarly, we have the same solutions (10.20) 
for non-zero A. But we note that the diherential equation (10.33) is not solvable for any 
value of n. In fact only for n = 3 (i.e., for hve dimension) we have the complete solution 


pR = + d2t^ 

for a = 2 and ga = | (l + -^1^^) with A = g{r) = 0. 


(10.35) 
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Further, the physical and kinematical parameters have the expressions as 


2nuR — K n 
P=—^ -+ 


2/ii? 


n 


'r]{R‘^i2 — njjiRR) — {n — l)jjfR 
nR^jji + rjR) 


p{t) 


(10.36) 


9 = 


Rfi + n^R + (u + 1)RR7] 
R{ia + r]R) 


(10.37) 


a2 = 


n 


2(n+ 1) 


Rfi — R^ 
R{fi + r]R) 


2 


(10.38) 


Finally, for simple dust case we have the solution for A = 0 in terms of hypergeo metric 
function as 


— to) 




.1 n + 1 3n — 1 
^^^^2’ 2n-2’ 2n-2’ 


{n + 1)2 
4nCi 


/(r)i?"-i] 


(10.39) 


But for non-zero A we can have solution only for u = 3 as 


(10.40) 

However, we can have an integral equation from (10.15) (with p = 0) as 


t - to = 


dR 

\j 7 ( 7 ^^^ + f{f) + 


and we have the following particular solutions : 


(10.41) 


(i) F{r) = 0 


(n + l)/(r) 
2A 


{t - to). 


2A 


n{n -|- 1) 


(ii) /(r) = 0 (i.e., iF = 0 in Case If) 


= JAUMpA Srah 
2A 


{t - to) 


{n -I- 1)A 


2n 
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(iii) A = 0 (n = 3) 


3 f{r) 


(iv) /(r)=A = 0 


ri+l_ Ti \ /TT7—T / \ 

R 2 = (t-to) 


(v) /(r) = F{r) = 0 


^ ^ g(i ^o)^J „(„+!) 


Also for the dust fi has the solution 


fx R = Cie V 3 -)- (726 


^gjr) 

2A 


for n = 3. 


(10.42) 


10.4 Asymptotic Behaviour 

We shall now discuss the asymptotic behaviour of the solutions presented in the previous 
section for both perfect fluid and dust model separately. The co-ordinates vary over the 
range : to < t < oo; — oo < r < oo] — oo < Xi < oo, i = 1,2, 


10.4.1 Perfect Fluid Model 

As p > 0, p 7 ^ 0 so we must have | < gi < 1. We shall hrst consider the case when 
a = 2. For large t {Case I i.e., (3' ^ 0) 

Tl 

p ~ 

6 ^ t ^ 
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^2 +-2 
(J ~ I 

For Case II, we have similar behaviour for large t together with ji ~ where qi is 
the value of qi for n = 3 and also we have \ < q\ < 1. Thus as f > cx), {p, p) fall off 
faster compare to {9, a), while the scale factor R (and p) gradually increases with time. 
So the model approaches isotropy along fluid world line as f —>• cx). 


10.4.2 Dust Model 

For the dust case with non-zero A (and n = 3) we have for large t, 

p = R ^ 

p = Po, a constant 

9 = 9q, a constant 

a = (To, a constant 

We note that for Case I, po = ao = 0 while for Case II we have po 7 ^ 0, (Tq 7 ^ 0. Thus 
universe will behave locally like de-sitter model in Case I though the global geometry 
will be different. However, in Case II the universe will not isotropize and for large t the 
measure of anisotropy and expansion scalar become hnite constant. 

10.5 Discussion 

In this chapter, we have investigated cosmological solutions for [n + 2)-dimensional 
Szekeres form of metric with perfect fluid (or dust) as the matter distribution. We can 
classify the solutions in two categories namely, (i) /3' 7 ^ 0 and (ii) j3' = 0. The first 
set of solutions are known as quasi-spherical solution while second class of solutions are 
termed as quasi-cylindrical type of solutions. 

However, if we assume the arbitrary functions A{r), Bi{r) and C{r) to have the con¬ 
stant values namely A{r) = C{r) = ^ and Bi{r) = 0 (for all i = l, 2 ,...,n) i.e., if we 
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have chosen the arbitrary function /(r) to be zero then using the transformation 

xi = SinOriSinOn-i . Sin92Cot^6i 

X2 = Cos6nSin6n-i . Sin62Cot\6i 

X3 = C os6n-iSin6 n-2 . Sin62Cot\6i 

Xn-i = Cos93Sin92Cot^9i 
Xn = Cos 92 Cot\ 9 i 

the Szekeres metric (10.1) with the solutions (10.11), (10.12), (10.16) reduces to the 
{n + 2)-dimensional spherically symmetric metric 

ds^ = df - R'^dr'^ - R^dnl 

It is to be noted that if the arbitrary functions A, Bi,C depend on r, then we can 

not get the above spherical form of the space-time. So these arbitrary functions play an 

important roll to identify the nature of the space-time. 

Finally, the study of asymptotic behaviour shows that some of the solutions will be¬ 
come isotropic at late time while there are solutions for which shear scalar remains 
constant throughout the evolution. Thus Cosmic ‘no-hair’ Conjecture is not valid for 
all solutions. This violation of Cosmic ‘no-hair’ Conjecture is not unusual because the 
Szekeres metric may not have always non-positive 3-space curvature scalar. Also in this 
higher dimensional Szekeres space-time we have solutions which expands as de Sitter in 
some directions but not in other directions. 



Chapter 11 


Shell Focusing Singularities in 
Quasi-Spherical Gravitational 
Collapse 

11.1 Prelude 

An extensive study [Joshi et al, 1994, 1999; Lake. 1992; Ori et al, 1987; Harada, 1998] of 
gravitational collapse has been carried out of Tolman-Bondi-Lemaitre (TBL) spherically 
symmetric space-times containing irrotational dust. Due to simplihcations introduced 
by the spherical symmetry several generalizations of this model have been considered. 
A general conclusion from these studies is that a central curvature singularity forms 
but its local or global visibility depends on the initial data. Also the study of higher¬ 
dimensional spherical collapse reveals the interesting feature that visibility of singularity 
is impossible in space-times with more than hve dimensions with proper choice of regular 
initial data [Goswami et al, 2002]. 

By contrast, there is very little progress in studying non-spherical collapse. The quasi- 
spherical dust collapse models, given by Szekeres metric [Szekeres, 1975] were analyzed 
by Szekeres himself [1975], Joshi and Krolak [1996], Deshingkar, Joshi and Jhingan [1998] 
and extensively by Goncalves [2001]. The solutions for dust and a non-zero cosmological 
constant were found by Barrow and Stein-Schabes [1984]. In this work, we study the 
gravitational collapse in the recently generalized (n -|- 2)-dimensional Szekeres metric. 
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As in four-dimensional space-time, this higher dimensional model does not admit any 
Killing vector and the description quasi-spherical arises because it has invariant family 
of spherical hypersurfaces. 

In this chapter, we study the occurrence and nature of naked singularities for a dust 
model with non-zero cosmological constant in (n -I- 2)-dimensional Szekeres space-times 
(which possess no Killing vectors) for n >2. We hnd that central shell-focusing singular¬ 
ities may be locally naked in higher dimensions but depend sensitively on the choice of 
initial data. In fact, the nature of the initial density determines the possibility of naked 
singularity in space-times with more than hve dimensions. The results are similar to the 
collapse in spherically symmetric Tolman-Bondi-Lemaitre space-times. We also extend 
our study for global characteristic of the singularity by studying both null and time like 
geodesic originated from the singularity using the formalism of Barve et al [1999]. The 
nature of the central shell focusing singularity so formed is analyzed by studying both 
the radial null and time-like geodesic originated from it. 


11.2 Basic Equations and Regularity Conditions 


The dust solutions have been obtained for (n -|- 2)-dimensional Szekeres’ space-time 
metric for which the line element is given in equation (10.1). Under the assumption that 
jS' 7 ^ 0, the explicit form of the metric coefficients are also given in equations (10.2) and 
(10.3), where the form of u is given in equations (10.16) and (10.17) and R satisfied the 
differential equation 




f{r) + 


F{r) 


+ 


2A 

n{n -f 1) 


R^ . 


( 11 . 1 ) 


Here A is cosmological constant, /(r) and F{r) are arbitrary functions of r alone. 

In the subsequent discussion we shall restrict ourselves to the quasi-spherical space- 
time which is characterized by the r dependence of the function v (i.e., F ^ R). From 
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the Einstein field eqnations we have an expression for energy density for the dust, as 


, , n F' + (n + l)Fz/' 

p{t,r,xi, ...,Xn) = - 


( 11 . 2 ) 


2 + RC) 

Thus a singularity will occur when either (i) F = 0 i.e., (3 = —oo or (ii) a = —oo. 
Using the standard terminology for spherical collapse, the first case corresponds to a 
shell-focusing singularity, while the second case gives rise to a shell-crossing singularity. 
As in a TBL space-time, the shell-crossing singularities are gravitationally weak, and 
we shall not consider them any further here. Hence, we shall restrict ourselves to the 
situation with a > —oo. 

Suppose that t = U is the initial hypersurface from which the collapse develops. For 
the initial data we assume that i?(U,r) is a monotonically increasing function of r. So, 
without any loss of generality, we can label the dust shells by the choice R(ti,r) = r. 
As a result, the expression for the initial density distribution is given by 


piir, xi,..., Xn) = p(U, r, xi,..., x„) 


n F' + {n + l)Fz/' 
2 r"'(l -|- rC) 


(11.3) 


If we started the collapse from a regular initial hypersurface the function pi must be 
non-singular (and also positive for a physically realistic model). Furthermore, in order 
for the space-time to be locally fiat near r = 0, we must have /(r) ^ 0 as r —> 0. Then, 
from equation (11.1), the boundedness of as r —>• 0 demands that F{r) O(r-) 
where m > n — 1. But, for small r, piir) ~ ^ and consequently, for regular Piir) near 
r = 0, we must have F{r) ~ 0(r"+^). Thus, starting with a regular initial hypersurface, 
we can express F{r) and Pi{r) by power series near r = 0 as 


and 


F(r) = ^ F,-(11.4) 

j=0 


OO 

f(^) = 5 ^ Pi rC 

i=o 


(11.5) 
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As v' appears in the expression for the density as well as in the metric coefficient, we 
can write 


E 

i=-i 


u'ir) = Vj 


( 11 . 6 ) 


where > —1. 

Now, using these series expansions in equation (11.2) we have the following relations 
between the coefficients, 

n , -I- 1 i n . ( 1 \ 

Fi, 


n(n +1) „ n f ^ 1 

Po-^—Fo, ft-2("+i + TTy; 


n 


P3 = 


n 

P2= 2 


n + 1 + 


n + 1 + 


Fiv. 




2F2Z/0 (1 + - Z /2 


l + I/_l/ (1 + Z/_ (1 j/_ J3 




(11.7) 




OR 


Po = 


n 


a 

.^0 


+ {n + l)Fo 


n 


Pi = 


2F2 I , N ^1 

— + i (" + 1) - 73 




n 

P2= 2 


3 F 3 

1^0 


2z/i 




P3 = 


n 


4 F 4 


+ < (n + 1) - 


3z/i 




F 3 + 2 


Z2t 


Z2n 




Z2n 


7^2 + 


/ 2z/it/2 _ ^ _ 

V ^^0 ^0 


z/r 


Fi 


( 11 , 8 ) 




according as ^ > —1 or z/_^ = —1. 

Finally, if we assume that the density gradient is negative and falls off rapidly to zero 
near the centre then we must have pi = 0 and p 2 < 0. Consequently, we have the 
restrictions that Fi = 0 and F 2 < 0 for z/_j > —1. 
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11.3 Formation of Singularity and its Nature 

In order to form a singularity from the gravitational collapse of dust we first require that 
all portions of the dust cloud are collapsing i.e., R < 0. Let us dehne tsf{r) and fsc(r) 
as the time for shell-focusing and shell-crossing singularities to occur occurring at radial 
coordinate r. This gives the relations [Szekeres, 1975] 

R{tsf{r),r) = 0 (11.9) 

and 

R'{tsc, r) + R{tsc, r)iy'{r, Xi, X 2 ,Xn) = 0 . ( 11 . 10 ) 

Note that ‘tsc may also depend on xi,X 2 , ■■■,Xn- 
As mentioned earlier, the shell-crossing singularity is not of much physical interest so 
we shall just consider the shell-focusing singularity for the following two cases: 


(i) Marginally bound case : /(r) = 0 with A 7 ^ 0 

In this case equation (11.1) can easily be integrated to give 


t — ti-\- 


2n 


{n + 1 )A 


Sink 


-1 


2Ar”+^ 


n{n + l)F{r) 


Sink 


-1 


n{n -|- l)F(r) 


(1 


. 11 ) 


Also, from equations ( 11 . 1 ) and (11.11), we have 


tsf{r) =ti + 


2n 


{n 1)A 


Sink 


-1 


2Ar"+i 


n(n -|- l)F(r) 


( 11 . 12 ) 


(ii) Non-marginally hound case with time symmetry : /(r) 7 ^ 0, R{ti,r) = 0, A = 0 
Here, the solution of equation (11.1) gives t as a function of r : 


2 

(n -h 1 )a/F 


r (6 -1-1) 
r (6 + i/2) 


n+l n+1 

r 2 — R 2 





n —1 


t — ti-\- 


(11.13) 
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and so the time of the shell-focusing singularity is given by 


tsf{r) 


^ ^ 2y/n r(6 -|-1) n+i 


(11.14) 


where 2 F 1 is the usual hypergeometric function with 6 = | -|- . However, for the 

hve-dimensional space-time (n = 3), R has the particularly simple form 


- EA(t _ i.)2, 


and therefore 


tsf{r) 


ti 


■/W)' 


(11.15) 


(11.16) 


11.4 Formation of Trapped Surfaces 

The event horizon of observers at inhnity plays an important role in the nature of the 
singularity. However, due to the complexity of the calculation, we shall consider a 
trapped surface which is a compact space-time 2-surface whose normals on both sides 
are future-pointing converging null-geodesic families. In particular, if (r = constant, t = 
constant) the 2-surface Sr,t is a trapped surface then it, and its entire future development, 
lie behind the event horizon provided the density falls off fast enough at inhnity. 

If tahir) is the instant of the formation of apparent horizon then we must have [Szek- 
eres, 1975] 

2AR"''^^(tah{r),r) - n{n -P l)R^~^(tah{r),r) + n{n -f l)F(r) = 0. (11-17) 

Thus, for the above two cases, the explicit expressions for tah are the following: 




2n 


{n 1)A 


Sink 



2Ar^+i 


nin 


+ 1)F{7 


Sinh-^ 


/ / 2Ai^»+i(U,r) y 
n(n + l)F(r) J 

(11.18) 
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for /(r) = 0, A 7 ^ 0 and 


, , 2 [v^rffe + l) ^ . ^,177 F 1 

(n + l)/F [ r(6 + l/2) 2 U 2 > > 


(11.19) 


for /(r) 7 ^ 0, R{ti, r) = 0, A = 0 and for n = 3 we get 


tah{r) =ti + r^^ -1. (11.20) 

From the expressions for tsf{r) and tahif), we note that the shell-focnsing singnlarity 
that appears at r > 0 is in the fntnre of the apparent horizon in both cases. Bnt since 
we are interested in the central shell-focnsing singnlarity (at r = 0 ) we reqnire the time 
of occnrrence of central shell-focnsing singnlarity to(= tsf{0)). From eqnations (11.12), 
(11.14), (11.16) and (11.17)-(11.20), taking the limit as r —0, we have that 


tahir) - to = 


1 

(n + 


Fir + 



(4A + 3n(n + l)Fo)Ff \ 
4Fo(2A + nin + l)Fo) / 


r^+ 



(4A + 3n(n + l)Fo)A’iF 2 (32A^ + 40n(n + l)AFo + 15?7,^(n + 1 )^ + F^)Ff 
4Fo( 2A + n(n + l)Fo) ^ 24F2(2A + n(n + l)Fo)' 


r 


3 


in 


l)^c 


n — 2 
n—1 


n + 1 
^pT-n-l 


n — 1 


2n 

Fir"-i + 


n — 1 


+2- 


n — 2 Fi 
2 (n - 1 ) Fo 


f n—1 


( 11 . 21 ) 


for fir) = 0, A 7 ^ 0 and 


tahir) - to = 


0F F(6 + l) 
(n + l)Fo"/^F(6 + l/2) 


Fir + 




9 F 1 F 2 

4Fo 





2 

n — 2 

in + l)Ff-^ 


77 + 1 1 277 

For"-i H-Fir"-i + ... 

n — 1 


( 11 . 22 ) 


for fir) 7 ^ 0, Riti,r) = 0, A = 0 and 
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for n = 3 (with /(r) 7 ^ 0, r) = 0, A = 0) 

1 


tahir) - to = 


2K 


3/2 


Fi r + I F 2 + Fq - 


3F2 

4 ^ 


+ 


(11.23) 


Here, to is the time at which the singularity is formed at r = 0, and tahir) is the instant 
at which a trapped surface is formed at a distance r. Therefore, if tahir) > to, a trapped 
surface will form later than the epoch at which any light signal from the singularity 
can reach an observer. Hence the necessary condition for a naked singularity to form is 
tahir) > to; while that for a black hole to form is tahir) < to. It is to be noted that this 
criteria for naked singularity is purely local. Hence, in the present problem it possible 
to have local naked singularity or a black hole form under the conditions shown in the 
Table I, for > —1. 


TABLE-I 



Naked Singularity 

Black Hole 

Marginally bound case: 
/(r) = 0 ,A ^0 

(i) Pi < 0 , Vn 

(ii) Pi = 0,p2 < 0,n = 2 
(hi) Pi = 0,p2 < 0 ,n = 3 , 

< -2-PoV^ + A 

(iv) Pi = 0,P2 = 0,P3 < 0 , 
n = 2, 

F, < -2FS'Qi + 4- 

(i) Pi = 0,p2 < 0 ,n = 3 , 

F2 > A + A 

(ii) Pi = 0,p2 < 0 ,n > 4 

(hi) Pi = 0, p2 = 0, p3 < 0, n = 2, 

F, > -2FPQl + 4r 

(iv) Pi = 0 , P2 = 0 , P3 < 0 , n > 3 

(v) Pi = 0,p2 = 0, ...,pj = 0, 

Pj+i < 0 , j > 3 ,Vn 

Non-marginally bound 
case with time 
symmetry: /(r) 7^ 0, 
RiU, r) = 0 , A = 0 

(i) Pi < 0 , Vn 

(ii) Pi = 0,p2 < 0,n = 2 
(hi) Pi = 0,p2 < 0 ,n = 3 , 

F2 < -F^ 

(iv) Pi = 0,P2 = 0,P3 < 0 , 
n = 2 ,F,<-f;;Fp 

(i) Pi = 0 , p2 < 0 , n = 3 , F2 > -F^ 

(ii) Pi = 0 , p2 < 0 , n > 4 

(ih) Pi = 0, p2 = 0, p3 < 0, n = 2, 

^3 > -^^0 

(iv) Pi = 0 , P2 = 0 , P3 < 0 , n > 3 

(v) Pi = 0,p2 = 0, ...,pj = 0, 

Pj+i < 0 , j > 3 ,Vn 


However, if we relax the condition on the density gradient (namely the condition that 
Pi = 0 , P 2 < 0 i.e., Fi = 0 , F 2 < 0 for > — 1 ) then it is possible to have a naked 
singularity in all dimensions. These results are identical to those obtained in the pres¬ 
ence of spherical symmetry. Hence the local nature of singularity is not affected by the 
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geometry of the space-time with respect to whether it is spherical or quasi-spherical. 
However, we note that the Szekeres space-times retain certain special geometrical prop¬ 
erties shared by spherical space-times in all dimensions, in particular they do not allow 
gravitational radiation to be present in the space-time. 


11.5 Radial Null Geodesics 


Now we shall consider A = 0, /(r) = 0 and initial time = 0 so that i?(0, r) = 0. So 
the solution for R of the equation (11.1) becomes 


R = 


n+1 77/ “1“ 1 / _ ,—r 

r 2 - ^/F{r) t 


n+1 


(11.24) 


Now with the help of equations (11.9) and (11.24), the singularity time can be ex¬ 
pressed as 


tsf{r) = 


2 r 2 


n + l 

Near r = 0, above singularity curve can be approximately written as 


(11.25) 


WW =to- ^"3/2 ( 11 - 26 ) 

(^+ 1)^0 

where m (> 1) is an integer and Fm is the first non-vanishing term beyond Fo in the 
series form (11.4) and R = 

The equation of the outgoing radial null geodesic (ORNG) which passes through the 
central singularity in the past is taken as (near r = 0) 


toRNG — to + (11.27) 

to leading order in t-r plane with a > 0,^ >0. Now to visualize the singularity the 
time t in the geodesic (11.27) should be less than ts/(r) in equation (11.26). Hence 
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comparing these times we have the restrictions 

F 

^ > m and a < -——^ (11.28) 

(n + 

Further, for the metric (10.1) an outgoing radial null geodesic should satisfy 

^ = R' + Ru' (11.29) 

dr 

We shall now examine the feasibility of the null geodesic starting from the singularity 
with the above restrictions for the following two cases namely, (i) .^ > m and (ii) ^ = m. 


When ^ > m then near r = 0 the solution for R in (11.24) simplihes to 


R = 


W \ n+l 

^ m 


2Fn 


f n+l 


+1 


(11.30) 


Now combining (11.27) and (11.30) in equation (11.29) we get (upto leading order 


m r 


a ^ = ( z/_i + 1 + 


2 m 


n + 1 / \ 2Fr 


Fm ^ ^ 

f n+l 


(11.31) 


which implies 




n + l 


2Fn 


(1 + iS) 

We note that if = 0 then the above results are same as TBL model in eq. (9.38), 
so we restrict ourselves to 7 ^ 0 . As ^ > m so from (11.32) 


n + l 

m < - 

n — 1 

Since m is an integer so we have the following possibilities: 

n = 2 {AD space-time) : m = 1, 2 , | 

n > 3 {5D or higher dim. space-time) : m = 1 , .^ = 1 -|— 
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Now for z/_j > —1, if we assume that the initial density gradient falls off rapidly to 
zero at the centre then we must have pi = 0, p 2 < 0 i.e., Fi = 0, F 2 < 0. Hence the least 
value of m should be 2. Thus for a self consistent outgoing radial null geodesic starting 
from the singularity we must have n = 2 only. So the singularity will be naked only in 
four dimension. However, if we relax the restriction on the initial density and simply 
assume that the initial density gradually decreases as we go away from the centre (i.e., 
Pi < 0) then Fi < 0 and m = 1 is the possible solution. Therefore, naked singularity 
is possible in any dimension. However for = —1, the assumption about the initial 
density gradient (pi = 0, p2 < 0) does not implying Fi = 0. Hence here m = 1 is a 
possible solution for all dimension with Fi < 0. Therefore naked singularity may occur 
in any dimension for this choice of the parameter z/_i. 

Now we shall study the second case namely ^ = m. Here also for = 0 we have 
identical results as in TBL model ((9.39), (9.40)). For 7 ^ 0, using the same procedure 
as above we get 

n + 1 


and 

{^-1 + + 1 ) a 

(11.33) 

Hence integral solution for m is possible only for n = 2 and 3. So consistent outgoing 
radial null geodesic through the central shell focusing singularity is possible only upto 
hve dimension i.e., we cannot have naked singularity for higher dimensional space-time 
greater than hve. It is to be noted that this conclusion is not affected by the value of 
the parameter 

Now we shall examine whether the restriction (11.28) for a is consistent with the 
expression a in equation (11.33). In fact equation (11.33) takes the form 


a = 


1 

m 


n -f 1 




2Fn 


1+n 


2 1 

2^hm = — [—{n 1)6 — C]~™ [(^^-i + ^)C + (^t- + l)(^_i + 1)^] 


(11.34) 
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with the transformation 

a = bF^-\ = (11.35) 

Since eqnation (11.34) is a real valned eqnation of b, so we must have 

(n + 1)6 + C < 0 , 

which using (11.35) gives us the restriction on a in equation (11.28). Hence the geodesic 
(11.27) will have consistent solution for a and m. So the above conclusion regarding the 
formation of naked singularity is justihed. Further, introducing the variable ej) by the 
relation 

0 = -(n + l)6-C (11.36) 

We have seen from equation (11.34) 


4r-'(C+ •#>)”+' = PC - (n- +1).#.]”+* (11.37) 

with the restriction 0 < 0 < 

Finally, we investigate whether there is only one null geodesic emanating from the 
singularity or it is possible to have an entire family of geodesic through the singularity. 
So we write the equation for the outgoing radial null geodesics to next order as 

f = to + 0.^^ + (11.38) 

where as before a, h, a are all positive. Then proceeding in the same way we have 
(i) when ^ > m, a and ^ have the same expressions as in equation (11.32) while 


a = 


m(l 


n) 


1 -|- n 


+ 1 


(11.39) 


and 


h 


a 


-a)^/Fo 


1 — n 

Fm \ 


(11.40) 
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Fig 11.3 Fig 11.4 

In Figs 11.1 and 11.2, we have shown the variation of (0 < (/> < —C) for the variations of 

—1) and C(< 0) for four and five dimensions (i.e., n = 2 and 3) respectively. However 
for v_, = —1, the dependence of 6 over C has been presented in Figs 11.3 and 11.4 for n = 2 
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We note that since here m can have values 1 and 2 for n = 2 and m = 1 for n > 3 so 
a is always positive. But for h to be positive dehnite, is restricted within the range 
— 1 < < —a i.e., is negative dehnite. 

(ii) for ^ = m, a and ^ have the same expressions as in equation (11.33) and a is 
obtained from the equation 



with expressions for b and ( from equation (11.35), while the other constant h is to¬ 
tally arbitrary. Thus it is possible to have an entire family of outgoing null geodesics 
terminated in the past at the singularity, provided a obtained from equation (11.41) is 
positive dehnite. 

11.6 Radial Time-Like Geodesics 

We shall now investigate whether it is possible to have any outgoing time-like geodesic 
originated from the singularity. For simplicity of calculation we shall consider only 
outgoing radial time-like geodesic (ORTG). Let us denote by a unit tangent 

vector held to a ORTG with r an affine parameter along the geodesic. Hence from the 
geodesic equation we have [Goncalves, 2001] 


{K'^y + 2aK^K'^ + RT" {K^)' + a' = 0 


(11.42) 


with 



(11.43) 


In order to satisfy the above two equations the simplest choice for (iF*, K'^) is 


R* = ±1, R’' = 0 . 


(11.44) 


and so we have the solution 


t — to = ±(r — To), r = ro = constant. 


(11.45) 
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Here + (or —) sign corresponds to ORTG (or ingoing RTG) and tq is the proper time 
at which radial time-like geodesic passes through the central singularity. 

Now similar to null geodesic let us choose the radial time-like geodesic near the sin¬ 
gularity to be of the form (to leading order) 

toRTcir) =to + crP (11.46) 

with c and p as positive constants. Further, consistent with equations (11.44) and 
(11.45) we assume 

K^{t,r) = Ait-to^r^ (11.47) 

where A (> 0), A and 6 are constants. Hence we have 

K^'itoRTGA) = ^ q = Xp + 6. (11.48) 

Also from the geodesic equation (11.46) using (11.48) we get 


= cprP-^ = — = + {R + RCy 

dr A’’ 


(11.49) 


Using the solution for R near the singularity in equation (11.49) and equating equal 
powers of r we have 

(i)p=l-g, if -^<q<l. 


(ii) p= 1 + 


^ C-2+}+^)_ 


2m _ _^ 

,1+1 ) ^ TTT^ 


(1+^) V 2Fo; 


q < 


2m 

n+1 




■4-^+ -ft + l + 


1/2 


if ^ = -^1 


Therefore equation (11.46) has consistent solution for c and p depending on the pa¬ 
rameters involved. So it is possible to have outgoing radial time-like geodesic originated 
from the singularity. 
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11.7 Local Visibility 

We have shown the existence of ontgoing both nnll and time-like geodesics which were 
originated in the past from the singnlarity. We shall now examine whether these 
geodesics are visible to any non-space like observer. Now for visibility of the singn¬ 
larity (local nakedness) any fnture directed geodesics (nnll or time-like) starting from 
the singnlarity shonld be ontside the domain of dependence of any trapped snrfaces 
both before and at the time of formation of apparent horizon (AH) which is the outer 
boundary of the trapped surface. For A = 0 and /(r) = 0 the time of formation of 
apparent horizon tah{r) satisfies (from eq. (11.17)) 

R^-\tah{r),r) = F{r) (11.50) 

The apparent horizon is given by the curve 


tah{r) =tQ-- 


in 




3/2 


{Frr 


+ 0(r™+i)) 


2 ^ 


n -|- 1 


/ n+1 _ . n+1 I 1 > \ , . 

+0(r^+^)j (11.51) 


As the apparent horizon and the singularity curve form at the same time at r = 0, so 
the visibility of the singularity is determined by the relative slopes of the curve tah{f') 
and the curve for outgoing radial non-space like geodesics (denoted by toRcir)). Hence 
the necessary and sufficient condition for singularity to be at least locally naked is that 


lim 
r —>• 0 


{( 


^)/( 




(11.52) 


Now we shall examine this condition for both ORNG and ORTG separately. 


11.7.1 Outgoing Radial Null Geodesic 


In this case the ratio of the slopes is 





CHAPTER m Shell Focusing Singularities in Quasi-Spherical Gravitational Collapse 149 


( dtqhjr) \ 
^ dr ^ 



dr 



So we have the following possibilities: 

(i) If m > then the above ratio of the slope is always negative, hence the condition 


is violated. So we always get black hole. 

(ii) If m < then the above ratio of the slopes approaches to +cxd as r ^ O’*" for 


m < S, and while the limit will be hnite and greater than unity for m = ^ provided 
® Thus the singularity will be locally naked for any dimension if m = 1 

and a < — while naked singularity appears only for four dimension if m = 2 

and a < — 

(iii) If m = then we have shown that naked singularity appears only upto hve 
dimension and the singularity will be locally visible for 



or 


Fm < -2Fo" - a{n + m = ^. 


11.7.2 Outgoing Radial Time-Like Geodesic 


In this case the ratio of the slopes is 


( dtah(r) \ 
^ dr ) 



dtoRTci'^) 

dr 



(11.54) 
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We shall now study the following possibilities for local visibility: 

(i) If m > then as before only black hole appears. 


(ii) If m < we have the conclusion as in ORNG except that instead of restricting 

a, we have here c < -%?. 

’ (n+l)p p3/2 

(iii) Similar is the situation for m = 

Hence the conditions for local visibility of naked singularity is consistent for both 
time-like and null geodesic. 


11.8 Geodesics and the Nature of Singularity 

Here we consider A = 0 and /(r) = 0. So the solution of R is given in equation (11.24). 
We now follow the geodesic analysis of Joshi and Dwivedi [1993] for TBL model. For 
simplicity of calculation we introduce the following functions: 


X = ^ 




rF' 


tQ' 

^=Q 


C 


(n-1) 


Q = e~ 


1_«_ 

O _ n+l 

U — (a-l)(„+l) 

r 2 


(11.55) 


where the constant a is restricted by a > 1. In order to determine the nature of the 
singularity we examine whether it is possible to have outgoing null geodesics which 
are terminated in the past at the central singularity r = 0. Suppose this occurs at 
singularity time t = to at which R{to, 0) = 0. In order to decide this we start with radial 
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null geodesic, given by 

dt , R’Q - RQ' 
dr ~ Q 

Next, we introduce the notation u = r“, so we write 


(11.56) 


where U (X, u) 


dR 

du 


U{X,u), 


(11.57) 



with 77 = ^ X + ^ . 


We now study the limiting behaviour of the function X as we approach the singularity 
at 7? = 0, M = 0 along the radial null geodesic identified above. If we denote the limiting 
value by Xq then 


Xo = lim f = Urn g = Urn U{X,u) =f/(Xo,0). 

7?—^-Om—*>0 7?—*>0 7?—>0m—>-0 

Thus if this polynomial equation in Xq has at least one positive real root then it 
is possible to have a radial null geodesic outgoing from the central singularity. More 
explicitly, Xq is the root of the following equation in X: 



eo 


n + 1 


X" + 


eo 


n + 1 


hoj ^ 


X"^ - ©oX"^ + ^©0 = 0 (11.59) 


where suffix ‘o’ for the variables rj, ( and © stands for the values of these quantities 
at r = 0. As exact analytic solution is not possible for X so we study the roots by 
numerical methods. Table II shows the dependence of the nature of the roots on the 
variation of the parameters involved. 
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Table II: Positive roots (Xq) of the eqn. (11.59) for different values of the parameters 
namely a, Co) ©o and in different dimensions, n. 


a 

ho 

^0 

Co 

©0 



Positive roots (Xq) 









4D 

5D 

6D 

7D 

8D 

lOD 

14D 

1 

-6 

.05 

.01 

-5 

.00997 

.09864 

.21131 

.30963 

.38979 

.50797 

.64774 

2 

2 

5 

3 

2 

3.3019, 

3 

1.732, 

1.633 

1.4422, 

1.3195 

1.1968, 

1.3161 

1.2457, 

1.1345 

1.1699, 

1.075 

1.105, 

1.0334 

3 

2 

5 

3 

2 

— 

— 

— 

— 

— 

— 

— 

5 

4 

.1 

3 

2 

1.7068, 

.77462 

1.2678, 

.80678 

1.1613, 

.82742 

1.1148, 

.84873 

1.0889, 

.86583 

1.0612, 

.89102 

1.0377, 

.92113 

1 

0 

9 

1 

1 

2.035 

1.6116 

1.5163 

1.5028 

1.5511 

— 

— 

1 

0 

4 

1 

1 

14.8239 

— 

— 

— 

— 

— 

— 

1 

.1 

1 

1 

1 

— 

— 

— 

— 

— 

— 

— 

1 

1 

1 

1 

1 

1.3104, 

1 

1.1547, 

1 

1.0934, 

1 

1.0627, 

1 

1.045, 

1 

1.0265, 

1 

1.0124, 

1 

4 

0 

5 

.1 

4 

— 

— 

— 

— 

— 

— 

— 

1 

0 

1 

1 

9 

2.25, 

1.5963 

2.5336, 

1.1725 

2.2120, 

1.1045 

1.9675, 

1.0752 

1.7967, 

1.0588 

1.5829, 

1.041 

1.3756, 

1.0256 

1 

0 

1 

5 

1 

— 

— 

— 

— 

— 

— 

— 

1 

0 

1 

.1 

1 

.90297, 

.10733 

.87007, 

.36226 

.86041, 

.54185 

.85671, 

.66031 

.85235, 

.74455 

— 

— 


From the dehnition of the functions in equation (11.55) it is clear that the parameters 
a,^o,(o are always positive while 770 , ©o can have positive as well as negative values. 
From the Table we see that as we decrease the value of a, keeping the other parameters 
hxed, then the formation of a naked singularity is more probable in higher dimensions. 
The situation for (^o is similar. But for ^0 and ijq the condition is reversed: black hole 
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formation becomes more probable as we increase the dimension of the space-time if we 
decrease the value of (or ho) and keep the other constants hxed. However, no such 
conclusions can be drawn about the variation of ©q. 


11.9 Strength of the Naked Singularity 

A singularity is called gravitationally strong, or simply strong, if it destroys by crushing 
or tidally stretching to zero volume all objects that fall into it; it is called weak if no 
object that falls into the singularity is destroyed in this way. A precise characterization 
of Tipler strong [Tipler, 1987] singularities has been given by Clarke and Krolak [1986], 
who proposed the strong focusing condition. A sufficient condition for a strong curvature 
singularity is that, for at least one non space-like geodesic with affine parameter A ^ 0 
on approach to the singularity, we must have 


lim X^RijK^K^ > 0 
A ^ 0 


(11.60) 


where 77* = ^ is the tangent vector to the radial null geodesic. 

Our purpose here to investigate the above condition along future-directed radial null 
geodesics that emanate from the naked singularity. Now equation (11.60) can be ex¬ 
pressed as (using L’Hopital’s rule) 


lim X^RijK^K^ = <o(Ho ^oXo)(go {n+i)m) 
A^O 2^0 

where = 77(Xo, 0), iVo = iV(Xo, 0). 

The singularity is gravitationally strong in the sense of Tipler if 


(11.61) 


) (’^ + 1)00 
4o - (n + l)riQ > max < 0,-— 


or 


,^0 - (’^ + l)?7o < { 0, 


^0^ 


(n -h 1)00 

n+1 
^0^ 
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If the above condition is not satisfied for the values of the parameters then 
JjXttI/ 

„ X^RijK^K^ < 0 and the singularity may or may not be Tipler strong. 


11.10 Discussion 

In this chapter, we have studied gravitational collapse in [n + 2)-dimensional space-time 
using a higher-dimensional generalization of the quasi-spherical Szekeres metrics with 
non-zero cosmological constant. We have examined the local nature of the central shell- 


focusing singularity by a comparative study of the time of the formation of trapped 


surface and the time of formation of the central shell-focusing singularity. If we assume 
the initial density gradient falls off rapidly and vanishes at r = 0, then naked singularity 
formation is possible only up to space-time dimension hve for > — 1. However, if we 
drop the above restriction on the initial density distribution then the Cosmic Censorship 
Conjecture may be violated in any dimension (n > 2). Following the approach of Barve 
et al [1999] for null geodesic we have shown that if we choose the parameter to be 
greater than —1, then it is possible to have a class of outgoing radial null geodesic for 
four and hve dimensions only with the restriction that initial density falls off rapidly 
to the centre. However for = —1, naked singularity is possible in all dimensions 
irrespective of the assumption on the initial density. Thus we deduce that the nature of 
the central singularity depends sensitively in these metrics on the choice of the initial 
data (particularly on the choice of initial density prohle). This is also conhrmed by 
our geodesic study following the approach of Joshi and Dwivedi [1993] where we have 
shown numerically that the nature of singularity depends on the value of the dehning 
parameters at r = 0. Finally, we have examined the strength of the naked singularity 
using the criterion introduced by Tipler [1987]. We found that the naked singularity 
will be a strong curvature singularity depending on the appropriate choice of the value 
of the parameters at r = 0. On the other hand, for time-like geodesic the results are 
very similar to the study of null geodesic. Regarding local visibility we have consistent 
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results for both ORNG and ORTG with some restrictions on the parameters involved in 
the equation of the geodesics. Therefore we conclude that it is possible to have at least 
local naked singularity for the given higher dimensional non-spherical space-time. These 
investigations provide some insights into the phenomenon of gravitational collapse in a 
situation without any imposed Killing symmetries. However, the collapses are special in 
other senses which permit exact solutions to be found. In particular, there is an absence 
of gravitational radiation in these space-times [Bonnor, 1976]. An investigation of its 
role is a challenge for future analytic and computational investigations. 


Chapter 12 

Shell Crossing Singularities in 
Quasi-Spherical Szekeres Models 

12.1 Prelude 

In the study of gravitational collapse, we always encounter with two types of singularities 
— shell focusing singularity and shell crossing singularity. In Tolman-Bondi-Lemaitre 
(TBL) dust model, these two kinds of singularities will corresponds to i? = 0 and R' = 0 
respectively. A shell focusing singularity (i.e., i? = 0) on a shell of dust occurs when it 
collapses at or starts expanding from the centre of matter distribution. The instant at 
which a shell at the radial co-ordinate r will reach the centre of the matter distribution 
should be a function of r. So different shells of dust arrive at the centre at different 
times and there is always a possibility that any two shells of dust cross each other at 
a hnite radius in course of their collapse. In this situation the comoving system breaks 
down, both the matter density and kretchman scalar diverge [Seifert, 1979; Hellaby et 
al, 1984, 1985] and one encounters the shell crossing singularity. If one treats it as the 
boundary surface then the region beyond it is unacceptable since it has negative density. 
It is therefore, of interest to find conditions which guarantee that no shell crossing will 
occur. 

Goncalves [2001] studied the occurrence of shell crossing in spherical weakly charged 
dust collapse in the presence of a non-vanishing cosmological constant. The positive 
cosmological constant model conceively prevent the occurrence of shell crossing thereby 



CHAPTEB\1‘A Shell Crossing Singularities in Quasi-Spherical Szekeres Models 


157 


allowing at least in principle for a singularity free ‘bounce’ model. Nolan [1999, 2003] 
derive global weak solutions of Einstein’s equations for spherically symmetric dust-filled 
space-times which admit shell crossing singularities. Recently, Hellaby et al [2002] in¬ 
vestigate the anisotropic generalization of the wormhole topology in the Szekeres model. 
We have studied the shell crossing singularity in Szekeres model of the space-time both 
from physical and geometrical point of view. We study the physical conditions and ge¬ 
ometrical features of shell crossing singularities. 

12.2 Physical Conditions for Shell Crossing Singularity 

Here we consider (n -|- 2) dimensional Szekeres space-time model with metric ansatz is 
given in equation (10.1). Under the assumption that (3' ^ 0, the explicit form of the 
metric coefficients are also given in equations (10.2) and (10.3), where the form of v is 
given in equations (10.16) and (10.17) and R satished the differential equation (11.1) for 
non-zero cosmological constant A. 

We shall now make a comparative study of shell focusing and shell crossing singularity 
time and hnd conditions in favour (or against) of formation of shell crossing singularity 
for the following different choices: 


(i) /(r) = 0, A = 0: 


In this case equation (11.1) can be integrated to give 



( 12 . 1 ) 


So from (11.9), R(tsf{r),r) = 0 gives the result 



( 12 . 2 ) 


Now to avoid the shell crossing singularity either all shells will collapse at the same 
time (i.e., U/(r) is independent of r) or larger shell will collapse at late time (i.e., U/(r) 
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is a monotone increasing function of r). These two conditions can be combined as 


t'.f(r) > 0 


or equivalently from equation (12.2) 

F'{r) ^ n + 1 
F{r) ~ r 

Now combining equations (11.10), (12.1) and (12.2) we have 


(12.3) 


- t,,(r) = 


9r^ I 


(12.4) 


(n + l) v'F(r) + (n + ly 

But if it is so happen that B! + RF = 0 is a regular extremum for /3, then we must 
have hnite p. This implies from equation (11.2) that F' + (n + 1)FF = 0. Hence, if 
there is no shell crossing singularity corresponding to equation (11.10) we must have two 
possibilities: 


F'ir) 

either (a) —^ + (n + 1)F = 0 
F[r) 


and 


n + 1 F'{r) 

r F{r) 


= 0 


(12.5) 


,, , F'(r) . . , , n + 1 F'(r) . . 

or (b) + (n + l)i/ =0 and- > 0 (12.6) 

F{r) r F{r) 

For the hrst choice tsf{r) is constant, so all shells collapse simultaneously while for 

the second choice tsf{r) is a monotonic increasing function of r and there is an inhnite 

time difference between the occurrence of both type of singularities. 

The value of i? at f = tsdr) is 


n+l 

{fi((,c,r)} " 



(12.7) 


Therefore as an complementary event, the conditions for occurrence of shell crossing 
singularity are R Rif = 0, p = cxd, 0 , > 0, t;/(r) < 0. 

As Ff{r) < 0 implies 


F'ir) n + l 

—— > - 

F(r) 


r 


( 12 . 8 ) 
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so -R > 0 demands 

^ + («+l)>''>0 (12.9) 

Hence we have 

F(r) ~ and ~ (12.10) 

for shell crossing singularity with / > (n + 1) and (n + l)p > —/. 

(ii) /(r) = 0, A ^ 0: 

This choice will give the solution to equation (11.1) as 


t + 


2n 


{n + 1)A 


Sink ^ 


2Ar’*+i 


n(n + l)F(r) 


— Sink 


-1 


2Ai?’"+i 
n{n + l)F(r) 


( 12 . 11 ) 


At the shell focusing time ts/(r), R = 0, hence we have 


tsf{r) =ti + 


2n 


{n + 1)A 


Sink 


-1 


2Ar^+^ 


n{n + l)F(r) 


( 12 . 12 ) 


Since to avoid the shell crossing singularity, ts/(r) should be an increasing (or a constant) 
function of time i.e., > 0, which implies the same condition (12.3) as in case (i) 

and is independent of A (whether zero or not). Further, substitution of equation (12.11) 
in equation (11.10) will give 


tsc{r) - tsf{r) = tank ^ 



(12.13) 


But if there is no shell crossing singularity corresponding to R' + RC = 0 (then it will 
correspond to an extremum of j3) then p must be hnite. This will be possible only when 
F' + {n + l)Fv' = 0. But from eq.(12.13) it is permissible only when t(j(r) = 0 i.e., 
tsf{r) is independent of r. 
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(iii) /(r) 7 ^ 0,A = 0, R(ti,r) = 0 {time symmetry): 

In this case explicit solution is possible only for five dimension (i.e., for n = 3) and 
the result as 

= (12.14) 


But the shell focusing condition R{tsf{r),r) = 

tsf{r) =ti + 


0 gives 


VW)' 


So t{f{r) > 0 will give 

F 4 

-F < - (12-15) 

F r 

Here the time difference between the two types of singularities is 


tsc{r) - Rf{r) = 

The r.h.s. of equation (12.16) always positive by the inequality (12.15). 




' + 1 


— + C 

2F' 


- 1 


(12.16) 


12.3 Geometrical Features of Shell Crossing Singularity 

Now we shall discuss the shell crossing singularity from geometrical point of view. In 
fact geometrically, a shell crossing singularity (if it exists) is the locus of zeros of the 
function R! + Rv'. Now writing explicitly the function R' + Rv' using the solution (10.16) 
for z/ we have 


R' + Rf' = e" 


{R!A - RA') - RB[) X, + {CR! - C'R) 


(12.17) 


2 = 1 2=1 

The discriminants of this quadratic function in Xj’s {i = 1,2, ...,?7,) with respect to 
these variables are respectively 


Ax,- = —4 


{R!A - RA'f xl + {R!A - RA') {R'Bk - RB'^) Xk 


k=i-\-l 


k=i+l 













CHAPTEB\1‘A Shell Crossing Singularities in Quasi-Spherical Szekeres Models 


161 


+{R'A - RA'){CB! - C'R) - ^ ^(R'Bi - RB[) 
and 


1=1 


, = (12.18) 


/S.Xn = {R!A — RA 


/\2 


R^ ( I - R!‘^ 


, i=l 


(12.19) 


Thus R! + Rv' will have the same sign for all XiS when Aa;„ < 0 (because then also 
Axi,..., Axj <0, i = 1, 2,..., n — 1 V Xj+i, ...,Xn). Hence R' + Rv' has the same sign for 
all XiS (i.e., no shell crossing) if and only if 


p/2 

— > Bf - M'C'= «i(r) (say) (12,20) 

i=l 

But if -^ = ipir) then Ax^ = 0 and so Ax* = 0 (i = 1, 2,..., n — 1) at just one value 
of Xn- At this value of R' + Rv' = 0 at one value of Xi {i = 1,2, ...,n — 1). So 
here shell crossing is a single point in the constant (t,r)-hypersurface (n dimensional). 
In other worlds, it is a curve in the f-constant {n + 1)-D hypersurface and a 2 surface in 
{n + 2)-D space-time. 

When ^ < ip{r) then Ax„ > 0 and from equation (12.17), R' + Rv' = 0 will represent 
the equation of a n-hypersphere in the n dimensional XiS plane. Further, for Ax„ > 0 
the two limiting values Xn at which Ax* (z = 1, 2,..., n — 1) changes sign are 


-{R!B^ - RB'J ± -4A'C')-i?'2 

2{R'A - RA') 


( 12 . 21 ) 


Then for every x„ such that Xm < x„ < Xn 2 , there are two values of x„_i (only one if 
Xn = Xn^ or Xnj) such that R! + Rv' = 0 and so on. These values of Xj’s {i = 1, 2,..., n —1) 
are 


^?'i o 


2{R'A - RA' 


-{R!Bi - RB'i) ± Ir^ - 4A'C" j - 


,k=l 


-4 {R'A - RA') 


n ^ n 


1/2 


( 12 . 22 ) 
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These values of Xj’s will lie on the hypersphere with centre 


and radius 


/ {R'Bi-RB[) 
\ 2{R'A-RA') 


{R!B^-RB'^) \ 
2{R'A - RA') ) 


2{R'A - RA') 


This hypersphere is different from the hypersphere with i/' = 0 i.e., 


(12.23) 


A'{r) ^ ^ B[{r)xi + C\r) = 0 (12.24) 

i=l i=l 

So the shell crossing set intersects with the surface of constant r and t along the line 
(curve) ^ = —z/'=const ant. 

Now for positive density we note that F' and R' + Rv' must have the same 
sign. We now consider the case where both are positive (when both are negative, we just 
reverse the inequalities). When both are zero then it can happen for a particular value of 
Xj’s {i = 1, 2,..., n) if ~ ~ ^ which can not hold for all time. This is possible 

for all Xi if F' = R' = F = 0. This implies that at some r^ F' = f = A' = C = B'^ = Q 
{i = 1, 2,..., n). Hence we choose 

F' R' 

- -— > —F and — > —F (12.25) 

{n+l)F R ^ ’ 

Also from the solution (10.16) we have 


_ / _ xf + Er=i B'i{r)xi + C'jr) 

"" A{r)E:=i^l + EtlB^{r)x, + C{r) 

Now writing in a quadratic equation in xi we have for real xi, 


So 


F^ +\2{A' + AF) Y^xk + Y,{B'^ + BkF) I < - 4A'C' 


V 


/2 


n 

= V - 4A'C" 


max 
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Hence from (12.25) we have 


F' 


> 


B? - 4H'C" , V 


(12.26) 


("+i)-p \| tr 

which implies F' > 0, V r. Now for R' + Ru' > 0, we shall study the three possible 
choices separately. 


(i) /(r) = 0, A = 0: 

Here the solution for R can be written as 

(t - a(r)) 

So as f —>• a, R^R' + R^v' —^F{r) a'{r) and as f —> cx), ■^ + z/' —>• + F. 

Hence for R' + Ru' > 0 we must have a' < 0 and B'j^ — 4A'C' . 


(ii) /(r) = 0, A ^ 0: 

The solution for R can be written as 


„ _ n(n + l)F(r) 
2A 


{n + 1)A 
2n 


(t-a(r)) 


In this case as t —>• a, i ?"2 R' + _R "2 F —>• —^JF{r) a'{r) and as t —oo, -^ + F —>• 
n{n+i) ■ Thus for R' + RF > 0 we must have a'{r) < 0 and 

^ a'(r) > ^/Y.UB?-^A'C' . 


(n+l)F 

F’ 

{n+l)F 


n(n+l) 


(hi) /(r) 0, A = 0, R{ti, r) = 0, u = 3: 

Here the solution for R is 


R^ = r^ 


Fir), ,2 

-Fit - *.) ■ 


The limiting value of ^ + F as f —> cxo will be ^ F F — F Hence for R' + RF > 0 
we should have g + ^ > ^/Eti B? - 4A'C" . 
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12.4 Discussion 

In this chapter, a details study of shell crossing singularity has been done for dust model 
with or without cosmological constant for Szekeres model of {n + 2)-D space-time. The 
physical conditions however do not depend on A (whether zero or not) and the form of 
the conditions are identical for the three cases presented there. For geometrical condi¬ 
tions the locus of shell crossing depends on the discriminant of the co-ordinate variables 
Xj’s (i = 1, 2,..., n). If both ^ and (^n+i)F greater than a/X] ~ 4A'C" then there 
will be no shell crossing singularity even if B! + Rv' = 0. Here p is hnite and R! + Rv' = 0 
will correspond to a real extrema for j3. On the other hand if ^ = a/X) ~ dA'C' 
then shell crossing singularity is a 2-surface in (n -|- 2)-dimensional space-time. For 
^ < a/X] ~ AA'C, the shell crossing set lie on a n-hypersphere and it intersects 
with constant (t, r) along the curve = —C =constant. 
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Short Discussions and Concluding Remarks 


The thesis consists of two parts - Part A and part B. Part A consists of seven chap¬ 
ters dealing with cosmological solntions both in Einstein gravity and in Brans-Dicke 
theory of gravity with Varying speed of light and also qnintessence problem has been 
discussed for these solutions. In the last chapter (chapter 7) of part A, validity of Cos¬ 
mic No-Hair Conjecture has been studied in Brane world scenario. Part B consists of 
Eve chapters dealing with spherical gravitational collapse in higher dimensional Tolman- 
Bondi-Lemaitre models and Quasi-spherical gravitational collapse in higher dimensional 
Szekeres models. 

In chapter 1, the cosmological implications of the theory of varying speed of light and 
a brief idea of Quintessence models have been discussed. 

In chapter 2, Brans-Dicke cosmology in an anisotropic Kantowski-Sachs space-time 
model have been studied in varying speed of light theory. The flatness problem has been 
solved, when Brans-Dicke scalar held 0 is constant or variable, specially for radiation 
dominated era both perturbatively and non-perturbatively. The other problems like 
lambda, quasi-lambda and quasi-flatness problems have been solved for varying speed 
of light and their asymptotic behaviour has been examined. 

Chapter 3 deals with Friedmann-Robertson-Walker space-time models of a perfect or 
causal viscous fluid distribution with a scalar held 0 having potential V (0) for isotropic 
space-time model, considering variation of the velocity of light. For perfect fluid distribu¬ 
tion and radiation dominated era, it has been shown that V is monotonically decreasing 
function of 0 and 0 increases with c. For viscous fluid distribution, V decreases expo¬ 
nentially with 0 when co-efficient of bulk viscosity rj ~ and V is monotonically 

decreasing function of 0 when f] p. Therefore, for perfect and causal viscous fluid dis- 
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tribution, the potential function V is monotonically decreasing with Brans-Dicke scalar 
held 0. 

Quintessence problem in anisotropic cosmological models for Brans-Dicke theory has 
been discussed in chapter 4. It has been shown that it is possible to have an accelerated 
expanding universe without any matter due to Brans-Dicke held. In order to solve the 
held equations, it has been assumed power-law form of scale factors. It has also been 
shown that the anisotropy character is responsible for getting an accelerated model of 
the universe. The hatness problem has been discussed in this context using conformal 
transformation in Jordan frame. 

Chapter 5 deals with anisotropic cosmological models for Brans-Dicke theory but with 
a scalar held which is self-interacting. Assuming the power-law form for the scale factors 
and brans-Dicke scalar held, it has been studied the accelerating or decelerating nature 
of the universe for various values of the parameters. For the present universe (i.e., for 
dust models) it has been shown in principle that perfect huid with barotropic equation 
of state can be considered to study the quintessence problem. In addition to the non¬ 
decelerating solution, hatness problem has been solved potentially. 

An extensive analysis of quintessence problem for Brans-Dicke theories with varying 
speed of light have been studied in chapter 6. Here also without any quintessence mat¬ 
ter, an accelerated form of universe is possible. For the power-law forms and exponential 
forms of the cosmological scale factors and scalar held, the velocity of light is in the same 
form and identihed the cases where the quintessence problem can be solved for some re¬ 
strictions on parameters involved. Along with providing a non-decelerating solution, it 
can solve hatness problem also. 

In chapter 7, Cosmic No-Hair Conjecture in Brane Scenarios have been studied. It has 
been examined that for perfect huid model and for general form of the energy momen- 
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turn tensor for any matter fields (scalar field, perfect fluids, kinematic gases, dissipative 
fluids, etc.) whether Cosmic No-Hair Conjecture in brane scenario is satished if it is 
valid in general relativity. 

In chapter 8, a brief review of spherical and non-spherical gravitational collapse for 
dust and perfect fluid models have been discussed. 

Chapter 9 deals with Tolman-Bondi-Lemaitre dust collapse model in (n -|- 2) dimen¬ 
sional space-time. The local nature of the central shell-focusing singularity by a compar¬ 
ative study of the time of the formation of trapped surface and the time of formation of 
the central shell-focusing singularity has been examined. When the initial density gra¬ 
dient falls off rapidly and vanishes at r = 0, then naked singularity formation is possible 
only up to space-time dimension hve for marginally bound case and any dimension in 
non-marginally bound case. If the above restriction on the initial density distribution 
has been relaxed then naked singularity is possible in any dimension both for marginally 
and non-marginally bound cases. The geodesic study also supports the above results for 
marginally bound case. Also the numerical results for geodesic study in non-marginally 
bound case favour the formation of naked singularity in any dimension. 

Cosmological solutions for (n -|- 2)-dimensional Szekeres form of metric with perfect 
fluid (or dust) as the matter distribution and cosmological constant have been studied 
in chapter 10. The solutions have been found in two cases namely, quasi-spherical and 
quasi-cylindrical. For specihc values of the arbitrary functions (namely, A{r) = C{r) = ^ 
and Bi{r) = 0, V f) the Szekeres metric reduces to spherically symmetric metric for a 
suitable transformation of co-ordinates. The study of asymptotic behaviour shows that 
some of the solutions will become isotropic at late time while there are solutions for 
which shear scalar remains constant throughout the evolution. Thus Cosmic ‘no-hair’ 
Conjecture is not valid for all solutions. 
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In chapter 11, gravitational collapse in (n + 2)-diniensional space-time using a higher¬ 
dimensional generalization of the quasi-spherical Szekeres metrics with non-zero cos¬ 
mological constant have been studied. The local nature of the central shell-focusing 
singularity by a comparative study of the time of formation of trapped surface and the 
time of formation of the central shell-focusing singularity has been examined. When 
initial density gradient falls off rapidly and vanishes at r = 0, then naked singular¬ 
ity formation is possible only up to space-time dimension hve for > —1 and any 
dimensional space-time for = —1. If the above restriction on the initial density dis¬ 
tribution has been relaxed then the Cosmic Censorship Conjecture may be violated in 
any dimension. Also both by geodesic study and numerical study, it has been deduced 
that the nature of the central singularity depends sensitively on the choice of the initial 
data. The radial null and radial time-like geodesics and their local visibility have been 
studied. Therefore, it is possible to have at least local naked singularity for the given 
higher dimensional quasi-spherical space-time. 

Finally, in chapter 12, a details study to shell-crossing singularity has been done for 
dust model with or without cosmological constant for Szekeres model in (n -|- 2) dimen¬ 
sional space-time. The physical conditions do not depend on cosmological constant and 
for the geometrical conditions, the locus of shell-crossing depends on the discriminant 
of the co-ordinate variables. Also the conditions for no shell-crossing singularities have 
been studied. 
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